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Abstract 


A  study  has  been  made  of  convective  heat  transfer 
between  a  rotating  cylinder  and  a  stationary,  concentrlo 
outer  cylinder  for  the  Taylor  number  range  2000  to  9000,  and 
for  rotor  surface  temperatures  In  the  range  80  P  to  430  P. 

Measurements  were  made  of  overall  heat  transfer, 
and  velocity  and  temperature  distributions  In  the  gap  between 
the  two  cylinders.  A  hot  wire  anemometer  was  used  to  determine 
whether  secondary  flows  exist  in  the  gap  in  the  Taylor  number 
range  of  the  Investigation. 

The  principal  results  are  as  follows: 

1.  That  turbulent  flow  can  exist  In  the  gap  without 
an  axial  flow  of  air  In  the  range  of  Taylor  number  3000  to 
9000. 


2.  That  the  correlation  offered  by  Kays  and  BJorklund 
for  the  heat  transfer  in  the  gap,  Nu/Nucond>  ■  0.175  (Ta)1/2 
can  be  extended  for  low  temperature  differences  from  Taylor 
numbers  of  2000  to  Taylor  numbers  of  9000.  At  higher  temperature 
differences  the  heat  transfer  Is  reduced  but  still  varies  as 
(Ta)1/2. 


3.  That  V^+  vs.  y+  follows  the  universal  curve  for 
pipe  flow  out  to  y+  ■  10,  and  levels  off. 

4.  That  the  friction  factor  follows  a  relation 

„  0.3342 

1 rotor  “  0.358 

Ta 


5.  a.  That  near  the  walls  and  out  to  y+  ■  JjO, 


two  curves,  one  for  the  rotor,  and  one  for  the  stator,  suffice 
to  describe  €  ^/z'f  vs,  y+  ^or  3000  ^  Ta  <  9000,  and  for 
temperature  differences  of  from  60  P  to  330  P. 

b.  That  the  eddy  diffusivity  In  the  middle  of  the 
channel  can  be  expressed  by  the  relation 


(1  +  <-/*)  •  0.0894  (frotor)^  1 - f  (i+Jstatop)Ta 

"  r  a/rrotor  T  rotor 

1.75  -*r"<  1.95 


,3/2 


6.  That  the  ratio  &  /€  equals  1.4  for  the  conditions 

H  M 


studied. 


7.  That  the  correlations  obtained  can  be  used  to 
determine  overall  friction  factors  and  heat  transfer  for  the 
apparatus  In  the  ranges  of  Taylor  number  and  temperature  dif¬ 
ferences  studied  when  given  the  temperature  of  the  rotor, 
the  temperature  of  the  stator,  and  the  rpm  of  the  rotor. 
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Nomenclature 

area  of  surface,  ft2;  *  rotor*  °*  rotor,  ft2 

specific  heat,  Btu/lb . -deg .  7 
gap  width,  ft. 

electromotive  force  between  metals  a  and  b 
friction  factor,  dimensionless 
defined  In  text,  equation  (9*7) 
temperature  gradient  In  quarts  windows,  deg/lnoh 
stagnation  enthalpy  Btu/lb  • 

unit  convective  conductance  on  surfaoe  of  heated 

O 

cylinder,  Btu/hr-Ft  -deg  7 

A 

thermal  conductivity,  Btu/^ir-f t  -deg  7/ft 
Oladstone-Dale  constant,  ft^/lb 
length  of  rotor,  ft 
mixing  length 

magnification  factor  on  lnterferometrlo  plates 
index  of  refraotlon 
total  pressure,  lbs/ft2 

static  pressure  lbs/ft2 
heat  transferred  by  convection,  Btiyfor 
heat  transferred  at  radius  r,  Btu/hr-f t^, 
in  r  direction 

radius,  ft;  rrotor,  rotor  radius;  r8tator,  stator 
radius;  rc,  radius  at  oenter  of  gap 
temperature,  deg  R 

mean  velocity  of  flow  in  x  direction 
velocity,  ft/seo 


viii 


v# 

, rotor 

V*,c 

Vr 

Vz 

y 

n 

yt,m 


velocity  in  ^  -direction 

velocity  of  rotor  surface 

velocity  at  center  of  gap 

velocity  In  r-dlrectlon 

velocity  In  z -direct Ion 

distance  from  wall,  ft 

distance  from  wall,  thousandths  of  Inches 

distance  from  wall,  as  Indicated  by  total  pressure 

probe.  V.  displacement  correction 


Greek  letter  symbols 


A 

4 

Fringe 

6 

<5*7 

* 

X 

A|4u 

M 

f 


ratio  between  eddy  dlffuslvlty  for  heat  and  eddy 
dlffuslvlty  for  momentum 
coefficient  of  thermal  expansion  for  quartz 
windows  inches/inch -deg  F 
fringe  shift.  Inches 
spacing  of  no -flow  fringes 

fringe  shift  measured  In  wavelengths  of  light 
eddy  dlffuslvlty  for  heat,  ft2/sec 
eddy  dlffuslvlty  for  momentum,  ft2/sec 
thermal  dlffuslvlty, 
wavelength  of  light.  Angstrom  units,  A 
wavelength  of  light  In  vacuum.  Angstrom  units 
viscosity,  (lbs)(sec^ft2 
kinematic  viscosity,  ft2/sec 
density,  lbs/ft^  or  slugs/ft^ 


>* 


lx 
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T '  shear  stress,  lbs/ft2 

Cd  angular  veloelty  radlans/sec;  ^  erltleal 

angular  velocity  at  which  Taylor  vortices  develop 


Dimensionless  Numbers 


^ft)6or 
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S&eti&r 
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Nu 


cond 


Pr 

Re 

Ta 
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friction  factor  on  rotor  surface, 

defined  by  V"  ,  m  ~  'fcotor — 

3  *n>6*r  0  ,/  z 

Sfbhor  V*>C 

friction  factor  on  stator  surface, 

defined  by  fMaiar  =  —^7~i 

/s&tor  fyc 

Nusselt  Number, 

Conduction  value  of  Nasselt  number  In  the  gap, 
defined  by 


Prandtl  numb e r 

Taylor  munber, 


Reynolds  number,  Z/d/v  ; 


friction  velocity, 


non-dlmenslonallzed  velocity, 
non -dlmens Iona llzed  distance  from  the  wall,  /<U 


Unless  otherwise  stated,  properties  In  dimenslonleaa 
numbers  are  evaluated  at  the  average  between  the  rotor  and 
stator  temperatures . 

The  barred  and  primed  quantities  refer  to  time  mea^ 
and  fluctuating  quantities  in  turbulent  flow  and  are  defined 
by  equations  (2-9)  to  (2-12)  and  (2-14)  to  (2-16). 


1 


Chapter  I 
Introduction 

1.1  Purpose  of  Investigation 

The  purpose  of  this  study  Is  to  determine  experi¬ 
mentally  velocity  and  temperature  distributions,  and  overall 
heat  transfer  In  the  annulus  between  a  rotating,  hot  cylinder, 
and  a  concentric,  stationary  outer  cylinder  under  conditions 
such  that  the  properties  of  the  separating  air  are  markedly 
variable.  Emphasis  is  placed  on  the  turbulent  regime  of  flow. 

These  data  are  useful  since  they  enable  the  prediction 
of  overall  heat  transfer  rates  In  turbulent  flow  from  a  knowledge 
of  the  velocity  field;  this  through  the  concept  of  the  eddy 
dlffuslvltles  for  momentum  and  heat  transfer. 

The  concentric  cylinder  apparatus  was  chosen  for 
this  work  because  it  provides  a  means  for  obtaining  high  rates 
of  energy  transfer,  and  relatively  high  temperature  differences 
In  a  small  apparatus.  Further,  in  the  absence  of  secondary 
flow  effects  this  apparatus  Is  capable  of  producing  one  of 
the  few  possible  examples  of  established,  compressible  flow. 
Analysis  of  the  data  is  then  simple  in  certain  respects. 

Other  problems  are  Introduced,  however,  such  as  the 
existence  of  Taylor  Vortices,  a  three-dimensional  secondary 
flow,  which.  If  persisting  at  ranges  of  speed  covered  by  the 
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present  experiment  would  complicate  the  analysis. 

1.2  Summary  of  Previous  Work: 

1.2.1  Flow  Between  Rotating  Cylinders 

A  discussion  of  previous  investigations  of  flow 
between  rotating  cylinders  can  be  conveniently  divided  into 
three  parts;  one  concerned  with  the  flow  regimes  within  the 
annulus;  the  second  with  the  heat  transfer;  and  the  third 
with  velocity  and  temperature  measurements  in  the  annulus. 

1.2.1. 1  Flow  Regimes 

The  breakdown  to  turbulent  flow  of  the  curved  Couette 
flow  existing  between  the  cylinders  at  low  speeds  of  rotation 
of  the  inner  cylinder  does  not  occur  in  the  same  manner  as  in 
pipes  or  on  flat  plates.  In  these  the  laminar  flow,  at  a 
critical  Reynolds  Number  Uk>  ,  becomes  unstable  to  disturb¬ 
ances  existing  in  the  flow.  These  disturbances  increase  with 
time  and  change  the  character  of  the  flow  from  the  well- 
behaved  laminar  to  the  random  fluctuations  of  turbulent  flow. 

The  curved  Couette  flow  on  the  other  hand  responds 
to  small  disturbances  in  a  different  manner.  At  a  critical 
value  of  the  parameter  called  the  Taylor  number, 

the  laminar  flow  breaks  down.  This  reakdown  is  not  into  the 
random  fluctuations  which  typify  turbulent  flow  but  into  a 
new  type  of  laminar  flow  containing  torroidal  cells  of  secondary 
motion  looping  the  inner  cylinder,  as  shown  in  Fig.  1-1. 
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*/rs  of  (fe>ur>fsr- 

TayJor-  U?r/tces 


In  a  classical  paper  0.  I.  Taylor1  predicted  not 
only  the  development  of  these  vortices  but  also  the  conditions 
of  cylinder  rotational  speeds  and  fluid  properties  at  which 
they  would  ooour.  Taylor  superimposed  on  the  velocity  compon¬ 
ents  for  laminar  flow  small  s/nusoldal  disturbances. 


Substitution  of  the  resultant  velocities  Into  the 
conservation  equations  and  the  solution  of  these  yielded,  for 
an  Incompressible  fluid,  and  for  t  a  critical  speed 

of  rotation  for  the  Inner  cylinder,  beyond  which  the  pertur¬ 
bations  would  grow. 

The  vortices  formed  are  called  Taylor  Vortices. 

This  flow  Is  laminar  In  nature  and  represents  a  secondary 
flow  superimposed  on  the  circumferential  flow. 


For  the  case  where  the  outer  cylinder  Is  stationary 
Taylor  found  that  the  laminar  vortex  motion  sets  in  at  a 


critical  speed  given  by 

*  - 


(1-D 


where 


For  the  case 


and 


P  «  a  057 t  (/-  o  CS10-  )  +  — 'n 


o}  p* 


0.0577 


.VV  —  /&37 
~4/.Z 


JS*  On 


(1-2) 


"Stability  of  a  Viscous  Fluid  Contained  between  Two  Rotating 
Cylinders,"  0.  I.  Taylor,  Philosophical  Trans,  of  the  Royal 
Society  of  London,  series  A,  vol.  223,  1923*  PP»  289-3^3 • 
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Taylor  alao  predicted  the  vortex  spaolng,  which  he  found  to  be 
equal  to  the  annulus  width.  Experiments  conducted  by  him, 
using  water  as  the  fluid,  confirmed  the  analytical  value  of 
to  within  2$.  Further,  the  vortex  spacing  measured  by  Taylor 
was  1.05  cm,  for  an  annulus  width  of  1.105  cm. 

J.  W.  Lewis1  also  investigated  experimentally  this 
type  of  motion  and  found  that  the  vortices  themselves  mere 
stable  to  arbitrary  small  disturbances,  such  as  changes  In 
the  speed,  much  as  In  pipe  flow  where  a  turbulent  flow  persists 
at  Reynolds  Numbers  below  the  critical  as  the  speed  Is  reduced. 
He  further  found  that  for  -  0.17,  thelamlnar  vortex  flow 

became  "Irregularly  turbulent"  at  a  speed  twice  the  critical 
value.  He  did  not  observe  this  breakdown  at  speeds  up  to 
five  times  the  critical  for  the  values  -0.43,  and 

°.71.  .  ’  /  ’•>  "• 

The  turbulent  flow  which  ensues  at  speeds  much  higher 

2 

than  the  critical  has  not  yet  been  analysed.  Pal  Shlh-I, 
however,  has  Investigated  experimentally  the  case  for  air  flow 
between  the  cylinders  at  Taylor  numbers  from  10-20,000.  He 
concluded  from  hot  wire  velocity  measurements  and  pressure 
distributions  that  secondary  vortices  persist  even  at  speeds 
several  hundred  times  beyond  the  critical.  He  also  determined 


lnAn  Experimental  Study  of  the  Notion  of  a  Vlsoous  Liquid 
Contained  Between  Two  Coaxial  Cylinders,"  Proo.  Soy.  Soo. 
London,  Ser.  A,  117,  pp.  388-407  (1927-1928). 

^"Turbulent  Flow  between  Rotating  Cylinders,"  Pal  Shlh-I, 
NACA  TN-892,  March  19*3. 
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that  the  vortices  occur  In  pairs  and  their  nuaber  decreases 
with  the  speed  of  the  inner  cylinder;  that  in  certain  areas 
the  initial  vortex  motion  breaks  down  into  smaller  eddlces 
which  degenerate  further  to  produce  regions  where  turbulent 
fluctuations  predominate.  Pal's  conception  of  the  flow 
pattern  is  shown  in  Fig.  1-2. 

Pal's  conclusions  were  based  on  two  experimental 
observations;  a  non-uniform  axial  pressuredistrlbutlon  and 
two  different  types  of  mean  velocity  distribution  observed 
at  different  axial  locations. 

He  further  observed  at  high  speeds  sudden  changes 
in  the  velocity  distribution  at  a  given  axial  location  from 
one  type  to  another,  typified  by  a  sudden  drop  in  mean  speed. 

Pai  associated  each  of  the  two  types  of  velocity 
distributions  with  a  region  in  the  flow;  one  type  in  those 
regions  where  the  secondary  flow  predominates,  the  other 
type,  the  region  In  which  the  random  turbulent  fluctuation 
predominate.  Thelatter  velocity  distribution,  he  found  was 
of  a  logarithmic  type* 

The  findings  of  Pal  were  extended  by  Kaye  and  Elgar^ 
who  studied  the  turbulent  vortical  flow  with  smoke  tracer  and 
hot-wire  techniques.  Their  experiments  were  performed  both 
with  and  without  axial  air  flow  through  the  annulus.  They 


^ "Modes  of  Adiabatic  and  Diabatic  Fluid  Flow  in  an  Annulus 
with  an  Inner  Rotating  Cylinder,"  J.  Kaye  and  E.  C.  Elgar, 
ASME  Paper  #57-HT-l4,  1957. 


D/skrki c/  Jay/or  Vortices 
/Jy/>o~/heS/2e<V  6y  &/ 


ffyisne  /  -  £ 
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suggest  that  four  nodes  of  flow  exist  In  tho  annulus,  depending 
on  tho  axial  flow  Reynolds  Number  and  tho  Taylor  Number.  Fig. 

1-3  shows  a  scheaatlo  ropresontatlon  of  their  results.  In 
whloh  the  four  regime s,  lanlnar,  laalnar-plus-vortloes,  turbulent, 
and  turbulent -plus -vortloes,  are  delineated.  Included  also 

as  dashed  lines,  are  the  separation  lool  conjeoted  by  Booker 

1  2 

In  his  dissertation.  ' 

Becker's  work  will  be  discussed  later  In  connection 
with  the  heat  transfer  results  of  the  present  experiment. 

Included  In  his  dissertation  Is  an  excellent  review  of  the 
literature  for  the  flow  under  discussion  together  with  an 
extensive  blbllograph  for  both  the  dlabatle  and  adiabatic 
cases. 


^"An  Experimental  and  Theoretical  Study  of  Heat  Transfer  In 
an  Annulus  with  an  Inner  Rotating  Cylinder,*1  So.  D.  Thesis, 

Kurt  Becker,  Massachusetts  Institute  of  Technology,  Aug.  1957* 

2 In  a  subsequent  paper  Becker  and  Kaye  suggest  a  different 
arrangement  of  the  regimes  for  the  diabetic  and  adiabatic  oase 
as  concerns  the  transition  from  laminar  flow  to  turbulent  flow 
and  laminar-plus -vortices  flow  to  turbulent -plus -vortices  flow. 
This  is  apparently  a  reevaluation  of  the  work  of  Kaye  and  Elgar 
In  this  region  since  Becker  and  Kaye  state  "Based  on  these 
[temperature]  profiles  [taken  by  Kaye  and  Becker  In  each  of 
these  regions]  and  upon  the  difficulties  of  Interpretation  of 
hot-wire  data  for  the  adiabatic  oase,  we  reoonmend  the  use  of 
a  nearly  horizontal  demarcation  line,  between  the  laminar  and 
turbulent  regimes,  ...  for  both  adlabatlo  and  dlabatlo  flows." 
See  "K.  M.  Becker,  J.  Kaye,  "Measurements  of  Diabetic  Flow  In 
an  Annulus  with  an  Inner  Rotating  Cylinder."  ASMS  paper  Mo. 
61-SA-19,  August,  1961 


Axial  Flow  Reynold's 
Huaber 
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Figure  1-3 

Soheaatio  Representation  of  Differ out  Regions 
of  Flow  in  the  Annulus  (dashed  line:  Beoker1!  Hypothesis) 
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To  summarize  this  section,  one  can  say: 

1.  The  breakdown  of  the  curved  Couette  notion 
under  dlsousalon  does  not  proceed  directly  from  laminar  to 
turbulent  flow  but  rather  from  laminar  to  a  vortical  flow, 
still  laminar,  superimposed  on  the  mean  circumferential  flow. 
This  transition  ocours  at  a  critical  value  of  the  Taylor  number. 

2.  At  speeds  much  higher  than  the  critical  these 
vortices  still  persist. 

3.  Indications  are  that  at  still  higher  Taylor 
Numbers  the  vortices  break  up  completely  and  a  random  turbulent 
motion  occurs. 

1.2. 1.2  Overall  Heat  Transfer 

Most  prior  studies  of  heat  transfer  from  a  rotating 
Inner  cylinder  to  a  concentric  outer  cylinder  have  been  limited 
to  a  determination  of  overall  heat  transfer  coefficients.  The 
earlier  studies  are  reported  by  Kaye  and  Elgar, ^  Becker,1* 

Qazley2  and  Kays  and  Bjorklund.3  These  latter  three  works  are 
the  most  extensive  to  date.  Qazley  measured  the  heat  transfer 


1Kaye  and  Elgar,  ASME  paper  #57-HT-l4. 

4 

Becker,  Sc.D.,  Thesis,  MIT,  1957* 

2 "Heat  Transfer  Characteristics  of  Rotational  and  Axial  Plow 
between  Concentric  Cylinders,"  Carl  Qazley,  Jr.,  A.S.M.E. 
Paper  #56-A-128,  1956. 

^"Heat  Transfer  between  Concentric  Rotating  Cylinders,"  I.  S. 
Bjorklund  and  V.  M.  Kays,  ASMS  Paper  #58-A-99,  1958* 
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in  the  gap  for  smooth  and  Blotted  rotors  and  statore  with  and 
without  axial  through  flow.  The  length  of  the  annulus  was 
10  Inches ^  and  air  gap  diameter  was  about  5  Inches.  Different 
rotors  were  used,  with  air  gaps  of  either  0.017'  Inohes  or 
0.24  Inches.  Rotational  speeds  varied  from  540  to  4700  rpn 
with  temperature  differences  of  about  100  to  150  deg  P.1 
The  Taylor  Number  range  was  about  20  to  80  for  ^r»4r  ■  0.0068 
and  from  about  250  to  3000  for  <^//p^r  -  0.095* 

Becker  performed  similar  studies  with  a  smooth  rotor, 
and  was  mainly  concerned  with  obtaining  the  lines  of  demarcation 
for  the  four  flow  regimes  mentioned  earlier. 

The  experiments  of  Kays  and  Bjorklund  involved  smooth 
surfaces  without  axial  airflow  but  with  provision  for  rotation 
of  the  stator  as  well  as  the  rotor. 


The  results  of  Kays  and  BJorklund,  and  Becker 
(considering  only  his  results  with  no  axial  air  flow)  for 
la  ^  90  can  be  placed  on  a  line  (/ 

-  0.175  7*  (1-3) 


where 


(1-4) 


and  Is  based  on  rrotor,  and  fluid  properties  evaluated 
at  a  temperature  halfway  between  the  temperatures  of  the  two 
cylinders . 


Personal  communication,  C.  Oazley,  Jr.  to  0.  S.  Longobardo, 
June  24,  1939. 
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The  data  of  Oazley  follow  the  sane  trend  but  are 
higher,  possibly  beoause,  as  Becker  suggests,  he  did  not  take 
into  aooount  radiation  between  the  two  cylinders. 

The  physical  dimensions  and  ranges  of  temperature 
difference  between  the  rotor  and  stator,  and  for  the  three 
works  cited  and  the  present  work  are  given  in  the  following 
Table : 


Table  1.1 


Investigator 

X 

A^degP 

dL 

' /rotor 

46 

Length 

L,  inches 

Oazley 

20-80 

100-150 

0.0068 

0.0017 

10 

250-3000 

0.095 

0.024 

Becker 

0-575 

30-62 

0.233 

0.0058 

56 

Kays  &  BJorklund 

0-2000 

60 

0.054- 

0.0067- 

18.25 

0.246 

0.0305 

Longobardo 

3000-7000 

60-350 

0.227 

0.084 

4.5 

As  shown 

later,  the 

overall 

heat  transfer  results 

for  low  temperature  differences  obtained  in  the  present  experi¬ 
ment,  although  at  higher  Taylor  Numbers  can  be  placed  on  the 
line  represented  by  eq.  (1-3 )• 


1,2. 1.3  Velocity  and  Temperature  Distribution  Data 


Detailed  time-mean  velocity  distributions  for  purely 
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rotational  flow  have  been  obtained  by  Taylor,1  Pai,^  and 
Wattendorf .3 

Taylor's  experiments  were  conducted  to  differentiate 
between  his  vorticity  transfer  theory  and  the  momentum  transfer 
theory  of  Prandtl .  The  results  for  only  one  rotational  speed 
are  reported  in  detail,  that  corresponding  to  a  Taylor  Number 
of  60,000.  Taylor 's^aW  pitot  tube  data  gave  the  result  that 
Increased  toward  the  stator  in  the  middle  of  the  channel. 
Taylor,  noting  the  small  magnitude  of  the  variation  (about  6$) 
tried  to  determine  how  much  of  the  variation  was  contributed 
by  pitot  tube.  He  then  reran  the  experiment  using  a  larger 
pitot  tube  (by  a  factor  of  2  on  outside  diameter)  and  obtained 
a  greater  variation  in  \J^K . 

Using  these  data  to  extrapolate  to  a  pitot  tube  of 
zero  diameter,  he  obtained  the  result  -  constant  in  the 
central  85$  of  gap  width.  Unfortunately,  Taylor's  data  do  not 
approach  the  wall  sufficiently  close  for  further  evaluation. 


^'Distribution  of  Velocity  and  Temperature  between  Concentric 
Rotating  Cylinders,"  G.  I.  Taylor,  Proc.  Royal  Society  of 
London,  ser.  A,  vol.  157,  no.  892,  Dec.  2,  1936,  pp.  494-512. 

2Pai  Shih-I,  NACA  TN  892. 

3 

"A  Study  of  the  Effect  of  Curvature  on  Fully  Developed 
Turbulent  Flow,"  F.  L.  Wattendorf ,  Proc.  Royal  Soolety  of 
London,  ser.  A,  no.  865,  vol.  148,  Feb.  1935*  PP*  565-598* 
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Pal  Shlh-I'a  experiments  were  conducted  to  Investigate 
the  validity  of  Taylor's  extrapolation,  and  have  already  been 
discussed  In  Section  1.2 .1.1.  Again,  no  data  were  obtained 
close  to  the  wall,  but  Pal  found  that  away  from  the  wall  the 
velocity  obeyed  a  logarithmic  law. 

1.2.2  Determination  of  Eddy  Dlffuslvltles  for  Momentum  and 
Heat 


Consider  a  laminar.  Incompressible  flow  with  a  mean 
velocity  £7  in  the  X  direction.  One  writes  for  the  shear 
stress 


T/f~ 


In  a  turbulent  flow  an  augmentation  of  the  transfer 
of  momentum  occurs  and  It  is  customary  to  write 


?/? ,  %)<& 


(i-M 


where  is  the  "eddy  diffusivity  for  momentum."  Similarly, 
for  the  heat -transfer  one  can  write  for  the  laminar  case 


dl 


(1-5) 


and  for  the  turbulent  case 


/V  c  Z/fc)  d<f 

where  is  the  "eddy  diffusivity  for  heat." 

Defining  ot  as 

<*•»  *«/<&, 


(1-6) 


and 


noting  that^?*- 


equation  (1-6)  can  be  written 


(1-7) 
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„  jL.  dl 

/*?  d<f 
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(1-5) 


(1-6) 


Defining  o(  as 

a  <•*/{&, 
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(1-7) 


15 


^  s' t /  '*l  dr  .  o » 

/*?  L  *scy  (1_8) 

The  value  of  of  now  takes  on  an  added  significance 
since  with  a  knowledge  of  of  and  the  velocity  field  the 
temperature  field  and  hence  the  heat  transfer  can  be  determined. 


Mechanistic  theories  of  turbulence  such  as  Prandtl's 

l  2 

Mixing  Length  Theory,  and  Taylor’s  Vorticlty  Transport  Theory 

have  been  offer  ed  to  arrive  at  a  useful  formulation  for 

and  subsequently, 


Prandtl,  on  the  assumption  that  a  "lump”  of  fluid 
In  moving  from  a  position  ^  in  the  fluid  to  a  position 
conserves  not  only  its  identity  but  its  momentum  as  well 
arrives  at 

€m 


where  ^  is  called  the  mixing  length  and  Is  analagous  to  the 
mean,  free  path  in  kinetic  theory. 

The  theory  has  its  limitations;^  one  such  being  that 


■^L.  Prandtl,  Essentials  of  Pluld  Dynamics  (Hafner  Publishing 
Company,  Inc.,  New  York,  195*0,  P •  Ilf  bo  121. 

2S.  Goldstein,  Modern  Developments  In  Fluid-  Mechanics  (Oxford 
at  the  Clarendon  Press,  193*0  vol.  I,  pp.  209-214. 

^R.  W.  Stewart,  "A  New  Look  at  the  Reynolds  Stresses," 
Canadian  Journal  of  Physics,  vol.  3^*  1956,  pp.  722-725. 
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it  requires  to  go  to  zero  at  points  of  maximum  velocity, 
suoh  as  at  the  center  of  a  pipe.  This  Is  not  borne  out  experi¬ 
mentally.  Another  Is  that  values  of  £  obtained  from  experi¬ 
mental  data  are  large  In  a  given  flow,  contrary  to  what  has 
been  assumed  In  the  theory. 

Hlnze3-  has  calculated  the  value  of  the  mixing  length 
from  measured  velocity  distributions  In  the  wake  of  a  circular 
cylinder.  He  found  that  Is  Infinite  at  the  axis  of  the 

wake,  zero  at  the  boundary,  and  obtained  an  average  value  of 
equal  to  1/5  to  1/4  the  wake  width. 


If,  using  this  theory.  It  is  assumed  that  heat  Is 
transferred  In  the  same  manner  as  momentum  It  can  be  shown 


that 


and 


The  assumption  p/  ■  1  has  been  used  by  Deissler, ^>3, 4, Ly0n5 


1J.  0.  Hlnze,  Turbulence  (New  York:  McOraw-Hlll  Book  Co.,  Inc., 
1958),  p.  394. 

2R.  G.  Deissler,  and  C.  S.  Elan,  "Analytical  and  Experimental 
Investigation  of  Fully  Developed  Turbulent  Flow  of  Air  In  a 
Smooth  Tube  with  Heat  Transfer  with  Variable  Fluid  Properties," 
NACA  TN  2629,  1952. 

3r.  G.  Deissler,  "Analysis  of  Turbulent  Heat  Transfer,  Hass 
Transfer  and  Friction  In  Smooth  Tubes  at  High  Prandtl  and 
Schmidt  Numbers,"  NACA  TN  3145,  1954. 

^R.  G.  Deissler  and  M.  F.  Taylor,  "Analysis  of  Turbulent  Plow 
and  Heat  Transfer  In  Non-Circular  Passages,"  NACA  TN  4384,  1958. 

^R.  N.  Lyon,  "Liquid  Metal  Heat  Transfer  Coefficients,"  Chemical 
Engineering  Progress,  vol.  47,  1951*  PP»  75-79* 
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and  Seban  and  Shlmlzakl,1  for  example ,  to  oaloulate  the  temper¬ 
ature  distribution  In  a  fluid  from  a  known  veloolty  distribution. 

However,  the  assumption  ^  «  1  la  not  borne  out  by 
experimental  evidence.2  In  a  pioneering  experiment  Isakoff 
and  Drew3  found  that  0^  was  a  function  of  Reynolds  Number  and 
radial  position  In  the  tube.  Further  there  seemed  to  be  a 
trend  toward  unity  as  the  turbulence  Increased. 


R.  A.  Seban  and  T.  T.  Shimlzakl,  "Heat  Transfer  to  a  Fluid 
Flowing  Turbulent ly  In  a  Smooth  Pipe  with  Walls  at  Constant 
Temperature,"  Trans.  ASHE,  Vol.  73,  1951,  pp.  803-809. 

2In  "The  Principles  of  Turbulent  Heat  Transfer,"  NACA  TM  1408, 
1957*  H.  Relchardt  states,  "In  this  connection  [the  ratio]... 
a  far  reaching  congruence  was  obtained  of  the  temperature  and 
velocity  profiles  In  friction  layers  [from  the  results  of 
F.  Ellas  "Heat  Transfer  from  a  Heated  Plate  to  Flowing  Air," 
NACA  TM  6l4,  originally  published  In  1930].  This  result  was 
looked  upon  as  a  confirmation  of  the  previously  held  conception 
of  the  Identity  between  momentum  exchange  and  heat  transport 

This  conclusion  from  the  measurement  results  of  Ellas  was, 
however.  In  error,  as  was  shown  by  the  author  [Relchardt]  In 
the  year  1940."  Ellas'  reported  results,  «<  ■  1,  were  based 
on  m  1  whereas,  as  Relchardt  pointed  out,  /£>  ■  0.72, 
leading  to  a  result,  W  ■  1.4. 

3S.  E.  Isakoff  and  T.  B.  Drew,  "Heat  and  Momentum  Transfer  In 
Turbulent  Flow  of  Mercury,"  General  Discussion  on  Heat  Transfer, 
The  Institution  of  Mechanical  Engineers,  London,  1951,  pp.  405- 
409. 
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Similarly,  Brown,  Amstead,  and  Short, *  again  using 
mercury  as  the  fluid  found  **<  between  0.7  and  0.9  for  200,000 
Re  4  800,000  In  the  center  of  the  pipe,  dropping  rapidly 
near  the  wall. 


Seban  and  ShlmlzaklP  using  air  in  a  heated  pipe 
found  <2*  to  be  about  1.2  for  90^  300,  and  about  unity 

y.  j. 

for  y  >300,  where  is  the  non -dimens  Iona  lized  distance 
from  the  wall. 


Corcoran,  Page,  et.  al.3  studying  the  flow  of  air 

at  small  temperature  differences  in  a  rectangular  channel 

found  ot*  ijp  »1.4  near  the  wall,  and  decreasing  with 
'  'Qt/r 

increasing  Reynold's  number  in  the  middle  of  the  channel. 

li 

More  recently,  Sleicher,  using  air  in  a  heated  pipe 
found  to  be  greater  than  unity  to  vary  with  Reynold's  number 
and  position, and  approach  1.4  near  the  walls. 


None  of  these  experiments  was  made  under  conditions 
of  property  variation  of  the  fluid. 


*H.  E.  Brown,  B.  H.  Amstead,  and  B.  E.  Short,  "The  Transfer  of 
Heat  and  Momentum  in  a  Turbulent  Stream  of  Mercury,"  Trans. 
ASME,  vol *  79,  1957,  PP.  279-285. 

2  ' 

R.  A.  Seban  and  T.  T.  Shiraazaki,  "Temperature  Distribution  for 
Air  Plowing  Turbulently  in  a  Smooth  Heated  Pipe,"  General 
Discussion  on  Heat  Transfer,  The  Institution  of  Mechanical 
Engineers,  London,  England,  1951>  pp.  122-126. 

^W.  H.  Corcoran,  F.  Page,  Jr.,  W.  G.  Schlinger,  and  B.  H.  Sage, 
"Temperature  Gradients  in  Turbulent  Gas  Streams,"  Industrial 
and  Engineering  Chemistry,  vol.  44,  1952,  pp.  410-430. 

^C.  A.  Sleicher,  Jr.,  "Experimental  Velocity  and  Temperature 
Profiles  for  Air  in  Turbulent  Pipe  Plow,"  Trana.  ASME,  vol.  80, 
1958,  pp.  693-704. 
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Under  conditions  of  high  property  variation, 

DelBsler  and  Elan*  report  total  pressure  profiles  taken  for 
air  flowing  In  a  heated  tube  at  tube  wall  temperatures  up  to 
1500  F,  and  8000  ^  Re  ^  500,000,  and  temperature  profiles 
taken  up  to  wall  temperatures  of  only  300  F  because  of  the 
delicacy  of  the  thermocouple  probe  used.  Based  on  the  experi¬ 
mental  data  shown  for  ^  vs.  (their  Fig.  9)  and  £+ve.  if*" 
(their  Fig.  11)  the  scatter  of  the  data  appears  too  great  to 
permit  an  evaluation  of  of  .  Cert a/nty  they  report  no  data 
near  the  wall 

"...because  the  distortions  were  measured  at  high  Reynolds 
numbers  where  the  severe  velocity  gradients  and  the  presence 
of  the  hole  In  the  tube  wall  make  the  accuracy  of  the  measure¬ 
ment  doubtful"  (TN  2629  p.  17);  and  for  temperature  —  "The 
distributions  shown  for  the  region  close  to  the  wall  may  be 
subject  to  errors  due  to  conduction  along  the  thermocouple 
probe  prongs"  (TN  2629  p.  18.). 

Values  of  °<  have  been  obtained  from  experiments 
with  heated  jets  of  air  Issuing  Into  still  air. 

Hlnze  and  VanderHeggeZl jnen^  report  an  average 
value  &  1.36. 

Corrsin^  obtained  a  ratio  of  1.43  Independent  of 
the  magnitude  of  the  temperature  differential  up  to  300  C. 

These  results  are  for  turbulence  not  affected  by  the 
presence  of  a  solid  boundary. 

^R.  G.  Delssler  and  C.  S.  Elan,  NACA  TN  2629. 

2 

J.  0.  Hlnze,  Turbulence ,  p.  426. 


3s.  Corrsln,  M.  S.  Uberol,  NACA  TN  2124,  1950 
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Chapter  II 

Development  of  Basic  Equations1 
2.1  General  Equations  of  Motion 


The  transient  equations  of  linear  momentum  are  written 


as 


£ 


(2-1) 


where  body  forces  are  neglected  is  given  by 

It  -  is*  fv-  &J 

where  is  given  as  „ 

»•=*#/ 


4 


and 

'  '>  *'/, 

The  continuity  equation  is 

U*Jkc*o~ 


When  these  equations  are  transformed  to  cylindrical  coordinates , 
iff  2  )  in  terms  of  the  shear  stresses,  defined  below, there 
results 


*The  equations  in  this  chapter  are  to  a  great  extent  developed 
from  course  notes  in  fluid  mechanics  and  convective  heat  transfer 
taught  by  Professor  H.  0.  Elrod  at  Columbia  university. 
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r  -  component 


£>  -  component 
/&  + 


a> 


_L  / 

/*  ^ 


- 


; 


(2-2) 


2-  -  component 


/  -?S 


/- 


is  given  by 


(2-3) 


(2-1 ) 


4 


22 


The  Continuity  Equation  la 

0  &//"*) + /i 

(2-5) 


2.1.1  Equations  of  Motion  for  Laminar  Flow 


For  steady  state,  incompressible,  Couette-type 
laminar  flow  In  the  gap  between  the  inner,  rotating  cylinder 
and  the  outer  stationary  cylinder 9 


A  £_  £ 


The  Equation  of  Continuity  becomes 


At  r 


r  V 

rotor*  r 


0,  and  vr 


■  constant 
-  0  everywhere. 


(2-6) 


and  since 


The  -component  of  the  Equation  of  Motion  gives 


there  results 


or 


£  * 

■  constant 

/  Jr 


(2-7) 


constant 


^hus Constant,  and  ^  which  Is  a  sta tenant 

of  Conservation  of  Angular  Momentum. 
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Solving  for  V} 

with  Boundary  Conditions 

r  "  vvotov*  V#  ■  ^rotor  rrotor 
r  -  r.t.tor<  V  -  0 


There  results 


%otn  ) 

r  X  /-  '«■  yT  ^ 


(2-8) 

A  plot  of  Vp.  for  the  present  apparatus  at J7“raflians/sec 
(30rpa)  Is  shown  In  Pig.  2-1  .  Here,  V^rotor  -  0.43  ft/sec,  and 

7a  ■  36,  Just  below  the 
critical. 

2.1.2  Equations  of  Motion  for  Turbulent  Plow 

The  equations  that  describe  turbulent  flow  are  derived 
froa  those  of  Section  3*1  through  the  substitution  of  a  tine 
aean  plus  a  fluctuating  component  for  the  fluctuating  quantities. 


The  substitutions  aade  for  incompressible  turbulent 
flow  are  of  the  fora 

ZZ  +U*'  (2-9) 


such  that 


(2-10) 
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Solving  for 

■+£- 

with  Boundary  Conditions 

r  “  protor#  ■  ^rotor  protor 
r  -  pstator'  V  '  0 


There  results 


r  + 


(2-8) 

A  plot  of  V^.  for  the  present  apparatus  at£^  7T rad  lane/ sec 
(30rpa)  is  shown  in  Pig.  2-1  .  Here,  V ^ rotor  "  ft/sec,  and 

7a  •  36,  Just  below  the 
crltloal. 

2.1.2  Equations  of  Motion  for  Turbulent  Plow 

!  The  equations  that  describe  turbulent  flow  are  derived 
from  those  of  Section  3*1  through  the  substitution  of  a  tine 
aean  plus  a  fluctuating  component  for  the  fluctuating  quantities. 


The  substitutions  wade  for  iacoapressible  turbulent 
flow  are  of  the  fora 

U  *  2Z  +u'  (2-9) 


such  that 


(2-10) 


2 
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fro*  which  It  follows  that 


/s  '  =  /  Je^-~  *  O 


(2-11) 


Sinilarly  then  __ 

£  H' 

t’  ?  *?' 


(2-12 ) 


The  result  when  these  are  substituted  Into  the 
Equations  of  Motion  Is  that  the  equations  for  the  turbulent 
case  are  exactly  as  those  for  the  laminar,  except  that: 

1.  V  •  Is  changed  to  and 

to  and 


2. 

which  Is  now 


_  ^  C*> 

Is  changed  to  fy 


Thus 


(2-13) 


(2-13 )a 


The  averaging  procedure  Is  rather  straight  forward  and  Is 
described  by  Goldstein.1  The  tens  act  In  the  equations 

Just  as  the  stresses  due  to  viscous  forces,  and  are  called 


^S.  Goldstein,  Modern  Developments  In  Fluid  Meohanlos  (1st  ed; 
Oxford  at  the  Clarendon  Press,  1938)*  I,  pp.  191*193* 
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Reynolds  stresses.  They  represent  the  transport  of  momentum 
aoross  a  surface  that  Is  due  to  velocity  fluctuations;  for 
example,  the  term^/j^'  represents  the  transport  of  momentum 
across  the  surface  r  «  constant,  by  In  equation 
(2-13  )a. 

For  a  turbulent  compressible  flow  the  result  Is  not 
In  so  neat  a  form. 

Adding  to  equation  (2-12) 

/v3V/' 

results  In  the  same  expressions  as  before  plus  many  new  ones 
resulting  from  double  and  triple  correlations  among  density 
fluctuations  and  turbulent  components  of  the  velocities.  These 
are  listed  by  Hlnze.1  The  resulting  equation  of  mean  motion 
can  no  longer  be  put  Into  the  same  form  as  the  equation  for 
incompressible  turbulent  flow. 


Further,  the  continuity  equation  becomes  with  the 
notation  used 

which,  when  averaged,  for_  steady  flow  becomes 

'© 


The  effect  of  the  density  fluctuations  Is  seen  to 


J.  0.  Hlnze,  Turbulence  (McGraw-Hill  Book  Company,  Inc., 
New  York,  1959),  PP.  2 T>  24. 
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/ 


be  an  apparent  aource  contributing  a  quantity 


To  overcome  these  difficulties.  Professor  K.  0*  Elrod  has 
suggested  In  his  lectures  at  Columbia  the  use  of  the  following: 
Define  a  new  tine  Bean  of  some  flow  property  U*cis 


such  that 


(2-14) 


U=  U 


(2-15) 


Integrating  equation  (2-14)  results  In 

&  *  -f-  J  C£j£2^Eta9a>£ 

A 


777" 


=  Z  +  &r 


(2-16) 


Prom  equations  (2-15)  and  (2-16)  it  Is  seen  that 
U  "  *  U 

_  _  J 

Although  U'sC^  except  for  the  case  of  P  ■  constant 

Applying  the  new  averaging  technique  to  the  Continuity 
Equation,  equation  (2-5) 

m  *  /  4/rri vJ 


^Ting-Yl  Li,  H.  T.  Nagamatsu,  "Effect  of  Density  Fluctuations 
on  the  Turbulent  Skin  Friction  on  a  Flat  Plate  at  High  Supersonic 
Speeds,"  Guggenheim  Aeronautical  laboratory,  California  Institute 
of  Technology,  Pasadena,  California,  Memorandum  No.  11,  Nov.  1, 
1952,  p.  4. 
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when  time  averaged  this  becomes,  for  steady  flow 


A 

/ 

A 


£[^  -  m- 


■///^ 


That  Is 


S/S' 


*  C  onstant  *  0  on  the  stator. 


Therefore  Vr  -  0 


If  this  substitution  Is  made  Into  the  equations  for  momentum , 
the  equations  resulting  are  Identical  In  form  with  the  original 
transient  form  except  that: 

1.  Is  changed  to  V^-  and  to  ,  and 

2 .  is  changed  to  ^  which  Is  now 


$*’  ?;  +/*  Y 

The  terms are  the  virtual  stresses  as 
contrasted  to  S’ ^ for  the  Incompressible  case. 


(2-17) 


In  the  case  of  a  fluid  of  uniform  density  f 

reduces  to 
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2. 1.2.1  The  Sheer  Stress,^. 


(A 


In  aooordanoe  with  the  previous  seotlon  ^  Is 


therefore 


(2-18) 


where  le  developed  below: 

Jr  £/&.%] 

,  -yt® 


The  last  two  terns  represent  Interaction  between 
■olecular  and  turbulent  transport .  They  are  of  higher  order 


so 


than  the  first  tern  and  ean  be  neglected 3 
Equation  (2-18)  becomes 

4^=  +  /iFZ*  <2-X8)a 

If  an  eddy  dlffusivlty  for  momentum  Is  defined  as 


- 


(2-19) 


Prom  a  knowledge  of  the  velocity  distribution  and 


^Ting-Yl  Li ,  and  H.  T.  Nagamatsu,  "Effect  of  Density  Fluctuations 
on  the  Turbulent  Skin  Friction  on  a  Flat  Plate  at  High  Supersonic 
Speeds,"  pg.  3. 

2 

R.  G.  Delssler,  and  A.  L.  Loeffler,  Jr.,  Turbulent  Flow  and 
Heat  Transfer  on  a  Flat  Plate  at  High  Mach. Number  with  Variable 
Fluid  Properties,"  ASME  Paper  Number  55 -A -33#  1955#  p.  20. 

3 

For  example,  if  characteristic  values  are  taken  for  the  velocity, 
length,  density  and  viscosity,  say  and  the  differ¬ 

ential  equations  non-diaenslonallzed  then  in  the  stream. 

On  this  basis  the  molecular  viscous  effects  are  of*  order 
and  interaction  between  fluctuations  and  viscosity  are  of 
order  )  and  ***  negligible. 


ft)  -  y 

the  shear  stress  ^  ^  oan  be  found. 


The  ^  component  of  the  Momentum  Equation,  equation 
(2-1 ),  for  the  present  ease  reduces  to 

-JZZ  &  r  /s» 
or. 


«£> 


constant 


That  Is 

^  ^  firth  '*  (2-21) 

A  knowledge  of  the  shear  stress  at  one  wall  therefore 
yields  the  shear  stress  throughout  the  flow. 


2. 1.2. 2  The  Energy  Equation 


The  differential  equation  for  the  stagnation 
enthalpy  H,  Is  given  as 


#2c 


where  the  stagnation  enthalpy  Is,  neglecting  body  forces. 


and 

Applying  the  definitions  given  In  equations  (2-14)  and  (2-15) 
for  the  time  averages  of  the  quantities  Involved  and  again 
neglecting  all  Interaction  between  moleoular  and  turbulent 
transport,  the  resulting  equation  Is 
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The  term  of  Interest  In  this  equation  Is 


A  •+.  /fa  since  In  cylindrical  coordinates  this  term  Is 

FC  '  —T/  Ta  ~ 

//•  -t-/0*  # 

and  Is  the  radial  heat  flux  from  the  rotor  to  the  stato^j2^^ 


Thus 


(2-24) 


To  formulate  an  expression  for  consider  the 
time  average  ^  iJL 


£ 


-fig 

=*  j&G/r  _  -/t'di' 


Neglecting  Interaction  between  molecular 
and  terms  of  higher  order  magnitude 


and  turbulent  transport 


If  an  eddy  dlffuslvlty  for  heat  Is  defined  as 

^  =  A_  fiEE7' 

u=r?  /■ 


(2-25) 


S3 


(2-26) 


(2-27) 

2.1.3  The  Determination  of 

With  a  knowledge  of  S^jJ**'*  i£^and  Tfr)  one 
can  determine  ^  and 

To  determine  *#/£*,  the  following  equations  are  used: 
Prom  equation  (2-20)1 

From  equation  (2-26) 

Also 

/f&f*  ®"  (2-21) 

and 

f-f-  - grnafr  (2-^) 

Dividing  the  latter  two  equations  results  in 

_  f*r"  if/+  *»/*)  JT 

^Henceforth  the  superscript  " (/  )"  will  be  omitted. 


we  have  for  a"' 

&  “ 

To  obtain  ^//^a  knowledge  of  ^  and  7” is  required. 

The  relation  between  ^  at  K  and  1® 
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OR 


Define  X  as 

r»_i _ L  -  /"%,*,  )dX 

/-1  ?  t  £,/&<*&• 


By  algebra 

&  £,6-*)+' 

Also  g  4 

Since  *  ? 

j,  dT/dXtX  g/L 

7^7^ 


Thus,  a  plot  of 


yields  X. 


(2-28) 


(2-29) 


(2-30) 


(2-31) 


Note,  more  Importantly  that,  since  close  to  the 
wall  fa/r)  -  4  _7/Cr.  Yfi/r  In  that  region  Is  a  straight 

line  then  ^  ■  1,  and  from  equation  (2-30) 

Since  for  air  -  1.4, 

<X  -  1.4 


To  summarize,  what  la  required  for  a  determination 

of  d  Is 
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1.  Proa  equation  (2-20) 


-JT-' 

total  pressure  traverse  in  the  gap  between  the  two  cylinders 


a. 


which  will  be  obtained  froa  a 


b.  which  will  be  obtained  through  a  knowledge 

of  the  shear  stress  at  the  wall  (see  equation  (2-21)).  The 
deteralnatlon  of  %/  Is  described  In  Chapter  VI,  Section  6.4. 


2.  Proa  equation  (2-26) 

a.  or  which  Is  obtained  froa  lnterferoaetrlc 
measurements  of  density  In  the  gap 

b.  fir  which  Is  obtained  froa  fir/  (equation 
(2-27)),  which  In  turn  Is  aeasured. 


Chapter  III 
Experimental  Apparatus 
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3*1  General  Arrangement 

The  apparatus  from  which  the  experimental  data  were 
obtained  consists  of  an  electrically  heated,  motor-driven  rotor 
and  a  water  cooled  stator.  Auxiliary  equipment  was  provided 
for  heater  power  and  cooling  water. 

Provision  was  made  for  the  measurement  of  the  fol¬ 
lowing  : 

1 .  Rotor  temperature 

2 .  Stator  temperature 

3.  Rate  of  flow  of  cooling  water 

4.  Temperature  rise  of  cooling  water 

5.  Heater  power  Input 

6.  Static  pressure  distribution  axially  along 
the  stator  wall 

7.  Radial  total  pressure  distribution  In  the  gap 
between  the  two  cylinders 

3.  Root  mean  square  values  of  the 

radial  and  stream-wise  turbulent  velocity 
fluctuation  components 

9.  Radial  density  distribution  In  the  gap  between 
the  two  cylinders. 

A  schematic  diagram  of  the  apparatus  is  shown  In 

Pig.  3-1. 


The  Instrumentation  of  the  experimental  apparatus 
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Is  discussed  In  Chapter  IV. 

3*2  Inner  Rotsting  Cylinder 

The  Inner  cylinder  Is  shown  schematically  In  Fig. 
3-2.  It  Is  composed  of  a  stainless  steel  shell,  a  heater 
element,  two  end  plates,  a  stainless  steel  shaft,  and  a  slip 
ring  assembly  which  Is  used  for  transmitting  power  to  the 
heater  element,  and  also  for  transmitting  thermocouple  from 
the  thermocouples  embedded  In  the  rotor. 

3.2.1  Rotor  Outer  Shell 

Because  of  the  temperature  requirements  (up  to 
800  F)  and  strength  requirements  (up  to  10,000  rpm)  the  outer 
shell  Is  made  of  austenitic  stainless  steel,  available  as 
commercial  plpe.^  The  shell  was  machined  to  a  length  of  4 
Inches,  an  outside  diameter  of  3*320  Inches,  ±0.0005  Inches, 
and  an  Inside  diameter  of  2.938  Inches  ±.001  inches.  The 
surface  finish  of  the  outside  surface  Is  5  microinches  rms, 
as  measured  with  a  Prof  Hornet  er.^ 


Four  thermocouples  are  located  at  ninety  degree 
Intervals  in  the  wall  of  the  outer  shell,  two  at  an  axial 
distance  of  2-1/4  inches ( */Ls  '/z)  and  two  at  an  axial  distance 
of  one  inch  /X/l  =  >  ®8  measured  from  the  end  of  the 


Austenitic  stainless  steel  is  typified  by  high  strength  at 
elevated  temperatures  and  has  excellent  scale  resistance  where 
conditions  of  high  temperatures  and  oorroslon  are  met. 

2 Manufactured  by  Physicists  Research  Corp.,  Ann  Axtoor,  Mlohlgan. 
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rotor  nearest  the  slip  rings. 

The  holes  in  which  the  thermocouples  sre  placed  were 
drilled  axially  to  minimize  errors  In  measurement  which  would 
be  csused  by  temperature  gradients  and  insufficient  length 
of  thermocouple  wire  In  the  parent  metal.  The  installation 
was  made  according  to  procedures  outlined  by  Baker,  Ryder  and 
Baker.1 


A  further  discussion  of  the  thermocouple  Installation 
and  the  slip-ring  arrangement  Is  found  In  Chapter  IV. 

3.2.2  Rotor  Heating  Element 

The  rotor  heating  element  Is  composed  of  a  stainless 
steel  tube  4-1/4  Inches  long,  outside  diameter  of  2  Inches  and 
a  3/32  inch  wall  thickness.  This  tube  Is  covered  with  two 
layers  of  mica  sheet  approximately  0.004  Inches  thick  and 
wound  with  heater  wire.  The  wire  used  was  "Alloy  K",  pur¬ 
chased  from  C.  0.  Jellif  Mfg.  Corp.,  Southport,  Connecticut, 
and  had  a  resistance  of  0.48  ohms/ft. 

Winding  was  done  on  an  engine  lathe,  holding  the 
tube  in  a  three -Jaw  chuck  and  the  wire  In  the  tool  post  where 
It  was  constrained  so  as  to  be  under  tension.  The  carriage 


XH.  D.  Baker,  E.  A.  Ryder,  and  N.  H.  Baker,  Temperature 
Measurement  In  Engineering  (New  York:  John  Wiley  and  Sons, 


1953),  Vol.  I 
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was  then  given  such  a  motion  aa  to  produoe  six  wraps  per  inch, 
corresponding  to  a  heater  resistance  of  13  ohms. 

Control  of  heater  power  during  a  run  was  by  means 
of  a  Powerstat1 .  Provision  was  made  for  measurement  of  heater 
power  and 

3.2.3  Rotor  Shaft 

The  rotor  shaft  is  shown  in  Pig.  3-3.  The  shaft 
material  is  stainless  steel. 

Axial  grooves  placed  about  the  circumference  of  the 
shaft  receive  thermocouple  leads  from  the  outer  shell  of  the 
rotor  and  thus  enable  the  shaft  to  pass  through  the  stator 
end  plates  (see  Section  3*4)  with  a  clearance  of  0.005  inches. 

The  shaft  was  ground  to  an  accuracy  of  0.0001  Inches 
as  required  by  the  ball  bearings  selected  for  the  shaft. 

The  ball  bearings  were  "D-B  Mounted",  a  type  of 
mounting  whereby  the  backs  of  two  bearings  are  preground  and 
then  placed  together.  This  mounting  eliminates  the  eccentricity 
in  the  bearings  and  is  generally  recommended  by  bearing  manu¬ 
facturers  for  high  speed,  high  accuracy  grinding  spindles. 

The  light  bearing  loads  (about  three  pounds  per 
bearing)  and  the  method  of  mounting  assures  high  rotational 
speeds  with  a  minimum  of  run-out  (  0.0001  inches  on  the  shaft). 

^Manufactured  by  the  Superior  Electric  Company  (Bristol, 
Connecticut). 
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Figure  3- 3b 

stainless  steel  Rotor  shaft 
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3.2.4  Rotor  End  Plates 

The  end  plates  provided  to  hold  together  the  rotor 
assembly  are  shown  In  Fig.  3-4. 

The  requirements  of  a  suitable  end  plate  material 
are  1.  That  It  have  sufficient  strength  at  high  temperatures 
to  withstand  the  centrifugal  stresses  at  high  speeds. 

2.  That  It  have  low  thermal  conductivity  to  Inhibit 
heat  losses  from  the  outer  rotor  shell  and  heater  to  the 
shaft. 

t 

3.  That  It  be  easily  machined,  since  holes  and  grooves 
were  required  for  the  passage  of  thermocouple  and  heater  leads 
from  the  rotor  outer  shell  and  the  heater  to  the  slip  rings. 

A  first  attempt  was  made  to  use  a  product  of  the 
American  Lava  Corporation  of  Chatanooga,  Tennessee,  called 
Lava  Grade  A,  an  aluminum  silicate.  This  product  was  easily 
machined  before  firing.  By  trial  and  error  the  decrease  In 
size  because  of  firing  was  determined  and  It  was  possible  to 
have  the  end  plate  In  final  form  to  size  within  0.001  Inches. 
This  material  was  rejected,  however,  because  of  cracks  which 
developed  during  assembly  because  of  loads  Imposed  by  fastening 
screws . 


A  second  material,  Marlnlte,  manufactured  by  the 
Johns -Manvllle  Corporation  of  New  York,  New  York  was  found  to 
be  satisfactory.  The  containing  ring  A  and  the  spider  like 
member  B  In  Fig.  3-4  were  provided  so  that  the  Marlnlte  would 
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not  be  a  load  bearing  member.  The  two  end  plates  and  their 
containing  rings  are  each  1/4  Inch  wide. 

Further,  only  one  end  plate  was  fastened  to  the 
shaft  by  means  of  the  Teflon  collar  C  .  The  second  end 
plate  made  no  connection  with  the  shaft  to  allow  for  axial 
thermal  expansion  of  the  rotor  assembly  to  be  taken  up  In  this 
direction  (much  as  the  floating  tube-sheets  In  commercial  heat 
exchanger  design). 

3.2.5  Slip-Ring  Assembly  for  Rotor  Temperature  Measurement 

The  slip-ring  assembly  for  rotor  temperature  measure¬ 
ment  Is  shown  in  Fig.  3-5.  The  brushes,  brush  holders  and  coin 
silver  slip  rings  were  purchased  from  the  Oraphlte  Metallizing 
Corporation  of  Yonkers,  New  York.* 

The  slip  rings  were  pressed  over  a  bakelite  bushing, 
as  shown  In  Fig.  3-5,  with  bakelite  spacers  between  the  rings. 
The  bushing  was  grooved  to  receive  copper  lead  wires  which  at 
one  end  were  soldered  to  the  coin  silver  slip  rings  and  at  the 
other  end  to  the  Iron  and  constantan  thermocouple  wires  at 
points  in  the  ambient  air  external  to  the  rotor  and  rotating 
with  It.  The  copper  lead  wires  from  the  silver  GRAPHALLOY 
brushes  were  returned  to  fora  compensating  Junctions  with  Iron 


*From  Graphite  Metallizing  Corporation,  catalog  No.  3»  General 
Brushes:  Silver  GRAPHALLOY,  5/32  Inches  square. 

Brush  Holders:  No.  11509 

Coin  silver  slip  rings:  No.  1165-27,  1-11/16  lnohes  O.D., 
machined  to  5/32  Inches  wide. 


A  39m by  Of&w/nj 
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and  with  constantan  at  points  close  to  the  path  of  rotation 
of  the  previously  mentioned  copper-and>lron  and  copper-and- 
constantan  Junctions. 

The  brushes  are  spring  loaded  In  such  a  way  that  so 
that  contact  with  the  slip  rings  was  only  momentary  during 
measurement.  In  this  manner  spurious  thermoelectric  effects 
are  minimized.  A  more  detailed  discussion  of  the  accuracy 
of  emf  measurement  with  the  slip  ring  assembly  Is  given  In  the 
next  chapter. 

3.3  Outer  Stationary  C Under  and  Cooling  Jacket 

The  outer  stationary  cylinder  and  cooling  Jacket 
Is  shown  In  Pig.  3 -6a  and  b. 

3«3»1  The  Outer  Cylinder 

The  outer  cylinder  is  made  of  austenitic  stainless 
steel  available  as  commercial  pipe  and  machined  to  an  Inside 
diameter  of  4.074  t.001  Inches.  It  was  ground  to  a  surface 
finish  of  7  microinches  rms. 

Seven  holes,  0.020  Inches  In  diameter  were  drilled 
through  the  stator  walls  as  shown  In  Pig.  3-6  to  allow  the 
Insertion  of  0.020  Inches  outside  diameter,  0.010  Inches  Inside 
diameter  stainless  steel  hypodermic  tubing'1,  to  serve  as  axial 
static  pressure  taps. 

^No.  25  Stainless  steel  hypodermic  tubing,  purchased  from 
Becton  and  Dickenson  and  Co.,  Inc.,  Rutherford,  New  Jersey. 
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Outer  stationary  Cylinder 


and 

Cooling  Jacket 


Figure  3-6  b 
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The  hypodermic  tubing  was  soldered  Into  4-40 
stainless  steel  screws. 

The  stator  holes  were  counter  bored  and  tapped  to 
receive  the  tubing  and  Its  fitting.  Non-hardening  Permatex 
(No.  2)  was  used  at  the  Joint  to  prevent  cooling  water  leakage. 
The  stainless  steel  tubing  was  counter  sunk  at  the  wall  of 
the  stator  0.010  Inches  to  eliminate  the  effect  of  disturbances 
to  the  static  pressure  at  the  stator  wall  caused  by  the  presence 
of  the  tap.  This  effect  Is  discussed  In  Chapter  IV. 

Two  holes  0.020  Inches  In  diameter  were  drilled 
through  the  stator  wall  at  a  radial  position,  22-1/2  degrees 
from  the  static  pressure  taps;  one  at  an  axial  distance  of  1-1/2 
Inches,  the  other  at  an  axial  position  of  2-1/4  inches.  The 
total  pressure  probes  were  later  inserted  Into  these  holes. 
Suitable  packing  glands  were  silver  soldered  to  the  outer 
surface  of  the  stator,  concentric  with  these  holes  and  are 
shown  in  Pig.  3 -6c.  The  purpose  of  these  glands  is  to  prevent 
cooling  water  leakage  into  the  gap  between  the  concentric 
cylinders . 

Two  chromel-alumel  thermocouples  were  embedded  Into 
the  stator  at  approximately  the  same  radial  position  but  at 
different  axial  positions;  one  at  an  axial  distance  of  2-1/4 
Inches,  the  other  at  an  axial  distance  of  1-1/2  inches.  The 
details  of  the  thermocouple  installation  are  given  In  the  next 
Chapter. 


/-/'Jcfre  3-6  <= 
'fetal  flkssunz /joke 
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3*3*2  Cooling  Jacket 

The  cooling  Jacket  was  made  from  6  Inch  Schedule  80 
brass  pipe.  Water  Inlet  and  outlet  are  through  3/8  Inch  copper 
tubing  to  pipe -thread  brass  elbows.  In  each  of  which  Is  arranged 
a  thermocouple  well  and  thermocouple  as  shown  In  Pig,  3-7. 

These  thermocouples  are  for  the  measurement  of  the  temperature 
difference,  entering  and  leaving,  of  the  cooling  water. 

As  a  precaution  against  heat  losses  from  the  Jacket, 
the  Jacket  Itself  was  oovered  with  Johns -Manville  85$  carbonate 
of  magnesia,  15$  asbestos  Insulation.  Losses  were  further 
minimized  by  having  the  average  ooollng  water  temperature 
approximately  equal  to  room  temperature. 

3*3*3  End  Plates 

Two  atator  end  plates  were  provided  to  enclose  the 
rotor-stator  assembly.  These  are  shown  In  Fig.  3-8.  These 
were  made  of  Afartnite . 

Fused  quarts  windows  are  seated  In  the  reotangular 
openings  shown  to  provide  a  light  path  for  Interferometric  study 
of  the  density  distribution  In  the  gap  between  the  rotor  and 
stator. 

3.4  Bearing  Housings  andSese  Plate 

Pig*  3-9  Is  a  drawing  of  the  bearing  housings* 

The  base  plate  has  levelling  sorews  on  whloh  the 
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Fig  3-3 
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stator  assembly  Is  seated.  The  levelling  screws  enable  ad¬ 
justment  of  the  stator  position  so  that  It  can  be  made  con¬ 
centric  to  the  rotor. 

3.5  Motor 

The  rotating  cylinder  Is  driven  by  a  D.  C.  motor. 
Alternating  current  Is  led  through  a  constant  voltage  trans¬ 
former  and  a  rectifier  to  the  motor.  Speed  adjustment  Is  made 
by  means  of  slide  wire  rheostates  In  the  armature  and  field 
circuits  of  the  motor. 

A  Strobotach1  was  used  to  determine  variations  from 
the  desired  speed.  The  Strobotach  was  not  used  as  a  speed 
measuring  device  because  of  its  drift.  A  chronometrlc  tachome¬ 
ter  placed  against  the  rotating  shaft  was  used  to  set  the 
Strobotach  to  the  desired  glimpse  frequency.  The  setting  of 
the  Strobotach  was  checked  through  speed  measurements  with  the 
Chronometrlc  tachometer  every  fifteen  minutes. 

3.6  Cooling  Water  Arrangement 

The  cooling  water  arrangement  is  shown  in  Pig.  3-10. 

Hot  and  cold  water  were  mixed  in  a  chamber  made  from  2  inch, 

Schedule-40  steel  pipe.  A  pressure  regulating  valve  was 

supplied  for  both  the  hot  water  and  cold  water.  After  mixing, 

the  water  was  led  through  another  pressure  regulating  valve 

to  a  rotameter.  The  rotameter  was  used  to  determine  variations 

from  the  desired  flow  settings.  Variations  from  the  settings 
were  infrequent,  generally  occurring  before  noon  and  at  4:30  P.M. 

■'■Manufactured  by  the  General  Radio  Company  (Cambridge,  Mass.). 
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when  excessive  demands  were  placed  on  the  Pupln  Building  water 
supply.  At  these  times,  adjustment  was  required  In  control 
valve  settings  to  maintain  the  desired  temperature  and  flow 
rate.  The  adjustments  were  minor  In  nature.  Prom  the  rota¬ 
meter  the  water  was  led  through  Tygon  tubing  to  the  test- 
section  and  returned  from  the  test-section,  again  In  Tygon 
tubing,  to  a  pall  placed  on  a  lever  scale.1  Here  the  flow 
rate  was  measured  by  weighing  the  water  flow  over  a  time 
Interval.  Thus,  accurate  calibration  of  the  rotameter  was 
not  required. 


^he  water  outlet  from  the  rotameter  is  through  an  expanded 
tube  fitting  from  which  only  several  inches  of  copper  tubing 
extends.  A  tight  Joint  was  obtained  between  the  Tygon  tubing 
and  the  piece  of  copper  tubing  by  selecting  the  Tygon  size 
about  1/8  inch  smaller  in  inside  diameter  than  the  outside 
diameter  of  the  copper  tubing.  The  Tygon  was  then  immersed 
in  Acetone  for  about  one -half  hour  causing  it 3  diameter  to 
increase.  It  was  then  placed  over  the  copper  tubing.  When 
the  Acetone  evaporated  from  the  Tygon  a  "shrink"  fit  resulted 
between  the  Tygon  and  copper. 
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CHAPTER  IV 
Instrumentation 

4.1  The  experimental  apparatus  described  in  Chapter  III  was 
instrumented  so  that  ultimately  the  following  were  known: 

1.  Overall  rate  of  heat  transfer  from  the  rotor  to  the 
stator  under  varying  conditions  of  rotor  speed  and  temperature. 
The  required  measurements  were  rotor  and  stator  temperature , 
rate  of  flow  of  cooling  water,  and  temperature  rise  of  cooling 
water. 

2.  Velocity  distribution  in  the  gap  between  the  rotor 
and  stator.  The  required  measurements  were  total  pressure 
distribution,  static  pressure  at  the  stator  surface,  and 
density  distribution. 

3.  Temperature  distribution  in  the  gap.  The  measurements 
taken  were  density  distribution  in  the  gap  as  determined  with 

a  Mach-Zehnder  interferometer. 

4.  Whether  Taylor  vortices  persisted  in  the  flow  at 
rotational  speeds  of  the  experiment.  The  required  measures 
were  made  with  a  hot-wire  anemometer. 

The  measurements  taken  can  then  be  grouped  as  temperature 
measurements,  pressure  measurement,  water  flow  rates,  fluctu¬ 
ating  velocities  and  density  measurements . 
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4 ,2  Temperature  Measurements 

Wherever  possible ,  thermocouples  were  used  to 
measure  temperature.  The  emf  outputs  of  the  thermocouples 
were  measured  with  a  null-balance  potentiometer.  Thermocouple 
Installation  and  error  analysis,  cold  junctions  and  Indication 
circuitry  were  generally  as  outlined  by  H.  0.  Baker  et.al.1 

4.2.1  Stator  Temperature 

The  stator  temperature  was  measured  with  number  30 
B  and  S  gage,  glass  insulated  chrome 1 -a lumel  thermocouple  wire 
obtained  from  the  Leeds  and  Northrup  Company  of  Philadelphia, 
Pennsylvania . 

Two  thermocouples  were  placed  In  the  stator  wall 
as  shown  In  Pig.  4-1,  one  at  an  axial  location Vz  the 
other  at  an  axial  location//^.  ^  .  They  were  displaced 
from  each  other  In  the  ■£  -direction  by  22-1/2  degrees. 

The  Installation  was  made  exactly  as  outlined  by 

o 

Baker  except  that  the  drill  size  used  was  0.052  Inches  rather 
than  0.046  Inches.  The  junction  on  these  thermocouples  as 


■^H.  D.  Baker,  E.  A.  Ryder,  and  N.  H.  Baker,  Temperature 
Measurement  In  Engineering  (New  Yorks  John  Wiley  and  Sons, 
1953),  Vol.  I. 

2H.  D.  Baker,  et.al..  Temperature  Measurement  In  Engineering, 
pp.  143-145. 
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on  all  those  used  was  made  by  the  arc  method.* 

4.2.2  Rotor  Thermocouples 

The  rotor  temperature  was  measured  with  number  30 
B  and  S  gage  glass  insulated,  duplex  lron-constantan  thermo¬ 
couple  wire  obtained  from  the  Leeds  and  Northrup  Company  of 
Philadelphia,  with  error  of  1%  of  indication.  The  installation 
was  made  as  for  the  stator  thermocouples. 

4. 2. 3.1  Slip  Ring  Assembly 

The  measurement  of  the  emf  of  the  rotor  thermo¬ 
couples  was  made  through  a  brush  and  ring  type  slip  ring. 

Other  arrangements,  such  as  a  mercury  slip  ring,  or  magnetic 
coupling  between  a  rotating  and  stationary  coll  were  rejected 
because  of  the  bulk  of  the  installation  for  the  former  and 
the  additional  investment  in  electronic  equipment  required 
for  the  latter.  The  brush  and  ring  type  had  the  advantage 
that  the  component  parts  were  available  as  stock  items  from 
the  supplier,  that  the  installation  was  mechanically  simple 
and  that  several  successful  installations  were  described  in 
detail  in  the  literature. 

One  such  installation  is  that  discussed  by 


D.  Baker,  et.al..  Temperature  Measurement  in  Engineering, 


p.  50 
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R.  Chaplin1  and  It  was  this  design  that  was  paralleled  in  the 
present  installation.  Chaplin  reports  that  his  design  made 
possible  temperature  measurements  with  errors  of  less  than 
1°C,  with  a  chromel-alumel  thermocouple. 

4. 2. 3. 2  Design  of  Brush  and  Ring  Type  Slip  Ring 

The  most  severe  problem  to  be  overcome  in  the 
design  of  a  brush  and  ring  type  slip  ring  assembly  is  that 
of  the  electrical  noise  caused  by  the  sliding  contact  between 
the  brush  and  the  ring.  This  "noise"  is  of  two  types;  one 
caused  by  the  variation  of  contact  resistance,  the  other 
caused  mainly  by  spurious  thermal  emfs  from  hot  spots  caused 
by  frictional  heating.  Since  a  null  balance  method  of  emf 
measurement  was  used  in  the  present  experiment  the  variation 
of  contact  resistance  is  of  negligible  consequence  so  long 
as  it  remains  small.  The  second  type  of  noise  was  studied 
in  detail  by  Horton^  who  measured  the  generated  noise  voltages 
of  different  brush  materials  such  as  natural  graphite,  silver- 
graphite,  with  both  high  and  low  silver  content,  and  rhodium- 
graphite  acting  on  rings  of  silver,  rhodium,  gold,  copper, 

1r.  Chaplin,  "Multichannel  Slip  Rings  for  Stress  and  Temperature 
Measurement"  National  Gas  Turbine  Establishment,  Pyestock,  Hants, 
Great  Britain,  Memorandum  Number  M.289,  (December,  1956),  p.  7« 

p 

B.  M.  Horton,  "Sliding  Contacts  to  Transmit  Signals,"  Rev.  Scl. 
Insts. ,  Vol.  20,  Number  12,  (December  1949),  pp.  930-932. 
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and  electrographltle  carbon. 

Horton  found  that  1)  "surfaoea  covered  with  corrosion 
films  give  much  higher  values  of  generated  noise  voltage 
than  do  clean  surfaces",  2)  "All  of  the  contact  confeinatlons. . • 
[had].. .a  generated  noise  voltage  that  Increased  with  speed", 
and  he  concludes  by  stating  "metal -graphite  brushes  operating 
with  low  frletlon  on  clean  metal  surfaces  can  yield  very  low 
values  of  generated  noise  voltage  In  the  frequency  range 
considered  [0.5  to  200  o.p.s,]". 

Por  a  clean  silver  ring  acted  on  by  2  brushes  of 
fine  texture,  high  oarbon  content  silver  graphite  brushes 
with  a  normal  force  of  50  grams,  and  a  sliding  speed  of 
35  cm/sec.,  (about  350  RPM  on  the  present  apparatus),  Horton 
reports  a  generated  noise  of  0.3  microvolts.  This  value  was 
the  lowest  reported  In  his  Table  1.  Silver  rings  with  silver- 
graphite  brushes  of  low  silver  content  were  therefore  chosen 
for  the  present  assembly. 

At  first,  the  Iron  and  constantan  of  the  rotor 
thermocouples  were  soldered  directly  to  the  silver  slip  rings. 
Large  errors  were  observed  In  the  output  emf  when  the  brushes 
and  rings  were  In  continuous  contact.  Allowing  compressed  air 
to  flow  over  the  rings  reduced  this  error  to  more  reasonable 
proportions.  When  the  brushes  were  allowed  to  contact  the 
rings  for  a  period  of  1  or  2  seconds  as  necessary  to  make 
measurements,  errors  In  emf  varying  with  the  time  of  contact 
were  still  observed. 

Consequently,  the  Iron  and  the  constantan  leads  from 
the  rotor  thermocouples  were  joined  to  copper  leads  at  points 
In  the  ambient  air  external  to  the  rotor  and  rotating  with  It. 
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These  form  two  new  thermocouples,  copper -constantan  and 
copper- Iron ,and  require  compensation.  This  compensation  Is 
achieved  by  having  the  copper  leads  from  the  brushes  returned 
to  form  the  necessary  Junctions  with  Iron  and  constantan  at 
points  close  to  the  path  of  rotation  of  the  copper-iron  and 
copper-constantan  Junctions.  This  arrangement  Is  shown 
schematically  In  Fig.  4 SL  and  comprises  the  final  arrangement. 

Copper  was  chosen  as  the  Intermediate  material 
between  the  rings  and  the  rotor  thermocouples  because  of  the 
low  thermal  emf  generated  between  copper  and  sliver  as  compared 
with  iron  or  constantan  with  silver.1  Wlth  coPPer  as  the 

^he  data  for  the  following  was  obtained  from  the  International 
Critical  Tables  (first  ed.,  McGraw-Hill  Book  Co.,  Inc.,  New 
York,  1929)  Vol.  VI,  pp.  214-215. 

For  example,  if  theeemf  of  a  thermocouple  composed 

of  metals  A  and  B,  with  one  Junction  at  0  deg.  C,  the  other 

at  t  deg  C  is  denoted  by  E  then  for  the  metals  mentioned 

a  6 

;  *  /3. 7s E  -  3.&J5 £\/o~2m  +2.8SC /.  « /oS 

*2S5i £ 

Thus  at  38  deg.  C  (about  100  deg.  F)  the  "thermo¬ 
electric  power"  El  for  each  of  these  is 

°/Au^  -  0-0"  -xvr/oeyC-  * 

“  0-0*1  w/<A*C  *  0021  -Hur/AgP 
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Intermediate  metal  between  the  rotor  thermocouples  and  the 
slip  rings  the  error  In  Indication  was  reduced  to  the  order 
of  three  degrees  F.  even  with  the  brushes  In  continuous  contact 
at  4000  to  6000  rpra. 

Details  of  tests  run  to  determine  the  stated  error 
are  reported  In  Appendix  2ZT . 

4.2.3  Reference  Junction 

The  reference  Junction  used  for  the  thermocouples 
was  a  quart  thermos  bottle  filled  with  ice  made  from  distilled 
water  with  enough  water  added  so  as  to  Just  fill  the  interstices 
between  the  pieces  of  ice.  The  two  thermocouple  Junctions 
one  lron-constantan  theother  chrorael-alumel,  were  each  placed 
in  a  6  mm.  glass  tube,  sealed  at  one  end  and  filled  with 
moisture-free  kerosene.  These  were  immersed  into  the  thermos 
flask  to  a  depth  of  about  6  inches  with  the  sealed  end  of 
the  tube  about  3  inches  from  the  bottom  of  the  flask. 

The  thermos  in  turn  was  surrounded  by  slag  wool  for 
insulation  purposes.  The  ice  was  replenished  before  every 
run. 

4.2.4  Temperature  Indication 

A  null-balance  Leeds  and  Northrup  No.  8662  Portable 
Precision  Potentiometer  was  used  as  the  indicating  Instrument 
for  the  rotor,  stator  and  cooling  water  thermocouples.  The 
Instrument  has  two  scales,  one  0  to  16  millivolts,  theoobher 
0  to  80  millivolts,  with  provision  for  manual  cold  Junction 
compensation.  This  latter  was  not  used  since  a  separate 
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reference  Junction  was  maintained  during  the  experiments. 

The  experiments  reported  here  required  the  use  of  only  the 
lower  scale.  The  smallest  division  on  that  scale  Is  0.01 
millivolts  and  can  be  read  with  confidence  to  1/2  division, 
or  0.005  millivolts,  or  about  1/4  deg.  P  for  a  chrome 1 -a lumel 
thermocouple . 

4.2.5  Thermocouple  Corrections 

Besides  the  error  caused  by  the  slip-ring  Installation, 
consideration  must  be  given  to  errors  caused  by  1)  conduction 
along  the  leads  and  2)  location  of  the  thermocouples  relative 
to  the  rotor  and  stator  surfaces. 

4. 2. 5.1  Conduction  Along  the  Leads 

This  error  Is  a  result  of  heat  transferred  from 
the  point  of  Interest  to  the  ambient  using  the  thermocouple 
wires  as  the  heat  path.  Since  the  transfer  of  heat  Implies 
a  temperature  difference  the  thermocouple  junction  may  not 
be  at  the  same  temperature  as  the  surrounding  metal  of  the 
rotor  or  stator. 

Baker1  gives  an  expression  for  the  heat  transferred 
from  the  leads  to  the  ambient  as 


D.  Baker, 


et.al..  Temperature  Measurement  In  Engineering, 


Vol.  I,  p.  71 
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£'7Z7C4A  j#  (4-D 

where  fa  Is  the  heat  transferred,  Btu/hr, 

C  Is  the  circumference  of  the  leads,  ft 
£  Is  the  mean  value  of  the  thermal 
conductivity  of  the  thermocouple  lead  materials, 
Btu/hr-ft2 -degP/f  t 

A  Is  the  cross-sectional  area  of  the  leads, 

ft2, 

U  Is  the  effective  value  of  the  surface 

2 

boundary  conductance  referred  to  C,  Btu/hr-ft  - 
deg  F  and 

Is  the  surface  temperature  of  the  body 
minus  the  ambient  temperature. 

The  value  for  Is  given  with  sufficient  accuracy 
by  averaging  ^A  for  each  of  the  wires. 


The  value  of  U  Includes  the  thermal  resistance 
of  the  glass  sheathing  about  the  wires. 


The  quantity  must  be  supplied  by  the  parent 
metal.  One  then  computes  the  temperature  difference  between 
the  point  of  interest  and  the  temperature  of  the  sensitive 
Junction.  This  difference,  the  error,  is  given  by 


Error 


each  m (,  +  coth  vlJL 


( 4-2) 
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where 

The  primed  quantities  refer  to  the  parent 

metal. 

Is  the  depth  of  Immersion,  ft,  and 
Is  the  thermal  gradient,  deg. /ft,  along 
the  leads  In  the  parent  metal. 

Looking  at  these  equations  qualitatively,  one 
sees  that  the  error  can  be  reduced  by 

1.  Inhibiting  heat  transferred  to  the  ambient,^ 
by  suitable  choice  of  small  wires,  low  thermal  conductivity 
metals,  and.  If  necessary,  additional  Insulation  over  the  leads. 

2.  Encouraging  that  heat  transfer,^  which  must 
occur  Inside  the  parent  metal,  by  assuring  good  thermal  con¬ 
tact  between  the  parent  metal  and  the  thermocouple,  and  by 
providing  a  large  depth  of  immersion,^  .  These  allow 

to  be  transferred  with  small  temperature  difference  and  hence 
small  error. 

4. 2. 5.1  Estimation  of  the  Conduction  Error 

4. 2. 5. 1.1  Stator  Thermocouple 

A  sketch  of  the  thermocouple  Installation  used 
given  in  Pig.  4-/. 

The  hourly  heat  leaving  the  installation  is  given 
by  eq.  (4-1). 
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Where  90  -  (Tstator  -  Tcooling  water)  “  20  &•***«*»  maximum, 

C  »  TT  (O.Ol4l)/2  »  0.0037  ft.  based  on  the  equivalent 
diameter  of  the  thermocouple  metals,  ^  A  -  0.4  x  10"4 
Btu/hr  -  degP/ft,  and  V  -  8  BtVhr-ft2-degF.  This  latter 
value  Is  obtained  by  using  a  value  of  &  ,  the  surface  con¬ 
ductance  between  the  leads  and  the  water  of  8.8,  referred  to 
the  diameter  of  the  wire,  and  considering  the  thermal  resistance 
of  the  glyptal-coated  glass  braid  Insulation  ( 0.2)  over 
the  leads.  The  thermal  resistances  of  the  latter  can  be 
calculated  by  the  conventional  equations. 

Using  these  values,  ^  -  20  8(0.0037) (0.4xl0"4) 

■  0.021  Btu/hr 


To  find  7r i  of  eq. 


(4-2), 


.  ok'/iW 


is  is  the  overall  heat  transfer  conductance  between 
the  thermocouple  leads  and  portions  of  the  body  sufficiently 
distant  so  that  the  Isotherms  there  are  not  disturbed  by  the 
Installation.  Consequently,  It  Is  composed  of  two  parts;  one, 
the  conductance  between  the  thermocouple  leads  and  the  wall 
of  the  hole,  this  volume  being  filled  with  copper  cement;  the 
other  part  between  the  wall  of  the  hole  and  those  regions 
mentioned. 

t 

V  Is  then  given  as 
For  a  value  of  Jce9Bnt  -  2, 


r(JL  4.  vfc 


''^'cement  ‘  body 
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//  / 

Prom  Baker's  Chart  A,  Pig.  7-6,  /0L  .  »  /» _ 

'  body  '  owo 

since  the  depth  of  Immersion  Is  1  Inch  for  the  stator^  Z7* s  /eoo 


am* 


~  -***•/ 


25.4. 


The  value  of  csch  m &  In  eq.  (4-2)  Is  negligible, 
and  the  error  In  the  absence  of  gradient ^^~is  therefore 
negligible. 


The  temperature  gradient  In  the  stator  from  the 
inside  wall  to  the  outside  wall  Is  given  by 

/Q  . 

The  maximum  /k  to  be  passed  Is  of  the  order 
2000  Btu/hr.  The  thermal  conductivity,  7$  for  stainless 
is  15 .  Thus , 

2000  -  133*3  degrees/ft. 

1 T 


and  the  temperature  difference  between  the  two  walls  Is, 
since  the  stator  cylinder  is  1/4  inch  thick, 

JZ  -  133  x  0.25  -  2.77  degrees 

The  gradient^,  in  the  direction  of  the  thermocouple 
is  therefore  2.77  degrees/inch  -  33.24  degrees/ft. 

With  this  value  of  g  ,  length  of  Immersion^ 

=  1  inch  -  0.0833  ft,  and>7 "  305, 


( 
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the  conduction  error  due  to  Is  0.13  P>  And  the  total 
conduction  error  is  0  +  0.13  •  0.13  F* 

4. 2. 5. 1.2  Rotor  Thermocouples1 

Por  the  thermocouples  embedded  in  the  rotor  the 
conditions  inside  the  parent  affecting  the  aceuracy  of  the 
thermocuples  are  the  same  except  that  f  ■  0  since  the 
thermocouples  were  placed  axially,  and  is  higher  because 
of  higher  1^  ,  and  V  .  However,  the  value  of  is  unchanged, 

and  is  1  inch  for  two  of  the  thermocouples  and  2-1/4  inches 
for  the  other  two.  Thus,  Trr>&  is  of  the  same  order  of  magnitude 
as  for  the  stator  thermocouple  and  the  conduction  error  is 
negligible. 

4. 2. 5. 1.3  Cooling  Water  Thermocouples 

Similar  considerations  show  this  error  to  be 
negligible . 

4. 2. 5. 2  The  Error  due  to  Location  in  the  Stator  or  Rotor  Wall 
For  the  stator  It  was  found  thataT  across  the  wall 
was  2.77  degrees. 

Measurements  show  that  the  sensitive  point  of  the 
thermocouple  was  approximately  0.040  Inches  from  the  wall  of 
interest.  The  correction  to  be  made  is 

2_j2£  x  0.040  m  o.4  degrees  at  the  most 

0.25 

extreme  conditions  encountered  during  the  experiment. 

This  correction  was  applied  throughout  the 

1In  the  computation  of  the  conduction  errors,  perfect  oontaot  was 
assumed  between  the  thermocouple,  copper  cement,  and  the  solid 
body.  This  assumption  may  be  incorrect  for  one  of  the  rotor 
thermocouples.  See  Section  6.1.3* 
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computations  for  overall  heat  transfer. 

Por  the  rotor  thermocouples  the  heat  transferred 
per  square  foot  of  surface  Is  of  the  same  order  of  magnitude. 
Since  their  radial  location  Is  In  the  middle  of  the  wall 
slightly  higher  errors  result  and  are  of  the  order  1.5  F  at 
rotor  temperatures  of  400  P.  A  correction  of  this  error 
was  made  In  the  computation  of  overall  heat -transfer. 

4.3  Pressure  Measurements 

4.3.1  Total  Pressure  Measurement 

Total  pressure  probes  were  used  to  determine  the 
streamwlse  velocity  In  the  gap  between  the  rotating  and 
stationary  cylinder. 

4.3. 1.1  Manufacture  of  the  Total  Pressure  Probe 

The  total  pressure  probes  were  made  as  follows: 

A  random  length  of  stainless  steel  hypodermic  tubing 
with  dimensions  0.020  inches  outside  diameter  and  0.010  Inches 
Inside  diameter  was  deburred  on  Its  end  with  a  fine  honing 
stone . 


The  tube  was  washed  In  acetone  and  blown  clear  with 
compressed  air.  One  end  of  the  tubing  was  then  heated  to  a 
dull  red  to  prevent  splitting  of  the  tubing  during  subsequent 
bending  operations.  A  sliver  of  0.002  Inches  shim  stock 
was  Inserted  to  a  depth  of  approximately  1/8  Inches  Into  the 
hole. 


77 


A  snail  orescent  wrench  was  clamped  shut  over  the 
end  In  which  the  shim  stock  was  Inserted ,  flattening  the  tube 
about  the  shim  stock  for  a  length  of  approximately  1/8  Inches. 
This  flattened  portion  was  then  bent  90  degrees  to  the  main 
stem,  and  the  wrench  was  removed. 

The  piece  of  shim  stock  was  removed,  and  the  probe 
tip  honed  to  the  desired  dimensions.  Measurements  were  made 
of  the  probe  tip  dimensions  with  a  toolmaker's  microscope, 
to  an  accuracy  of  0.0001  Inches.  Pig.  4-4  shows  a  probe  used 
for  the  present  experiment. 

The  pressure  probes  were  mounted  In  the  apparatus 
shown  In  Pig.  4-5.  Pressure  In  the  probe  was  transmitted 
through  the  manifold  shown  to  the  manometer. 

A  1/4  Inch  steel  ball  was  partially  enfolded  In  the 
top  of  the  manifold.  The  manifold  was  spring  loaded  to  bear 
on  a  micrometer  splnder  and  barrel  which  was  mounted  In  a  pedes¬ 
tal  bolted  to  the  base  plate.  Rotary  motion  of  the  barrel  caused 
an  axial  motion  of  the  spindle  which  moved  the  manifold  through 
contact  at  the  steel  ball. 

• 

4. 3. 1.2  Determination  of  the  Radial  Position  of  the  Total 
Pressure  Probe  In  the  Gap 

To  determine  the  probe  distance  from  the  rotor  or 
stator  at  a  given  micrometer  setting  the  probe  was  brought 
close  to  the  rotor  or  stator  until  contact  was  made  using  a 
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0.001  Inch  piece  of  cigarette  paper  as  a  feeler  gage.  The 
center  of  the  probe  opening  relative  to  the  probe  height 
was  known,  the  measurements  having  been  taken  with  the 
toolmaker '8  microscope.  In  this  manner  It  was  possible 
to  determine  the  distance  from  the  rotor  or  stator  surface  to 
the  center  of  the  probe  opening. 

The  probe  was  used  to  determine  the  change  In 
rotor  radius  with  temperature.  With  the  rotor  stationary, 
and  heated,  the  following  relationship  was  determined: 
rrotor  -  1*660  +  0.000563(MV)  -  0.000294 

where  MV  Is  the  millivolt  reading  from  the 
rotor  thermocouple  connected  to  the  rings  furthest  from  the 
rotor . 


A  plot  of  the  data  taken  to  obtain  this  relation 
is  shown  in  Pig.  4-6. 

4. 3. 1.3  Total  Pressure  Probe  Corrections 

Total  pressure  probes  are  subject  to  several  effects 
which  may  cause  errors  In  measurement. 

4. 3. 1.3.1  Viscous  Effects 

So  long  as  the  Reynolds  number  based  on  a  probe 
dimension  such  as  radius  or  height  Is  sufficiently  large 
(for  example  >  100)  the  Indication  of  a  total  pressure  probe 
will  not  be  Influenced  by  viscous  effects. 


&</?<// j/v 
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For  Reynolds  numbers  between  30  and  100,  the 
Indication  will  be  about  0.99  the  true  total  pressure. 

For  Reynolds  numbers  less  than  30,  viscous  forces 
play  a  more  dominant  role  and  the  measured  Impact  pressure 
Increases  as  the  Reynolds  Number  decreases. 

l  2 

This  effect  has  been  studied  analytically  ' 
and  experimentally^*  4 • 

Hurd,  Che sky  and  Shapiro  recommend  from  experiments 

/f>c6cei£cd 

1 

(&</<>) 

The  results  of  MacMillan,  for  probes  of  different 


Barker,  M.  "On  the  Use  of  Very  Small  Pitot  Tubes  for  Measuring 
Wind  Velocity",  Proc.  Royal  Soc.  London,  A101,  (1922),  pp.  435-445. 

2R.  W.  Ladenburg,  B.  Lewis,  R.  N,  Pease,  H.  S.  Taylor,  Physical 
Measurements  In  Qas  Dynamics  and  Combustion  (Princeton.  New 
Jersey:  Princeton  University  Press,  1$$4),  pp.  112-118. 

^C.  W.  Hurd,  K.  P.  Chesky,  and  A.  H.  Shapiro,  "Influence  of 
Viscous  Effects  on  Impact  Tubes",  J.  Appl.  Mech.,  Vol.  20, 

(1953),  PP.  253-256. 

^F.  A.  MacMillan,  "Viscous  Effects  on  Flattened  Pitot  Tubes 
at  Low  Speeds",  Aeronautical  Research  Council,  Fluid  Motion 
Sub -Committee,  FM208la,  Oct.  1954. 
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width  to  height  ratios  Indicate  a  Reynolds  Number  effect 
starting  at  about  Q,  >20.  The  data,  of  MacMillan  are  shown 
In  Pig.  4-7.  For  a  width  to  Internal  height  ratio  of  about 
4  to  1  as  for  the  present  probes,  the  data  show  a  correction 
of  2#  for  a  Reynolds  Number  of  9-  This  corresponds  to  a 
velocity  of  about  4  ft/second  in  atmospheric  air.  No  velocities 
this  low  were  measured  during  the  present  experiment. 

4. 3. 1.3. 2  The  Effect  of  Transverse  Shear  and  Wall  Interference 

When  a  total  pressure  probe  is  placed  in  a  flow  in 
which  there  is  a  transverse  velocity  gradient  the  indication 
given  by  the  probe  Is  the  total  pressure  not  at  the  geometric 
center  of  the  probe  but  at  some  "effective  center",  displaced 
from  the  geometric  center  toward  the  higher  velocity  because 
the  presence  of  the  probe  causes  a  distortion  of  the  probe 
streamlines  toward  the  region  of  lower  velocities. 

• 

In  a  boundary  layer^wall  interference  effects 
further  confuse  the  situation. 

Some  experimental  work  has  been  done  on  the  "dis¬ 
placement"  correction  applicable  to  total  pressure  probes  in 
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turbulent  boundary  layers.* 

Llvesey^*  3  has  apparently  studied  the  displacement 
corrections  for  flat-end  pitot  tubes  In  a  turbulent  boundary 
layer  after  It  had  been  thickened  by  a  wire  cloth  screen. 

Llvesey  states  that  the  displacement  correction  Is 
a  function  of  the  solidity  of  the  cross-section  and  the  aspect 
ratio  of  the  probe  tip.  He  further  comments  "So  far,  work  on 
compressible  sub-  and  supersonic  flow  has  not  given  any  consist¬ 
ent  indication  of  that  [displacement]  error”. 


A.  D.  Young  and  J.  N.  Maas,  in  "The  Behaviour  of  a  Pitot 
Tube  In  a  Transverse  Total  Pressure  Gradient,"  Aeronautical 
Research  Council,  R.and  M  1770  1936;  I.  M.  Hall  in  "The 
Displacement  Effect  of  a  Sphere  in  Two  Dimensional  Shear 
Plow",  Journal  of  Fluid  Mechanics.  Vol.  1,  Part  2,  (July,  1956), 
142-162,  and  M.  J.  Llghthill  in  "Contributions  to  the  Theory 
of  Pitot -Tube  Displacement  Effect",  Journal  of  Fluid  Mechanics, 
Vol.  2,  Part  5,  (July,  1957) ,  pp.  493-512,  concerned  themselves 
with  the  displacement  corrections  applicable  to  flow  In  a  wake. 
The  corrections  described  are  not  applicable  here.  This  is 
mentioned  because  In  the  literature  one  finds  the  Young  and 
Maas  experimentally  determined  displacement  correction 
Incorrectly  applied  to  boundary  layer  total  pressure  traverses. 

2 

J.  L.  Llvesey,  "Design  of  Total  Pressure  Probes  for  Minimum 
Interference  with  Measured  Flow,"  J.  Aero.  Sci.,  Vol.  21, 

(Sept.,  1954),  p.  641. 

3 

Unpublished  work  by  Llvesey  is  referred  to  In  the  work  by 

P.  0.  A.  L.  Davies,  "The  Behaviour  of  a  Pitot  Tube  In  Transverse 
Shear",  Aeronautical  Research  Council,  (December,  1957)*  P.M. 
2619,  Fluid  Motion  Sub -Committee,  (December,  1957) • 
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F.  A.  Macmillan1  working  with  round  pitots  reports 
displacements  of  the  order  of  0.15D,  (where  D  Is  the  probe 
height).  In  a  pipe  and  a  boundary  layer  on  a  flat  plate. 

Davies^  found  experimentally  that  for  a  flat -ended 
round  pitot  of  conventional  form  no  displacement  corrections 
were  necessary  In  a  turbulent  boundary  layer  even  when  the 
pitot  was  touching  the  wall. 

Davies  also  suggests  that  because  of  the  wire  cloth 
screen  Llvesey  used  to  thicken  his  turbulent  boundary  layer, 
the  conditions  Llvesey  studied  were  more  nearly  like  those  In 
a  wake  rather  than  In  a  turbulent  boundary  layer. 

Deissler^  made  velocity  profile  measurements  In  a 
tube  with  both  round  and  flattened  probe  tips  In  the  vicinity 
of  the  wall.  Because  no  appreciable  difference  could  be 
noticed  between  the  measurements  made  with  the  two  types  of 
tips,  he  concluded  that  the  aerodynamic  and  geometric  centers 
of  the  probe  coincided. 


^F.  A.  Macmillan,  "Experiments  on  Pitot  Tubes  In  Shear  Flow", 
Aeronautical  Research  Council,  R  and  M  3028,  (Feb.,  1956). 

2Davies,  ARC,  F.M.  2619,  December,  1957. 

3 

R.  G.  Deissler,  "Analytical  and  Experimental  Investigation  of 
Adiabatic  Turbulent  Flow  In  Smooth  Tubes",  NACA,  TN  2138, 
(July,  1950). 
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The  Staff  of  the  Aerodynamics  Division,  National 
Physical  Laboratory  of  the  Ministry  of  Aviation1  In  establishing 
a  calibration  curve  for  skin  friction  determination  on  a 
flat  plate  with  a  Preston  tube  (a  round  pitot  attached  to 
the  surface)  show  a  curve  (their  figure  11)  for  tyt/d vs.  & 
(where  j f/e  is  ^effective  center  of  the  pitot  and 
the  diameter  of  the  tube).  They  show  that  t/cfo ,  although 
varying  with  Reynolds  Number,  has  a  value  about  0.62,  and 
also  that  ye  Is  displaced  toward  the  higher  velocities. 

o 

Matting,  Chapman,  Nyholm  and  Thomas  ,  from  measure¬ 
ments  of  velocity  profiles  and  separate  measurements  of  skin 
friction  on  a  flat  plate  at  high  Mach  numbers,  report  an 
"Increasing  discrepancy  [with  Increasing  Mach  numbers]  between 
the  slope  deduced  from  skln-frlctlon  measurements  and  the 
apparent  slope  of  the  measured  profiles",  and  attribute  this 
to  a  displacement  effect  caused  by  the  probe.  They  note  that 

the  apparent  distortion  is  of  the  same  order  of  magnitude  as 

•a  Ji 

reported  by  Taylor-'  and  Von  Doenoff  . 


^Staff  of  the  Aerodynamics  Division,  N.P.L.,  Ministry  of 
Aviation,  "On  the  Measurement  of  Local  Surface  Prlction  on  a 
Plat  Plate  by  Means  of  Preston  Tubes",  Aeronautical  Research 
Council,  R.  and  M.  No.  3185,  May  1958,  Published  1961. 

2F.  W.  Matting,  D.  R.  Chapman,  J.  R.  Nyholm  and  A.  0.  Thomas, 
"Turbulent  Skin-Frlctlon  at  High  Mach  numbers  and  Reynolds 
Numbers  in  Air  and  Helium",  NASA  Technical  Report  R-82,  1961. 

^G.  I.  Taylor,  "Measurements  with  a  Ha If -Pitot  Tube",  Proc. 

Roy.  Soc.  (London),  Ser.  A,  Vol.  166,  1938,  pp.  476-481. 

4 

A.  £.  Von  Doenhoff,  "Investigation  of  the  Boundary  Layer  About 
a  Symmetrical  Airfoil  In  a  Wind  Tunnel  of  Low  Turbulence, 

NACA  WR-L-507,  1940. 
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The  velocity  profiles  for  the  present  experiment 
also  Indicate  a  displacement  effect.  The  method  used  for 
correcting  the  date  for  this  effect  Is  given  In  Section  6.4. 

i 

4. 3. 1.3 .3  The  Effects  of  Turbulent  Fluctuations 


The  lag  of  the  total  pressure  probe  and  micromanometer 
causes  the  probe  to  indicate  some  average  value  of  the  total 
head.  The  total  head  Includes  not  only  the  dynamic  head  of 
the  mean  velocity  V&  but  also  the  effect  of  the  turbulent 
velocities  ^  l Vg  . 


Goldstein1  for  an  Incompressible  fluid,  suggests 
that  the  turbulent  terms  Included  In  the  reading  of  the  total 
head  probe  be 


4U,  -  A*; 


Hlnze  and  Van  der  Hegge  Zlynen2  as  reported  by  Hlnze3 
neglected  the  effects  of  the  lateral  turbulence  components 
because  of  lack  of  reliable  data. 


Alexander,  Baron,  and  Comings1*’  observed  that  the  total 


*S.  Goldstein,  Modern  Developments  In  Fluid  Mechanics  (1st  ed; 
Oxford  at  the  Clarendon  Press,  1938)*  I*  p.  253. 

2 

J.  0.  Hlnze  and  B.  G.  Van  Der  Hegge  Zlynen;  Applied  Science 
Research,  Vol.  1A,  p.  435  >  1949* 

O 

j.  0.  Hlnze,  "Turbulence,"  (1st  ed;  McGraw-Hill  Book  Company, 
Inc.,  New  York,  I'95£77pp.  132-137. 

^L.  G.  Alexander,  T.  Baron,  and  E.  W.  Comings,  University  of 
Illinois  Experimental  Station,  Technical  Report  No.  8,  1950. 
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head  tube  reading  decreased  with  an  Increase  In  relative 
turbulent  Intensity  (contrary  to  expectations),  so  that  they 
neglected  the  effect  of  turbulence  completely  and  used 

Hinze1  states  that  "Neither  of  these  two  methods 
Is  entirely  correct  but..»at  present  no  better  method  can  be 
suggested".  In  Section  6-2  Is  presented  an  analysis  which 
shows  that  the  terms  neglected  in  the  equation  above  for  the 
present  case  are  1.  Interaction  between  molecular  and  turbulent 
terms  and  2.  terms  such  as  which  Include  second 

and  third  order  correlation  and  are  much  smaller  than 
everywhere  In  the  flow.  The  equation  used  as  representing  the 
reading  of  the  total  head  probe  Is 

4.3.2  Static  Pressure  Measurement 

For  the  purpose  of  determining  the  axial  static 
pressure  distribution  at  the  stator  wall,  seven  0.010  Inches 
static  pressure  taps  were  provided  as  described  In  Section 
3.3.1. 


4. 3. 2.1  Errors  In  Indication  Caused  by  Pressure  of  Tap 

The  stainless  steel  tubing  from  which  these  taps 


J.  0.  Hlnze,  Turbulence . 


p.  136. 
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were  made  were  countersunk  0.005  Inches  to  eliminate  the 
effects  of  disturbances  to  the  static  pressure  at  the  atator 
wall  caused  by  the  presence  of  the  tap,  whereby  the  flow 
streamlines  tend  to  dip  into  the  tap  hole. 

l  2 

Experiments  conducted  by  R.  E.  Rayle  *  show  that 
a  tap  diameter  0.010  Inches  to  0.040  Inches  countersunk  to 
a  depth  of  1/2  the  diameter  will  substantially  eliminate  tap 
errors . 

4.3.3  Micromanometer  for  the  Indication  of  Total  and 
Static  Pressures. 

The  measurements  of  static  and  total  pressures  to 
within  0.001  inches  of  water  were  made  with  a  micromanometer 
from  a  design  by  A.  M.  0.  Smith  and  James  S.  Murphy^.  The 
pressure  to  be  determined  lsled  through  Tygon  Tubing  to  an 
inclined  sight  tube  which  forms  one  leg  of  a  U-table  manometer. 
The  other  leg  is  brought  to  a  reservoir  which  Is  mounted  on 
the  arm  of  a  Vernier  height  gage. 


^■R.  E.  Rayle,  "An  Investigation  of  the  Influence  of  Orifice 
Geometry  on  Static  Pressure  Measurements",  MIT,  SM  Thesis, 
Dep't.  of  Mech.  Eng'g.,  (1949). 

p 

R.  C.  Dean,  Jr.,  Aerodynamic  Measurements  (Boston;  Eagle 
Enterprises,  1953)*  P*  od. 


^A.  M.  0.  Smith,  and  J.  S.  Murphy,  "Micromanometer  for  Measuring 
Boundary  Layer  Profiles",  Rev.  Scl.  Inst's..  Vol.  26,  Number 
8  (Aug.,  1955).  —  — 
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During  operation,  the  reservoir  Is  moved  up  or  down 
so  that  the  meniscus  In  the  sight  tube  Is  brought  to  Its  zero 
position.  This  zero  position  Is  determined  with  a  40  power 
microscope  fitted  with  a  reticule.  The  Inclined  tube  allows 
the  vertical  height  of  the  meniscus  to  be  determined  with 
the  microscope  to  well  within  0.001  Inches  of  fluid.  The 
base  of  the  Instrument  Is  fitted  with  levelling  bubbles  and 
levelling  screws,  and  the  reservoir  Is  placed  as  close  to 
the  sight  tube  as  possible.  Thus,  two  Important  sources  of 
error  in  U-tube  manometers  are  eliminated^  /hat  caused  by 
Inaccuracy  in  the  determination  of  the  meniscus  level  and  that 
resulting  from  tilt  between  the  legs  of  the  manometer. 

Experiments  by  Smith  and  Murphy  to  determine  a 
proper  manometer  fluid  with  regard  to  a)  adhesion  of  the 
manometer  fluid  to  the  walls  of  the  tube,  b)  low  volatility 
and  c)  stable  density,  i.e.,  minimum  tendency  to  absorb  moisture 
from  the  air,  resulted  in  a  choice  of  DC-200  (J>CS)^  a  silicone 
fluid  manufactured  by  the  Dow-Corning  Corporation. 

MacMillan^  has  suggested  that  to  minimize  manometer 
zero  shift  caused  by  ambient  temperature  change  the  expansion 
of  the  reservoir  and  manometer  fluid  be  balanced  against  the 
change  in  the  meniscus  position  caused  by  change  In  surface 
tension. 


^Dow-Coming  DC2OO(0CS);  density  0.818  gms/cm?,  boiling  point 
152  deg.  C,  surface  tension  17.0  dynes/cm  and  viscosity  0.818 
centpolses. 

2P.  A.  MacMillan,  Liquid  Manometers  with  High  Sensitivity  and 
Small  Time  Lag",  J.  Scl.  Instr.,  V.  31,  1954#  p.  17 • 
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This  was  accomplished  by  Smith  and  Murphy  by  proper 
design  of  the  component  tubing  and  reservoir. 

Smith  and  Murphy  found  that  the  micromanometer  had 
an  accuracy  of  0.0002  Inches  due  to  random  adherence  of 
fluid  to  the  tubing  walls.  However,  the  accuracy  of  the 
present  manometer  is  0.001  Inches,  the  accuracy  of  the  vernier 
height  gage  used  to  position  the  reservoir. 

4.4  Temperature  Distribution  in  the  Gap  Between  the  Two 
Cylinders 

4.4.1  Introduction 

The  temperature  distribution  In  the  gap  between 
the  cylinders  was  determined  from  the  density  distribution  as 
measured  with  a  Mach  Zehnder  Interferometer. 

An  interferometer  was  chosen  for  this  purpose 
because  as  an  optical  method  it  has  several  advantages  over 
mechanical  means  such  as  a  thermocouple  probe. 

The  advantages  of  an  interferometer  are  as  follows: 

a)  its  use  does  not  disturb  stream  conditions 
of  temperature  and  flow,  as  would  the  use  of  a  probe 

b)  its  use  permits  the  entire  flow  field  to 
be  studied  from  a  single  interferogram. 

It  should  be  noted  that  in  addition,  a  probe  would 
have  been  difficult  to  calibrate  because  of  radiation  effects 
at  the  higher  temperatures. 
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The  disadvantages  of  the  Interferometer  are  that  It 
Is  costly,  difficult  to  build,  adjust  and  align.  Also,  Its 
use  In  the  present  experiment  Introduced  a  small  error  due 
to  quartz-window  heating  effects. 

The  theory  and  application  of  the  Mach-Zehnder 

12^ 

Interferometer  has  been  extensively  covered  In  the  literature  *  ,J 
In  Appendix  I  Is  given  a  brief  description  of  the  fundamentals 
of  Its  operation  and  adjustment. 

4.4.2  The  Design  of  the  Mechanical  Parts 

The  principle  mechanical  components  of  the  Mach- 
Zehnder  Interferometer  under  discussion  are 

1.  The  base  on  which  the  mirror  holders 

are  mounted 

2.  The  mirror  holders 

3 .  The  holders 

4.  The  light  source,  collimating  lens,  and 

camera. 

^See  for  example  the  extensive  bibliographies  in  Ladenburg,  et  al, 
pp.  75-78  and  Dean,  pp.  270-271. 

2L.  H.  Tanner,  "The  Optics  of  a  Mach-Zehnder  Interferometer", 
Aeronautical  Research  Council,  R  and  M  3069,  Oct.  1956* 

^L.  H.  Tanner,  "The  Design  and  Use  of  Interferometers  In  Aero¬ 
dynamics",  Aeronautical  Research  Council,  R  and  M  No.  3131,  Sept. 
1957. 
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A  photograph  of  the  entire  assembly  Is  shown  in  Pig.  4-F. 

4. 4. 2.1  The  Base 

A  photograph  of  the  base  plate  on  Which  the  mirrors 
are  mounted  Is  shown  In  Fig.  4-  $. 

The  base  consists  of  a  1/2  Inch  thick  steel  plate 
approximately  30  inches  square  to  the  back  of  which  have  been 
welded  I-beam  sections  In  a  U-arrangement  as  shown. 


Pour  blocks  of  steel  1  Inch  by  6  Inches  by  4  Inches, 
their  centers  at  the  comers  of  a  rectangle  20  Inches  by  10 
Inches,  and  oriented  at  45  degrees  to  the  legs  of  the  rectangle 
were  welded  to  the  1/2  Inch  thick  steel  plate  and  machined 
flat  on  a  Jig -borer.  The  Interferometer  mirror  holders  were 
bolted  onto  these  blocks.  No  stress  relief  was  provided  other 
than  that  of  allowing  the  base  to  stand  at  room  temperature 
for  about  a  month. 

The  rigidity  of  the  base  so  designed  has  proved 
satisfactory.  Indeed,  the  only  precaution  needed  against 
distortion  of  the  base  Is  that  it  not  be  subjected  to  temper¬ 
ature  gradients.* 

The  base  is  bolted  at  three  points  to  a  triangular 
plate,  and  this  plate  is  hung  at  its  comers  through  vibration 

*Dean,  p.  198  reports  that  for  a  similar  interferometer  at 
M.I.T.  a  temperature  variation  of  0.1  deg.  P  in  one  arm  of  the 
base  displaces  the  bands  about  one  band  space,  but  observes 
that  the  principle  precaution  required  is  not  to  allow  the 
base  to  be  subject  to  direct  sunlight. 
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figure  4-0 

Interferometer  Assembly 


Figure  4-9 


Interferometer  Base-Plate 
and 

Mirror  Holders 
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Isolators  from  a  f raise  made  of  2  Inch  pipe.  This  frame  has 
been  sufficient  to  prevent  major  disturbances  to  the  mirrors 
and  plates  once  they  have  been  adjusted  to  produce  fringes. 

The  pipe  frame  Is  mounted  on  wheels. 

The  arrangement  described  has  been  satisfactory 
since  from  day  to  day  and  month  to  month, only  minor  adjustments 
to  fringe  spacing  are  necessary.  Indeed,  the  fringes  have 
disappeared  but  once  accidentally,  when  a  mirror  holder 
Inadvertantly  was  given  a  direct,  sharp  blow. 

4. 4. 2. 2  The  Mirror  Holders 

The  mirrors  must  be  mounted  so  that  they  are  free 
to  rotate  in  any  plane  and  also  so  that  they  are  stress  free. 

A  gimbal-type  mounting  is  provided  to  meet  the  first  condition. 
Such  a  mounting  is  shown  in  Pig.  4-10.  The  inner  mount,  in 
which  lies  the  2  inch  by  4  inch  mirror  is  made  3  or  4  thousandths 
of  an  inch  larger  on  a  side  than  the  mirror.  Scotch  drafting 
tape  is  used  as  a  gasketing  material  between  the  mirror  and  the 
frame.  Scot^A drafting  tape  is  ideal  for  this  application  since 
its  wrinkled  surface  allows  for  expansion  and  contraction  of 
the  mirror  under  dulmal  temperature  change  and  yet  provides 
sufficient  drag  so  that  the  mirrors  do  not  shift  about. 

On  each  side  of  the  mirror  is  placed  a  rectangular 
metal  frame  which  fits  over  the  mirror  holder  and  is  attached 
thereto  with  screws.  Bach  metal  frame  has  three  protrusions 
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which  extend  oyer  the  mirror  and  serve  as  supports.  Turning 
the  screws  by  whioh  the  frame  is  attaohed  to  the  mirror  holder 
also  allows  rough  positioning  of  the  mirrors  as  needed  during 
adjustment.1  The  small  rotations  required  are  accomplished 
on  the  present  apparatus  by  having  a  flat  spring  attached  to 
the  frame  loaded  by  a  screw  of  small  lead  (0.025  Inches  per 
revolution).  This  load  creates  a  moment  which  tends  to  twist 
the  pieces  on  which  the  frame  is  supported.  The  system  is  so 
designed  that  one  revolution  of  the  screw  rotates  the  frame 
30  x  10  radians.  Minute  motions  are  provided  very  simply, 
as  contrasted  with  complicated  gear  trains. 

4. 4. 2. 3  Mirrors 

The  plates  and  mirrors  used  are  made  from  selected 
optical  quality  plate  glass  and  were  purchased  from  The  Labo¬ 
ratory  Optical  Company,  Plainfield,  New  Jersey.  They  are  2 

Inches  by  4  inches  by  3/4  Inches  thick  and  are  flat  to  within 

2 

0.1  wavelength  of  mercury  green  light. 

4. 4. 2. 4  The  Light  Source, Collimating  Lens,  and  camera 

The  light  source  used  is  a  General  Electric  B-H6 

*If  the  spacing  between  the  fringes  is  to  be  one  millimeter, 
so  that  they  can  be  seen  with  the  naked  eye,  and  if  mercury 
green  light  is  used  then  the  angle  of  rotation  of  a  mirror 
to  produce  the  fringes  is 

i  -  ?461  x  10  8  a  2730.5  x  10“7  radians 

2d  2  x  10"1 

^Verbal  communication,  from  Mr.  S.  Brower,  Laboratory  Optical 
Company. 
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air-cooled  Mercury  vapor  lamp  in  continuous  operation  or 
flashed,  as  needed.  The  rating  of  this  lamp  is  as  follows: 
Watts,  900;  Lumens,  60,000;  starting  voltage  (AC),  1200; 
Operating  voltage,  900;  starting  amperes,  2.6;  operating 
amperes,  1.25. 

The  percentage  of  the  total  energy  output  of  the 

* 

lamp  in  the  wavelength  range  5388-5536  A  is  16.8.  The 
principle  line  in  this  region  is  5461  A. 

A  short  distance  from  the  lamp  was  mounted  a  piece 
of  brass  shimstock  0.0003  Inches  thick  in  which  was  drilled 
a  hole  0.004  inches  in  diameter.  This  hole  was  placed  at  the 
focal  point  of  an  f/2.5  Aero  Ektar  lens,  and  thus  served  as 
the  point  source . 

Between  the  source  and  the  lens  was  placed  a  multi¬ 
layer  interference  filter1  which  passed  the  mercury  green  line, 

5461A,  with  approximately  70$6  transmission  and  half -band 

# 

width  of  approximately  70  A. 

With  this  arrangement  of  light  source,  operating 
continuously,  0.004  inch  hole  filter,  and  collimating  lens, 

1/20  second  exposure  time  was  required  to  expose  a  103  Q 
Eastman-Kodak  spectroscopic  plate  with  1/2  magnification  of 
the  test  section. 

For  flashing  operation  of  the  lamp,  a  2  microfarad 

^-Manufactured  by  Baird  Associates  -  Atomic  Instrument  Company, 
33  University  Road,  Cambridge  38,  Massachusetts.  See  Technical 
Circular  HD  503. 
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capacitor  at  3000  volts  was  discharged  through  the  lamp.  A 
circuit  diagram  of  the  arrangement  Is  shown  In  Pig.  4-11. 

One  flash  of  unflltered  light  was  sufficient  to  expose  an 
Eastman-Kodak  Super-Ortho  press  photographic  plate,  when  an 
opening  of  0.010  Inches  was  used  as  the  point  source.  For 
filtered  light  four  flashes  of  the  lamp  were  required  to 
expose  a  similar  plate. 

To  photograph  the  Interference  bands,  another  Aero- 
Ektar,  f/2.5,  7  Inch  lens  was  mounted  on  a  4  Inch  x  5  Inch 
Speed  Graphic  camera  box. 

4.4.3  Corrections  Required 

4. 4. 3.1  The  Refraction  Error 

Because  a  density  gradient  exists  In  the  gap  a 
light  ray  passing  axially  through  the  gap  Is  refracted  away 
from  the  heated  rotor. 

On  the  photographic  plate  a  refracted  ray  appears 
to  originate  at  the  same  position  In  the  gap  as  an  unrefracted 
ray  which  would  have  struck  the  photographic  plate  at  the 
same  point. 

There  is,  therefore,  an  error  In  fringe  shift  as 
read  on  the  photographic  plate  occasioned  by  the  difference 
in  optical  path  lengths  between  the  refraeted  ray  and  the 
unrefracted  ray. 


figure.  4^f! 

fileofrica/ 6'recttt  £r  Sn^UFfash 
of-  8M~C*  fiercctry 
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1  2 

Wachtell *  has  shown  that  to  a  first  approximation 
the  error  vanishes  If  the  camera  Is  focused  at  a  plane 
located  at  an  axial  dlstanoe  */*-  -  1/3,  measured  from  the  exit 
window;  that  Is  the  window  nearest  the  camera.  Further,  the 
error  In  using  the  first  approximation  was  found  to  be  less 
than  5  per  cent  of  the  density  difference  between  the  wall 
and  free  stream  conditions  for  a  boundary  layer  on  a  flat 
plate  If 

$  9  for  a  turbulent 

boundary  layer. 

In  this  equation  -  span  of  the  tunnel, 

&  m  boundary  layer  thickness 
“  density  difference  between 
the  free  stream  and  the  wall. 

For  the  most  severe  situation  In  the  present  study, 
7000  rpra,  Trotor  -  890  F,  conservatively  estimating  £  to  be 


l  « 

R.  Ladenburg,  and  D.  Bershader,  "Optical  Studies  of  Boundary 
Layer  Phenomena  on  a  Flat  Plate  at  Mach  Number  2.35,"  Final 
Technical  Report  for  Contracts  NR061  and  N6ori-105  Task  II, 
Palmer  Physical  Laboratory,  Princeton,  New  Jersey,  Dec.  15, 

1952.  See  Chapter  VIII,  The  Refraction  Error  by  Q.  P.  Wachtell. 


G.  P.  Vachtell,  "Refraction  Error  In  Interferometry  of  Boundary 
Layer  In  Supersonic  Flow  Along  a  Flat  Plate",  Fh.D.  Dissertation 
Princeton  University,  1951 • 
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0.040  inches,/^  -  0.0547,  fi,  -  0.0446  and/C  -  0.3635  x  10"2. 

Thus _ 

C-.  ffefrfij’ 

■=0.68  ^ 

and  the  criterion  Is  satisfied.  In  fact,  the  error  decreases 

4 

as  C  so  that  the  error  here  Is  negligible. 

The  camera  was  therefore  focused  In  the  plane 
X/l.  =  1/3  during  the  present  experiments. 

Another  approach  to  the  refraction  error  problem 
is  given  by  DeFrate*  and  is  useful  when  the  criterion  of 
Wachtell  is  not  satisfied. 

DeFrate  suggests  an  iteration  procedure  as  follows: 

A  density  distribution  is  measured  from  the  inter¬ 
ferometric  plate.  Rays  are  traced  through  this  density  dis¬ 
tribution.  Noting  the  angle  at  which  they  leave  the  test 
section,  one  can  calculate  the  y-position  in  the  gap  from  which 
they  appear  to  come  to  the  photographic  plate.  Further,  for 
each  ray  the  change  in  optical  path  length  from  a  no-flow 
condition  can  be  calculated.  For  each  ray  one  obtains  a 
density  and  a  y-position. 


1L.  A.  DeFrate,  "Application  of  the  Interferometer  to  the 
Study  of  Boundary  Layers,  Sc.D.  Dissertation,  M.I.T.,  1950* 


By  traoing  rays  a  density  distribution  oan  be  obtained. 

This  Is  compared  with  the  measured  distribution.  If  no  agree¬ 
ment  is  obtained,  rays  are  traced  through  the  caloulated  density 
distribution  and  the  process  repeated  until  agreement  Is 
achieved  between  the  assumed  density  distribution  and  the  one 
calculated  from  ray  tracing. 

4.4. 3.2  The  Error  Due  to  Temperature  Gradient  In  the 
Quartz  Windows 

Consider,  In  Fig.  4-12,  a  drawing  of  one  of  the 
quartz  windows  over  the  ends  of  the  gap.  These  windows  pro¬ 
vide  a  path  for  the  Interferometric  test  beam. 

During  the  experiment  the  window  experiences  a 
temperature  gradient  since  one  end  "sees"  the  ends  of  the  hot 
rotor,  and  the  other  Is  In  contact  with  the  cooling  water 
Jacket.  The  quartz  window  under  goes  a  varying  change  In 
thickness  across  the  light  path  as  well  as  a  change  in  Its 
index  of  refraction.  Thus,  the  optical  path  through  a  region 
In  the  quartz  In  the  center  of  the  gap  differs  from  that 
through  a  region  near  the  walls  and  an  error  In  Indicated 
density  results.  The  error  can  be  evaluated  as  follows: 

Consider  the  quartz  window  show  In  Fig.  4-12.  Under  a  temper¬ 
ature  gradient  expansion  occurs  as  shown  In  the  figure. 

The  optical  path  length  at  (1)  In  the  oold  condition 
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hot  condition 

*  mAjH  w  c  Cf  +*?  6  f 

A%" *  ^ 

same  can  be  written  as 

=  'dz^c^c 

and 

*^0*- *?&)(**,* Xi+p-  **0 

The  difference  in  optical  path  lengths  between  the 
locations  (1)  and  (2)  in  the  hot  condition,  the  optical  path 
length  of  interest  is 

AO  PL  =  d2j^f /  -L 

.*"2 vc^/  -t*  c*jAr7i)(7l<,f^ f  A /,^ 

Since 

/  ' 

and 

^•a^c  * 

-  4  -Trffo 


and  In  the 

CO^Ma- 
At  (2)?the 
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/ 

For  the  values 


Jn 

dr 


-  0.32  lnchel,  X  10-6 
lnch-deg  F 

■  6.272  per  deg  F  x  10"^ 

«  0.125  Inches 


71c.  «  1.46019  at  75  deg  F. 


AOfiL  -  O.res  fA 460/9) fa* )J/ax+£Vz)/o‘+^7for, -ATj/o'J 


The  term  2.007  x  10"12  {AT-A71)  is  negligible. 
The  result  Is 

_5cy) 

Since  /J  »  5460  A  =  21.496  x  10"6  inches  y 


^/degree  temperature  difference  ■  0.055964 
If  the  temperature  difference  between  points  (1)  and  (2)  in 
the  hot  condition  is  known,  the  error  in  fringe  shift  as 
measured  on  an  interferogram  is  known,  and  can  be  corrected. 


In  Appendix  II  is  outlined  the  test  procedure  used 
to  determine  the  temperature  gradient  in  the  glass.  Tenpllaq 

International  Critical  Tables  (first  ed.,  McGraw-Hill  Book 
Company,  New  York,  I929J  VoTT  IV,  p.  21,  Vol.  VI,  p.  341,  p.  342. 
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'was  uaed  to  determine  the  Isotherms  on  the  glass.  The 
experiments  showed  the  Isotherms  to  be  approximately  con¬ 
centric  circles  except  for  the  corners  of  the  glass. 

The  final  corrections  applied  (see  Section  6. 3. 1.2.1) 
are  small  over  small  distances  from  the  rotor  (the  region  of 
most  Interest)  but  are  applied  In  the  Interest  of  consistency 
of  the  data  obtained  from  the  Interferometer  and  the  rotor  and 
stator  thermocouples. 

4.5  The  Use  of  the  Hot-Wire  Anemometer 

A  hot-wire  anemometer  purchased  from  the  Plow 
Corporation  of  Cambridge,  Massachusetts  was  used  to  study 
the  turbulence  characteristics  of  the  flow. 

A  probe  consisting  of  only  one  wire  was  placed  in 
the  stream  so  as  to  measure  the  &  or  r-component  of  the  turbu¬ 
lent  fluctuations. 

Operation  and  use  of  the  hot-wire  anemometer  is 

1  2 

described  In  detail  elsewhere  *  and  specifically  In  the 
literature  distributed  by  the  Plow  Corporation. 


R.  W.  Ladenburg,  B.  Lewis,  R.  N.  Pease,  and  H.  S.  Taylor, 
Physical  Measurements  in  Gas  Dynamics  and  Combustion.  (1st  ed; 
Princeton  University  Press,  Princeton,  195^),  Article  P,  2  by 
L.  S.  G.  Kovasznay,  pp.  219-276. 


2  . 

See  also  the  bibliography  on  pp.  283-285  In  the  same  book 
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Chapter  V 

Experimental  Procedures 


5.1  Introduction 

The  apparatus  was  operated  at  3000,  5000  and  7000 
rpm,  with  rotor  temperatures  of  about  5^0  R,  600  R,  720  R,  and 
890  R.  These  conditions  correspond  to  Taylor  numbers  of  from 
2500  to  9000,  and  temperature  differences  of  0  P  to  330  P. 

This  represents  an  extension  in  Taylor  number  and  temperature 
difference  over  previous  investigators.  See  Tabe  1-1. 

The  experimental  procedures  followed  are  given 

below . 


5.2  Overall  Heat  Transfer 

For  the  measurement  of  overall  heat  transfer,  the 
data  required  are  rpm,  rotor  temperature,  stator  temperature, 
cooling  water  flow  rate,  and  temperature  rise  of  the  cooling 
water. 


After  the  start  of  rotation  and  the  application  of 
heater  current,  equilibrium  conditions  were  reached  in  about 
one  hour  at  lower  temperatures  and  two  hours  at  higher 
temperatures. 

Equilibrium  conditions  were  assumed  to  be  established 
when  three  successive  measurements  of  rotor  temperature  taken 
five  minutes  apart  were  within  1  deg  P . 


Ill 


The  cooling  water  flow  rate  and  teaperature  rise 
were  measured. 

Cooling  water  flow  rates  were  in  the  neighborhood 
of  30  lbs.  per  hour.  The  difference  in  teaperature  of  the 
cooling  water  ranged  froa  5  to  10  deg  F»  depending  on  the 
heat  input.  Water  flow  froa  the  apparatus  was  weighed  for  a 
period  of  about  twenty  minutes.  During  this  time  the  tempera - 
ture  difference  between  the  cooling  water  entering  and  leaving 
was  measured  every  two  minutes.  In  the  intervals,  water  flow 
rate  as  indicated  by  a  rotameter,  heater  power  and  rpa  were 
checked  for  constancy.  At  the  end  of  the  water  measurement 
period  rotor  and  stator  temperatures  were  again  measured. 

5.3  Measurement  of  Total  Pressure  Distribution 

For  the  runs  during  which  total  pressure  distribution 
in  the  gap  was  measured,  the  procedure  was  as  for  the  overall 
heat-transfer  runs.  When  equilibrium  conditions  were  reached 
and  the  rotor  temperature  determined,  the  graph  of  Fig.  4-6  was 
used  to  establish  the  radial  position  of  the  ppobe.  The  stator 
diameter  was  assumed  not  to  vary. 

During  the  early  velocity  distribution  runs,  at 
3000  and  especially  3000  rpm,  it  was  noticed  that  there  could 
exist  two  axial  pressure  distributions,  one  uniform  to  within 
2-1/2J6,  the  other  saw-toothed,  much  as  those  reported  by  Pal.1 


1  Pal,  Shih-I,  NACA  TN  892 
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A  typical  distribution  Is  shown  In  Pig.  5-1.  Tests  were  made 
with  a  hot-wire  anemometer  to  determine  the  characteristics  of 
the  flow.  It  was  found  that  there  could  exist  two  types  of 
flow  at  the  same  speed;  one  a  vortical  flow  of  the  type  reported 
by  Pai,  the  other  a  random  turbulent  flow.  These  results  are 
discussed  in  Chapter  VIII. 

The  uniform  axial  pressure  distribution  was  related 
to  the  conventional  random,  turbulent  flow  and  the  non-uniform 
axial  pressure  distribution  to  a  flow  in  which  Taylor  vortices 
as  discussed  by  Pai  predominated. 

At  5000  rpm  it  was  found  that  if  the  probes  were 
placed  near  the  rotor  surface  at  startup,  generally  a  uniform  ^ 

axial  pressure  distribution  would  result.  However,  at  times 
the  flow  would  "flip"  from  the  uniform  to  the  non-uniform  type 
during  a  run,  in  which  case  mechanical  disturbances  to  the 
flow  such  as  moving  the  probes  down  to  the  rotor  and  then 
to  the  stator  would  cause  a  return  to  the  uniform  type  of 
flow.  These  experiences  coincide  with  those  reported  by  Pai. 

At  3000  rpm,  simple  manipulation  of  the  probes  usually 
did  not  flip  the  flow  from  the  non-uniform  to  the  uniform  type. 

A  more  effective  method  was  a  sudden  increase  in 
rotor  speed  to  about  8000  rpm,  followed  by  a  gradual  reduction. 

This  did  not  succeed  at  all  times.  In  general,  obtaining  and 
maintaining  a  uniform  axial  pressure  distribution  was  difficult 


at  3000  rpm. 


i 
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At  7000  rpm,  the  non -uniform  type  of  distribution 
we s  never  encountered. 

There  wee  a  possibility  of  a  change  In  flow  during 
a  run,  which  could  go  undetected,  especially  while  the  probes 
were  In  the  middle  of  the  gap,  where  the  total  heads  for  the 
two  types  of  flow  were  almost  equal.  Therefore,  It  was  necessary 
to  measure  an  axial  pressure,  move  the  probe  to  a  new  location, 
remeasure  the  axial  pressure,  measure  the  total  pressure,  re¬ 
measure  the  axial  pressure  and  then  move  the  probe.  This  was 
especially  necessary  at  3000  rpm,  where  velocities,  axial 
pressures  and  total  heads  were  low.  At  5000  rpm,  two  or  three 
total  head  readings  were  taken  between  axial  pressure  checks 
when  the  probes  were  In  the  center  of  the  gap.  Closer  checks 
were  made  when  the  probes  were  near  either  wall. 

The  time  constant  of  the  probe  and  manometer  was  of 
the  order  of  eight  minutes,  and  of  the  axial  pressure  tubing 
and  manometer  of  the  order  of  five  minutes.  The  method  used 
for  measurement  was  to  bracket  the  reading  to  10.004  inches 
and  then  setting  the  well  to  the  mean  value. 

If  there  was  no  meniscus  motion  for  3  to  4  minutes 
the  reading  was  accepted.  If  there  was  motion  up  or  down,  the 
well  was  set  at  the  new  mean  value  and  the  procedure  repeated. 

The  average  time  for  one  reading  was  about  twelve 
minutes.  Including  the  static  pressure  ohecks. 


At  5000  rpm  the  average  time  for  a  run  Including 
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start  up,  was  about  eight  hours.  At  3000  rpm  the  average  time 
was  about  12  hours,  because  of  the  difficulties  of  maintaining 
a  uniform  flow. 

3. 4  Determination  of  Density  Distribution  with  a  Mach- 
Zehnder  Interferometer 

5.4.1  Alignment  of  the  Interferometer  with  the  Test  Piece 

The  method  for  obtaining  fringes  in  the  Interferometer 
is  given  in  Appendix^.  Before  use  further  alignment  is  required 
to  make  the  test  beam  of  the  Interferometer  parallel  to  the 
rotor  surface.  Methods  of  making  this  alignment  are  given  by 

12  7 

DePrate,  Wachtell  and  Howes  and  Buchele 

DePrate  and  Wachtell  suggest  using  the  diffraction 
pattern  caused  by  the  edges  of  the  test  piece.  Symmetry  of 
the  diffraction  pattern  indicates  that  the  test  beam  is  parallel 
to  the  rotor  surface.  Howes  and  Buchele  suggest  using  a  form 

4 

of  interference,  similar  to  Lloyds  mirror  interference,  which 


■^L.  A.  DePrate,  "Application  of  the  Interferometer  to  the  Study 
of  Boundary  Layers,"  Sc.D.  Dissertation,  M.I.T.,  1950* 

o 

0.  P.  Wachtell,  "Refraction  Error  in  Interferometry  of  Boundary 
Layer  in  Supersonic  Plow  Along  a  Plat  Plate,"  Ph.D.  Dissertation, 
Princeton, University,  1951* 

^W.  L.  Howes,  and  D.  R.  Buchele,  "Practical  Considerations  in 
Specific  Applications  of  Qas  Plow  Interferometry,"  HACA  TN  3507* 
July,  1955. 

^F.  A.  Jenkins  and  H.  E.  White,  Fundamentals  of  Optics  (2nd  ed; 
McGraw-Hill  Book  Company,  Inc.,  195°;,  pp.  230-239. 
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occurs  between  the  Reference  and  the  reflection  of  the  test 
beam  from  the  mls-allgned  surface.  The  results  of  this  Inter- 
ference  are  fringes  Inclined  to  the  main  Mach-Zehnder  fringes 
and  superimposed  on  them  at  the  surface.  When  the  surface  Is 
aligned  to  the  test  beam  this  Interference  disappears. 

The  method  of  symmetry  of  the  diffraction  pattern 
was  used  here  with  an  a  posterlorlcheck  using  the  methods 
suggested  by  Howes  and  Buchele.  The  method  was  applied  as 
follows : 


To  align  the  test  section  Initially  to  a  better 
approximation  than  obtained  by  eye,  a  frosted  100  watt. 
Incandescent  bulb  was  placed  at  the  end  of  the  test  section 
away  from  the  cameng  so  that  Its  rays  could  pass  through  the 
test  section.  Thus,  there  appeared  on  the  ground  glass  screen 
of  the  camera  a  black  and  white  image  of  the  gap  as  delineated 
by  the  rotor  and  stator  surfaces.  The  camera  lens  was  focused 
on  the  edge  of  the  rotor  nearest  to  It.  Then,  only  the 
test  beam  of  the  mercury  vapor  lamp  (with  filter)  was  allowed 
to  pass  through  the  test  section.  This  formed  a  green  Image 
on  the  ground  glass  screen.  When  the  white  Image  and  green 
Image  were  made  to  coincide  by  moving  the  test  section, 
alignment  was  near.  The  test  section  was  mounted  on  screw- 
jacks  for  this  purpose. 

The  camera  and  the  white  light  were  then  removed. 

The  diffraction  pattern  around  the  gap  was  Inspected  with  a 
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telescope  set  so  that  the  edge  of  the  rotor  was  out  of  focus. 
Adjustment  was  facilitated  by  the  fact  that  the  number  of  lines 
In  the  diffraction  pattern  changed  rapidly  with  small  motions 
of  the  test  piece.  When  the  diffraction  pattern  observed  was 
uniform  over  the  rotor  and  stator  (five  or  six  lines  were 
observable)  the  test  section  was  aligned. 

The  use  of  the  100  watt  light  bulb  for  the  Initial 
positioning  facilitated  alignment  especially  so  far  as  yaw 
of  the  test  piece  was  concerned. 

A  further  check  on  the  alignment  was  made  after  each 
run  from  the  "no-flow"  lnterf erograms .  In  these,  the  "picket - 

_  i 

fence  effect  mentioned  by  Howes  and  Buchele,  and  typifying 
alignment  of  the  test  section  were  clearly  in  evidence. 

5.4.2  Obtaining  "Plow"  and  "No-Plow"  Photographs 

After  alignment  of  the  test  piece  the  conventional 
Mach-Zehnder  interference  pattern  was  adjusted  for  fringe 
spacing  and  orientation;  further,  with  the  filter  removed,  the 
zeroth  order  interference  fringe  was  centerally  located  in  the 
gap  by  moving  the  translation  mount. 


XW.  L.  Howes  and  D.  R.  Buchele,  NACA  TN  3507,  p.  45.  Since 
the  light  source  is  always  of  finite  size,  not  a  geometric 
point,  non-parallel  rays  do  enter  the  test  section  and  cause 
interference  as  described  above.  The  "picket -fence"  effect 
is  caused  by  the  combination  of  two  opposed  patterns  of 
diagonal  fringes,  and  is  shown  in  Howes  and  Buchele's  Pig.  6. 

A  "No-Plow"  plate  from  the  present  study  is  shown  in  Pig.  5-2. 
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Figure  5-2 

No  Flow  Tnterferogram  showing  "Picket  Fence"  effect 
and  Dlffraotlon  Pattern 


< 
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The 

1. 

1/3  as 


2. 

3. 


interferometer  was  now  ready  for  use. 
procedure  followed  was: 

The  camera  was  focused  at  an  axial  distance 
recommended  by  Wachtell.^'2 

A  "no-flow"  photograph  was  taken. 

The  equipment  was  started  up,  and  equilibrium 


achieved. 


4.  A  "flow"  photograph  was  taken. 

For  the  runs  with  rotor  temperatures  up  to  and 
including 700R  *  the  "flow"  and  "no-flow"  photographs  were 
taken  on  separate  plates  with  the  filter  in  place.  For  the 
runs  at  rotor  temperatures  of900R  both  were  taken  on  the 
same  plate,  without  the  filter.  The  filter  was  removed 
because  a  single  high  intensity  flash  from  the  B-H6  lamp  when 
it  wa3  operated  by  capacitor  discharge,  as  it  was  for  these 
runs,  was  not  stiff icient  to  expose  a  plate  using  green  light. 


Since  only  a  few  interference  bands  resulted  because 
of  the  removal  of  the  filter,  and  since  the  bands  were  translated 
as  well  as  distorted  by  the  density,  the  "no-flow"  and  "flow" 
fringes  did  not  overlap  and  evaluation  was  possible. 


*G.  P.  Viachtell,  "Refraction  Error  in  Interferometry  of 
Boundary  Layer  in  Supersonic  Flow,"  Ph.D.  Dissertation, 
Princeton  University,  1951* 

2 

See  Section  4. 4. 3.1. 
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The  photographic  plates  so  exposed  were,  after 
development,  ready  for  evaluation. 

5.5  Use  of  The  Hot-Wire  Anemometer 

The  hot-wire  anemometer  described  In  Chapter  IV  was 
used  to  show  that  two  different  types  of  flow  can  exist  In 
the  annulus.  This  was  accomplished  both  with  measurements  of 
and ,  and  with  photographs  of  oscilloscope 
traces  of  the  signals  from  the  hot-wire.  During  these  runs, 
axial  traverses  were  made  with  a  single  probe  In  the  center  of 
the  gap  at  several  different  rpm  and  only  under  adiabatic 
conditions. 
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Chapter  VI 
Processing  the  Data 

In  this  chapter  are  presented  the  methods  used  to 
process  the  experimental  data. 

6.1  Overall  Heat  Transfer 

The  overall  heat  transfer  coefficient  ^  Is  defined 
as 

<PC 

where  ^  *»  Btu/hr-Ft2-deg  P 

2 

A  =  Area,  ft 
T  »  Temperature  deg  P 
Qc  o  Btu/hr.  transferred  from  rotor  to 

stator  by  convection 

6.1.1  The  Determination  of  Qq 

The  rate  of  heat  addition  to  the  cooling  water, 

Q  Btu/hr,  includes  radiation  transfer  between  the  rotor  and 
stator  as  well  as  convective  heat  transfer,  as  follows 

Q  »  Qc  +  ^radiation 

Q  was  determined  from  cooling  water  flow  rate  and  temperature 
increase. 


The  cooling  water  flow  rate  was  measured  by  weighing 
the  water  passing  through  the  cooling  jacket  over  a  time 
interval. 
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The  cooling  water  temperature  difference  was 
obtained  from  lron-conatantan  thermocouples  placed  aa  de¬ 
scribed  In  Section  3-6,  so  that  the  difference  In  millivolts 
was  Indicated  directly. 

As  the  cooling  water  temperature  averaged  80  deg  7, 
a  value  0.0288  mllllvolts/degree  was  used  to  determine  the 
temperature  difference.  The  readings  of  cooling  water 
were  taken  at  two  minute  Intervals  for  approximately  twenty 
minutes.  The  At  used  for  computing  Q  was  the  time  average 
of  these  10  odd  readings. 


Thus 

a.Jpis/Qfal ^  ^  (6-*) 


^radiation  was  determined  from  Stef an -Boltzmann's  Law  which 
for  an  annulus  is 


/A  t  ,  fats  f-A —  -  fl 


a 

Using  a  value  of  emlsslvity 

*  ^S+Jn.  *■ 


^radiation  for  the  preBent  *Pparat«8  becomes^ 

OmdHm  -$£JJ 


1W.  H.  McAdams, 
Company,  Ino., 


Heat  Trasmission,  (2nd  ed: 
New  York,  1942),  Table  Xlll 


MoOraw-Hlll  Book 
,  p.  59* 


2Ibid,  eq.  (14),  p.  3S4 
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Therefore  Qconvection  can  be  determined. 

6.1.2  The  Determination  of  A^*.  _ 

rotor 


Prom  the  cylinder  dimensions 
L  m  4.5  inches 

drotor  m  3*320  inches  to  3*332 
range  caused  by  thermal  expansion  at  the  higher 
temperatures 


/**  f<+ 

^326  u  76 


The  value  used  was 

Arotor  “  °»326  ft2 

6.1.3  The  Determination  of  (Trotor-Tstator) 

Trotor  was  taken  as  the  average  of  the  readings 
from  the  three  rotor  thermocouple  taken  immediately  before, 
and  after  measuring  the  water  flow.  At  first,  readings  were 
taken  every  3  minutes  but  the  variation  in  temperature  over 
these  short  periods  was  found  to  be  negligible.  When  changes 
did  occur  to  Trotor  they  were  generally  accompanied  by  changes 

ln  Tstator  80  that  (Trotor_T stator)  was  essentially  constant. 
The  change  in  (Trotor-Tstator)  during  a  run  was  about  3  degrees 
and  this  at  the  higher  temperatures. 
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During  these  runs,  three  thermocouples  were  used 
two  located  180  degrees  apart  at  X/L  *  1/2,  the  other  at 
X/l  ■  1/4,  and  90  degrees  from  either  of  the  first  two. 
Readings  between  one  of  the  thermocouples  at  X/l  »  1/2  and 
that  at  X/L  ■  1/4  agreed  to  within  1-2  degrees  Indicating 
that  axially  the  rotor  temperature  was  uniform*  The  third 
couple  however,  always  read  high,  and  was  flux  dependent. 
Indicating  perhaps  a  faulty  mechanical  Installation.  Since 
no  reason  was  found  for  rejecting  the  reading  from  this  couple 
Its  Indications  are  Included  In  the  average.  The  deviation 
from  the  reading  of  the  other  two  couples  was  never  more  than 
10  degrees,  this  at  the  highest  flux  rate,  which  when  averaged 
with  the  other  two  could  result  In  an  error  of  at  most  1 %  In 
overall  heat  transfer. 

6.2  Velocity  Distributions 

The  velocity  probe  Indicates  total  pressure  that  Is 
local  static  pressure  plus  velocity  pressure,  modified  In  some 
way  still  undetermined  analytically  or  experimentally  by  the 
turbulent  fluctuations.-1-  Unfortunately,  direct  computation  of 
the  velocity  is  Impossible  even  If  the  turbulent  fluctuations 
are  Ignored,  since  the  local  static  pressure  Is  unknown. 

A  graphical  Integration  procedure  is  required  to  determine  the 
velocity  In  the  gap  and  Is  developed  below.  This  method  has 


*See  Section  4. 3. 1.3*3 
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been  used  by  0.  I.  Taylor,^-  and  P.  L.  Wattendorf.^ 


Prom  the  r  -  Component  of  the  equations  of  .motion 
one  can  show  that 


where  Interactions  between  molecular  and  convective 
fluctuation  terras  have  been  neglected. 

The  last  term 

fa  f  Y*  everywhere  In  the 

flow  and  can  be  Ignored. 

At  the  probe,  it  is  difficult  to  state  what  effects 
the  skin  friction  inside  the  probe  mouth  and  the  lags  of  the 
manometric  system  will  have  regarding  the  effects  of  the 
turbulent  fluctuations  on  the  indication  at  the  manometer. 


AI30,  as  discussed  in  Section  4. 3. 1.3. 3  various 
investigators  have  treated  the  turbulent  fluctuations  at  the 
probe  in  several  ways,  from  ignoring  them  completely  to 
including  in  the  total  pressure  the  total  head  of  the  total 
velocity  vector. 


*G.  I.  Taylor,  "Distribution  of  Velocity  and  Temperature  Between 
Rotating,  Concentric  Cylinders,"  Proc.  Royal  Society  of  London, 
ser.  A,  vol.  157*  no.  892,  Dec.  2,  1936,  p.  498. 

2 

p.  L.  Wattendorf,  "A  Study  of  the  Effect  of  Curvature  on  Fully 
Developed  Turbulent  Plow,  Proc.  Royal  Society  of  London,  ser.  A, 
vol.  148,  no.  865,  Peb.  1935,  p.  570. 
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as 


If  the  total  pressure  at  the  probe  mouth  Is  taken 


(6-4) 


then,  equation  (6-4)  together  with  equation  (6-3)  yields 

-0  rfV"] 

No  statement  can  be  made  regarding  the  magnitude 
of  the  last  term.  It  has  been  neglected  by  other  Investigators. 


Therefore , 

0  (6-5 ) 


The  value  of  jS 
stator* 


Therefore , 


at  the  stator  Is  known  and  is  a  constant. 


Starting  at  the  stator,  this  equation  can  be  numerically 
Integrated  to  give  the  distribution  of  static  pressure  In  the 
gap.  The  trapezoidal  rule  was  used. 


Knowing  ■&  and  ,  using  equation  (6-5)  one  can 


find 


127 


6.3  The  Determination  of  Temperature  Distributions 

6.3*1  Density  Differences 

As  shown  in  Appendix  X  ,  from  point  to  point  on 
an  interferogram 

A 


_  6644444 


(6-8)a 


If  the  "no-flow"  and  "flow"  interferograms  in 
Fig*  6-1  a)  and  b)  are  superimposed.  Fig.  6-1  c)  results. 

i.u 

If  there  were  no  density  change  the  N6  fringe  would  go  through 
point  B.  Instead,  the  1st  fringe  passes  through  B.  The  fringe 
shift,  £  is  given  by 

- 

y\  -  5460  A 
/  “  ft 

1 

^  -  0.3635  X  10"2  fjp 


That  is. 


(6— 8)b 


and 

thus 


Afm  13.1414  x  * 

6. 3. 1.1  The  Measurement  of  Fringe  Shift 


y t'A 


(6-8)c 


The  Interferograms  were  evaluated  in  a  Jones  and 
Lamson  Pedestal  Comparator,  of  optical  magnification  20x. 


Q.  M.  Edelman,  and  M.  H.  Bright,  "The  Specific  Refractivity 
of  Oases  for  Various  Wavelengths  of  Light,"  Gas  Turbpie 
Laboratory  Report  No.  6,  Cambridge  Massachusetts,  M.I.T., 
•May*  1948,  p.  5. 
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The  lnterferograms  had  a  Magnification  factor  of  l/2x,  so  that 
the  overall  Magnification  fron  the  gap  to  the  soreen  of  the 
coaparator  was  lQx.  The  comparator  was  capable  of  measurement 
of  linear  distances  to  0.0001  Inch. 


Two  Ink  lines  were  drawn  on  the  outside  of  the 
quartz  window  nearest  the  camera,^  In  order  to  align  the 
plates  and  also  to  serve  as  a  reference  mark  to  locate  the 
surface  of  the  rotor  (see  Sec.  6. 3 .1.2. )  Therefore  the 
angle  &  In  Pig.  (6-1  )c  could  be  measured,  as  could  ^Fringe 
from  the  "no-flow"  plate. 


To  determine  A  ,  let  the  measured  coordinates  of 
point  B  be  (x,  y)  and  those  of  point  A  be  (x^,  y^) 

From  Fig.  (6-1 )c,  A  is  seen  to  be 

A=x-xa-  ^->6 

Uccfi 


or 


This  technique  was  used  by  DeFrate,  Sc.D.  Dissertation, 

MIT,  195LO,  and  0.  P.  Wachtell,  Fh.D.  Dissertation,  Princeton 
University,  1951.  Although  the  camera  was  not  focussed  on 
the  reference  mark,  the  Image  of  the  reference  mark  on  the 
photographic  plate  was  symmetrical  since  the  Image  of  an 
out-of -focus  symmetrical  object  is  symmetrical,  depending, 
of  course,  on  the  quality  of  the  camera  lens.  The  camera 
lens  used,  an  Aero-Ektar  f/2.5,  7  Inch  lens  Is  of  excellent 
quality. 


““  Alo~F/ow 
Fringe 

rotor  Distribution 


r tor 

_  F/g  £>-/& 

~~  Distortion  of 
Fringes  during 

F/Ow 


F/g  6-f  c 
Superposition 
of 

Flow  and  No-Flow 

Inter  Aerograms 
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6.3.1 .2  The  Detemlnatlon  of/^  Reference  Density 

The  Ideal  method  of  determination  of  Is  one  In 
which  the  Interferometer  alone  la  used.  This  can  be  accomplished 
If  unflltered  mercury  light  Is  used  as  a  means  for  assigning 
numbers  to  the  filtered  light  fringes  through  means  of  no¬ 
flow  and  flow  interferograms  taken  with  the  unflltered  light. 

In  this  manner,  rotor  temperature  could  be  determined  and 
checked  against  that  measured  from  the  rotor  thermocouples 

This  was  not  done  for  two  reasons:  First,  there 
existed  a  temperature  gradient  in  the  window  (See  Sections 
4. 4. 3. 2,  6.3. 1.2.1,  and  AppendixjT)  which  would  prevent  the 
numbering  of  the  fringes,  and  secondly  the  location  of  the 
rotor  surface  on  the  "flow"  lnterferogram  Is  only  apparent } 
that  Is,  the  rotor  Is  not  where  It  seems  to  be,  because  of 
refraction  and  diffraction  In  that  region. 

Because  of  the  effects  mentioned,  the  closest  meas¬ 
urement  to  the  rotor  surface  was  of  the  order  of  0.001  -  0.005 
Inches  in  the  real  plane,  the  low  figure  at  low  density  gradients 
the  high  one  at  high  density  gradients. 

Since  near  the  wall  is  constant.  It  was  calculated 
from  overall  heat  transfer  measurements  and  the  rotor  temperature 
as  determined  from  the  rotor  thermocouples. 

That  is 

For  ^7*,(^r~‘Srf»/)can  be  determined. 
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Por  the  present  experiments,  AT  between  0  and 
0.005  Inches  corresponds  to  at  most  ^  between  0  and  6.5. 

In  this  region,  and  particularly  for  ^  between  0  and  4, 

^  ^  to  the  best  yet  experimentally  determined.  That 

is,  this  region  constitutes  the  so  called  "laminar  sublayer". 


In  essence  then  the  effects  of  turbulence  are  being 
Ignored.  This  has  been  checked  a  posteriori,  and  there  results 
little  or  no  error.1 

Thus,  ft  is  determined  and  is  obtained  from 
equation  (6-8)c. 

6. 3. 1.2.1  The  Determination  of  the  Glass  Correction  to 

The  value  of  the  temperature  gradient  existing  in 
the  quartz  windows  has  been  determined  experimentally  for  the 
different  operating  conditions.  The  results  are  presented  in 
Appendix  TL,  Equation  (4-3)  is  used  to  determine  the  correction 
made  to  the  measured  densities. 

The  change  in  optical  path  lengths  In  the  quartz 
window  between  points  1  and  2  caused  by  a  temperature  gradient 
in  the  quartz  Is 

,AZ  OK  +£/***)  (4.3) 

Usually,  there  exists  on  the  flow  lnterferogram  some  area,  and 
therefore,  point  whose  density  is  known  accurately,  such  as  the 
free-stream  density  for  the  boundary  layer  on  a  flat  plate. 
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Prom  Seotlon  4. 4. 3. 2 


a.  -  1.46013  at  18  C  (-  64.4  F) 

-  1.46019  at  75  P 

C<.  -  1.7  x  10”6  /deg  C 

*  3 

-  .32  x  10"6  /deg  P 

^T-  11.29  x  10*  /deg  C  -  6.272  x  10"6  /deg  P 

Since  ^  •  0.125  Inches 

-  1.203  x  10"6  (Ttfi  -  T,jM  ) 

•  €  correction  S 

^  -  5460  x  10*8  inches  -  21.496  x  10"^  Inches 
2.54 

Therefore,  /degree  temperature  difference  •  1*203  x  10  ? 

21.496  x  10 


-  0; 055964 


If  the  gradient  In  the  glass  Is  expressed  as  6^ 
degrees/inch,  then 

^/inch  -  ^(0.055964) 
since  ^*.(13.1414  x  10“*)#/ft3 

-  dfe  (0.055964)  13.1414  x  10”4 

^  /0.001  Inch  -  ^  (0.055964)  13.1414  x  10“7 

-  <^(.73541,  10-7 )#/n3/0-001ln(ih 


The  correction  Is  small  for  small  differences  In 
gap  radius.  Por  the  case  7000  rpm,  Trotor  ■  890  R,  at  a  dls 
tance  of  0.030  Inches  from  the  wall,  the  correction  amounts 


to  2  deg  P. 
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6.3*2  The  Determination  of  Radial  Poetlon  on  the  Inter¬ 
ferometric  Plate 


Consider  the  drawing  of  an  Interferometric 
plate  shown  In  Pig.  6-2. 

The  reference  lines  can  not  be  drawn  absolutely 
parallel  and  perpendicular  to  the  radius  going  through  the 
reference  point  (A,  F);  hence  the  angularity  effect  when  the 
long  reference  mark  Is  placed  on  a  horizontal  center  line  in 
the  comparator. 

The  point  (A,  F)  Is  determined  from  average  meas¬ 
urements  of  the  points  (1,2),  (3,*0,  (5*6),  and  (7,8)  and 
for  all  the  runs  measured  the  location  of  (A,  F)  relative  to 
those  points  never  varied  by  more  than  0.0002  Inches  in  about 
0.012  inches. 


The  points  (B,  E),  and  (C,  D)  are  likewise  measured. 
On  both  the  flow  and  no-flow  interferograms  the  area  through 
the  gap  about  the  reference  point  was  blocked  off  In  the 
reference  path  of  the  Interferometer  so  that  the  interference 
fringes  would  not  obscure  the  surfaces  of  the  rotor  and  stator. 
The  angle  is  measured.  This  is  of  the  order  0°9'.  With 

these  quantities  the  following  can  be  determined: 


since  cos 


The  accuracy  of  measurement  of  pr  is  of  no  concern 
9'  -  0.99999*  and  can  be  measured  to  within  l1 


of  arc 


*9:  £'*- 

De^rmioa'/'/on  e>-/ 
/Q?d/a/ ^}s//~/or)  on  &n 
-Zn  -/e/~/er&jr<rm 
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The  radial  location  of  the  reference  point  /Vef ^ 
determined  from 

E  -  D  1 

M  -  Magnification  Factor  ■  - . 

cos  (rstator-rrotor'True 

-  E  -  D 
cos  I?  (.377) 

^ref  “  1*66°  M  +  (^ref "rrotor)piate 


Since 


j%  <M|  \ 

ref "rrotor 'plate 


P  -  D 
cos  ^ 


f ref 


1.660  M  +  L  D 
cos^ 


Since  the  reference  point  lies  In  an  area  of  very 
small  density  gradients  (and  hence  Is  not  subject  to  refraction 
effects),  and  since  the  camera  location  Is  unchanged  from  the 
no-flow  plate  to  the  flow-plate  one  can  write,  on  the  flow 
plate  that 

Zp  -  rpef  sln^  +  (A  -*) 
and  rref  cos  &  -  (F-^  ) 

where  X  and  y  are  the  location  of  the  point  of  Interest. 
Consequently, 

rplate  “  ^ ^ 


and 


r 


rplate 

M 


In  practice,  and  ££  are  taken  as  being  on  the  centerline 
of  a  fringe. 
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6.4  The  Determination  of  VjT  vs.  J/'7* 
Vjf  1b  defined  as 

9  A. 

*  Up*  sC*£n 


for  the  stator  and  as 


6  7^/a. 


Ylv+tx. 


for  the  rotor. 

j/ior  the  stator  and  rotor  are  defined  cuJ 
at  the  stator,  and 

rr79% 

at  the  rotor. 

To  calculate  these  quantities,  and  %  are 

required. 


Direct  measurements  of  shear  stress  were  not  made. 

If  ^  Is  evaluated  by  using  the  V*  measured  closest  to 
the  wall  and  forming  ^  an  error  Is  made.  First, 

the  closest  to  the  wall  that  V^.  was  measured  was  0.005  Inches 
and  0.008  Inches  for  the  stator  and  rotor  respectively.  The 
value  of  d(Ve/r)/  dr  at  these  points  Is  different  from 
the  value  at  the  wall.  This  Is  easily  shown  by  using  the 
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formulation  by  Deissler1  for  near  the  wall  which  la 

and  can  be  put  In  the  form  t  x  . , 

where  la  an  experimental  conatant  with  the 

value  1/8. 


Thua 


where 


a  *■  4- 


Using  the  data  presented  by  Sleicher2  for  ^  vs.  V  ,  and 
estimating  that  for  the  experimental  run  at  5000  rpm, 

no  heat  transfer,  one  finds  that  at  a  value  of  (Jr  -8. 

»  8,  U?  »  7  and  has  Increased  to  the  value  0.5* 

Let  £  represent  the  shear  stress  at  a  small 
distance  from  the  rotor.  Then 


Further 


and 


s  Jfo&r  frtrtrr 


1C.  C.  Lin,  Turbulent  Flows  and  Heat  Transfer  (Princeton  Univer¬ 
sity  Press,  Princeton,  New  Jersey,  1959)*  Section  E  by  R.  G. 
Deissler,  p.  230. 

2C.  A.  Slelcher,  Jr.,  Trans.  ASME,  vol.  80,  1958,  Table  1, 
p.  699. 
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Since  Jr+tn  ~  r 


This  Is  only  an  average  case.  The  situation  Is 
more  serious  at  7000  rpa. 

Secondly,  a  plot  of  vs.  ^  although  smooth 

did  not  close  to  >0,  but  Indicated  that  a  probe  displace¬ 
ment  effect  was  present.  See  Pig.  6-3. 

Accordingly,  the  following  was  done: 

Slelcher  measured  his  temperature  distributions  more 
accurately  than  his  velocity  distributions.  Unwilling  to  give 
up  the  quality  of  his  temperature  measurements  he  synthesized 
a  velocity  distribution  from  the  data  of  Relchardt,*  Delssler,^ 
Laufer^  (but  only  used  those  points  for  Re  ■  500*000)  and  his 

■^H.  Reichardt,  "Heat  Transfer  Through  Turbulent  Friction  Layers," 
NACA  TM  1047,  1953. 

2R.  G.  Delssler,  "Analytical  and  Experimental  Investigation  of 
Adiabatic  Turbulent  Flow  In  Smooth  Tubes,"  NACA  TN  2138,  1950. 

O 

J.  Laufer,  "The  Structure  of  Turbulence  In  Fully  Developed 
Pipe  Flow/'  NACA  TN  2954,  1953. 
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own  data.  The  synthesized  values  are  presented  In  the  table 
below,  and  Is  taken  from  Slelchers  previously  mentioned  ASME 
paper. 

Table  6-1 

Experimental  Distribution  of  Velocity  Near 
a  Pipe  Wall. 

2  3  4  5  6  8  10  12  14  16  18  20 

t£h  2  3  3.9  ^.8  5.6  6.9  8.1  9  9.9  10.6  11.2  11.7 

These  data  are  mostly  from  adiabatic  runs  but 

2 

analysis  shows  that  under  greater  rates  of  energy  transfer 
and  higher  temperatures  than  encountered  In  the  present 
experiment  ^^vs,  y^near  the  wall  does  not  vary  to  any  extent 
greater  than  the  experimental  error  with  which  one  can  measure 
the  velocity  there.  Deviations  do  not  occur  until  11. 

A  plot  of  these  values  Is  given  In  Pig.  6-4.  For 
the  present  experiment  these  have  been  used  to  determine  the 
value  of  the  shear  stress  at  the  wall  and  the  probe  displace¬ 
ment  by  the  following  scheme: 


The  data  available  are  the  measured  valuesof 
the  experimentally  determined  velocities. 


1C.  A.  Slelcher,  Jr.,  Ph.D.  Dissertation,  University  of 
Michigan,  1955. 

2R.  Q.  Deissler  and  C.  S.  Elan,  NACA  TN  2629,  Pg.  33,  Pig.  4 
entitled  Predicted  generalized  velocity  distribution  for  flow 
of  gases  with  heat  addition  at  Prandtl  number  0.73. 


lia 
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1.  Quess  a  value  of 

2.  Flna  ^  '  Jjt- 

3.  Proa  the  plot  In  Pig.  6-4  find  the  value  of 

4.  U+  “*</, 

*7/  Jft, 

5.  y,~" tyr>  /  the  dlaplacement 

6.  Plnd  y/l  *■ 

7.  Prom  the  plot  find  l/  + 

8.  Plnd 

9-  54-^2  -4 

10.  If  tDf  repeat  steps  1-9. 

11.  When  ^=Z},using  the  value  of  determined  In 
Step  1,  find  the  value  of  at  the  other  wall  from 

12.  Plnd  2),  for  the  first  data  point. 

13.  Find  .  This  must  be  the  same  as  If  the 

functional  relationship  given  In  Pig.  6-4  Is  to  hold  for  the 
present  case. 

For  the  present  data  one  finds  that  >  and 

Indeed,  44  *4  on  both  walls  except  at  the  lower  temperature 
high  speed  runs 


This  equality  shows  that  the  data  follow  the 
functional  relation  given  by  the  values  In  Table  1.  They 
have  been  In  no  way  forced  since  there  are  more  known 
quantities  than  there  are  equations. 

The  displacements  D2  never  varied  from  by  more 
than  0.0005  Inches,  except  for  2  cases  at  5000  rpm. 

The  displacements  found  over  the  rotor  and  stator 
were  In  the  same  general  direction  as  Indicated  by  the  refer¬ 
ences  cited  In  Section  4. 3. 1.3. 2,  that  Is  the  aerodynamic 
center  of  the  probe  Is  shifted  toward  the  regions  of  higher 
velocity.  Thus,  velocities  measured  over  the  stator  are 
further  from  the  wall  than  and  those  over  the  rotor  are 

closer  than  tyro  • 

1/  ^  V* 

Thus  is  known,  Uf  is  known  and  V&  vs.f/ 

can  be  calculated. 

Consider  the  application  of  the  method  to  the  data 
from  a  run  at  7000  rpm,  with  no  heat  transfer. 

For  the  stator 

5  24.53 

9  32 .66 

Try  Uf  -  2.9.  The  value  of  is  calculated, 

is  obtained  from  Fig.  6-4,  and  ^7-  Is  calculated. 
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Displacement  in 

Try  u#  =  2.9  thousandths  of 


yT,m 

u+ 

y+ 

yT 

inches 

5 

8.46 

10.7 

7.53 

2.53 

9 

11.27 

18.25 

12.84 

3.84 

Try  u#  =  3.3 

5 

7.43 

8.7 

5.33 

+  0.38 

9 

9.9 

14.1 

8.72 

-0.28 

Try  u#*  3  •  1 

5 

7.91 

9.6 

6.32 

1.32 

9 

10.54 

15.75 

10.36 

1.36 

This  is  satisfactory. 

Prom  the  principle  of  Conservation  of  Angular  Momentum 

..  .  rstator 

*, rotor  rrotor  *, stator 


The  friction  factor  is  defined  as 


145 


where  yg^  Is  determined  from 


I be 


\!or  " 


14  rotor  /'rotor 

2  rc 

i4 rotor  /* rotor 


■  —L. 

J  (  /stator  f /'rotor) 


and  is  approximately  the  velocity  In  the  center  of  the  gap. 


Therefore , 

* 

and 


ii<L 


6.6  The  Determination  of 

Prom  equation  (2-20) 


(2-20) 


All  the  quantities  in  this  equation  have  been  previously 

determined  except  - - 

dr. 

To  find  l/*/f  was  plotted  against  r.  Over  the 

rotor,  data  points  were  at  »  8,  10,  15,  20,  25,  35*  ^5, 

/ 

60,  85,  and  185.  Over  the  stator  **  5*  9*  12,  15*  20, 

30,  40,  65,  90  and  140.  These  are  modified  by  the  probe 
displacement.  The  points  are  not  evenly  spaced.  A  curve  was 
drawn  through  the  experimental  points  and  every  0.0025  inches 
of  r  a  value  of  was  taken  from  the  curve. 
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The  derivatives  were  obtained  as  outlined  In 
Lanczos.^  This  method  places  a  least-square  parable  through 
five  evenly  spaced  points.  The  derivative  is  then  taken  at 
the  central  point. 


The  least-squares  parabolas  used  In  the  differ¬ 
entiation  procedure  go  through  the  experimental  points  with 
negligible  differences  except  In  a  few  isolated  instances. 


One  such  Instance  is  presented  below 
For  example  at  5000,  A  £  * 


Rotor 


(yj>„ 

<Vr>m 

10 

344.06 

15 

297.79 

20 

273.36 

25 

264.58 

35 

250.92 

45 

243.54 

Least  squares 
parabola 

343.04 

298.25 

273.74 

263.1 

250.89 

243.63 


60  R 
Stator 

H*  (Vrl 

5  77.27 

15  135.13 

20  154.74 

30  168.55 

40  178.03 


(W 

Least  squares 
parabola 

73.9 

136.1 

154.07 

168.56 

176.4 


Although  evaluation  otG»/js  requires  a  knowledge  ot  T  , 

wa  11 

the  value  of  (1  +€"/&)/£  can  be  found  from  quantities  measured 
experimentally.  See  Section  8.5,  where  (1  +£*4/3/)//'  is  given 


as 


0  ~  - 


fjM#  )  / 

2  2/#  y*  '  rld(Vt/r) /Jr 


(8-3) 


1C.  Lanczos,  Applied  Analysis  (Prentice-Hall,  Inc..  Ensclewood 
Cliffs,  New  Jersey,  1956),  PP.  315-324. 
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6.7  The  Determination  of  Dimensionless  Quantities 

6.7.1  Taylor  Number 

The  Taylor  Number  Is  defined  as 

where  JSG  Is  the  kinematic  viscosity  at  temperature 
^  rotor  is  determined  from 

6-7-2 

The  quantity  Is  formed  from  the  product  of 

the  friction  factor  and  the  Taylor  number  for  flow  at  Taylor 
numbers  less  than  the  critical.  It  Is  given  by 

/ 

**  ti&der /fearfor- y&>r) 

where  the  friction  factor  is  defined  on  the  basis  of  rotor 
velocity. 

6.7.3. 1  Nusselt  Number,  Nu 

The  Nusselt  number  Is  given  as 

A f(4  *= 

.here  A  1.  evaluated  at  trotor  \  T,tl>t°r 


148 


6. 7. 3. 2  (Nu5so.lt  Number)  .  .  =  Nu, 

conducting  cond 


The  value  obtained  for  the  Nusselt,  number  In  laminar 
flow  for  the  present  apparatus  io  Nucc.]Ui  and  is  given  as 
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Chapter  VII 
Estimates  of  Errors 

In  this  chapter  an  attempt  Is  made  to  estimate  the 
errors  In  the  experiment. 

A  correction  for  a  determinate  systematic  error,1 
that  caused  by  the  temperature  gradient  In  the  quartz  windows, 
has  already  been  made.  See  Section  4. 4. 3. 2. 

There  Is  at  least  one  systematic  error  which  can 
not  be  estimated,  that  caused  by  neglecting  some  of  the  turbu¬ 
lent  fluctuating  velocities  In  calculation  of  the  velocity 
from  probe  measurements  (See  Sections  4. 3. 1.3*3  and  6.2). 

Another  error.  In  which  the  personal  Judgement  of  the  Investi¬ 
gator  Is  Involved,  cannot  be  evaluated. 

This  Is  the  case  In  calculation  of  for  50 

100.  In  the  central  portion  of  the  gap,  velocity  measure¬ 
ments  were  taken  at  radial  intervals  of  from  0.020  inches  to 
0.100  inches,  as  contrasted  with  Intervals  of  0.005  Inches  near 
the  wall.  Since  the  determination  of  required  differenti¬ 

ation  of  the  velocity  data,  as  outlined  in  Section  6.6,  a  smooth¬ 
ed  curve  drawn  through  the  data  near  the  central  portion  of  the 
gap  necessarily  involved  personal  Judgement. 

The  effect  of  the  scatter  In  the  data  in  this  region 
on  subsequent  calculations  Is  small. 

1Y.  Beers,  Theory  of  Errors  (2nd  ed.,  Addison-Wesley  Publishing 
Company,  Reading  Mass.,  1957 )»  P*  4 • 


It  Is  mentioned  that  In  the  evaluation  of  the 
conduction  errors  perfect  contact  Is  assumed  In  the  Instal¬ 
lation  between  the  thermocouple ,  the  cement  and  the  sides  of 
the  hole.  As  already  stated  (see  Sec.  6.1.3)*  the  Indication 
of  one  thermocouple  In  the  rotor  differed  from  the  other  two 
by  an  amount  which  was  flux -dependent.  Indicating,  perhaps  an 
air  pocket  In  the  hole  In  which  the  thermocouple  Is  Immersed. 

The  estimated  conduction  errors  computed  do  not 
Include  these  effects. 

Finally,  the  errors  calculated  are  to  be  regarded 
as  reasonable,  order  of  magnitude  estimates  of  the  errors  In 
the  determined  quantities. 

In  evaluating  the  errors  It  Is  assumed  that  the 
error  of  an  experimental  result  can  be  broken  Into  components 
determined  by  the  functional  relation  between  the  final  result 
and  those  experimentally  determined  quantities  which  enter 
into  the  calculation  of  the  final  result. 

Wilson*  shows  the  following: 

If  the  final  result  Is  related  to  the  components 

J as 

(/•  ■  ■■■*») 
then  the  square  of  the  error  In  ^  Is 

*E.  B.  Wilson,  Jr.,  An  Introduction  to  Scientific  Research 
(McGraw-Hill  Book  Company,  Inc.,  New  York,  1952 )  p.  272. 
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7.1  Overall  Heat  Transfer 

Equation  (6-1)  defines  the  overall  heat  transfer 
coefficient  Al  as 

Qc  /H ro£or  (^rbier  —  T^tatcr)  *  / Amtor  AT 

and 

(&h/h)  86  (%&c/0^)  +  C&Afbt&r/Arviorf  +  { Sw/A  T) 

7.1.1  The  Error  in  the  Determination  of  Qc 

The  quantity  tyOc.  in  turn  can  be  broken  into 

components 

-  f ^S/^r  'Cooling  ^ter  +  i(^7^  cooling  water 
(lbs/hr)cooiin”  water  aT  coollng  water 

Considering  first  the  error  in  the  measurement  of  the  cooling 
mater  flow  rate, 

^  ( lbs/hr )eoonng  Water  _  b  (lbs)  +  b  (time  Interval) 
( lbs/hr )COQllng  water  “  (lbs)  time  interval 

The  scale  used  mas  capable  of  measuring  to  mithln  at  least  2  oz; 
that  is  a  2  oz  melght  placed  on  the  scale  mhen  balanced  caused 
a  correct  deflection  in  the  lever  are.  The  melght  of  mater 
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measured  was  of  the  order  of  20  lbs.  -  320  oz. 

A  1/5  second  stop  watch  was  used  to  determine  the 
time  Interval.  Comparison  of  the  stop  watch  with  an  electric 
clock  over  a  thirty  minute  time  Interval  revealed  no  dis¬ 
cernible  error.  On  the  assumption  of  an  error  in  starting 
and  stopping  of  0.4  seconds  each,  one  second  Is  taken  as 
the  error  in  stop  watch  reading.  The  time  Interval  over 
which  the  weight  of  water  was  measured  was  15  minutes. 

Thus 

&  (IbsAr)  _  2  1 

lbs/hr  320  900 

The  error  In  the  measurement  of  the  cooling  water 
r.yi&perature  difference  can  be  divided  Into  three  parts.  They 

1.  Conduction  errors 

2.  Calibration  errors 
and  3.  Indication  errors 

Recall  that  this  measurement  was  made  by  Joining  the 
thermocouples  so  that  the  difference  In  emf  was  Indicated. 

The  conduction-error  was  found  to  be  negligible 
(Section  4. 2. 5.1). 

The  calibration  error  Is  of  no  consequence  also 
since  the  two  thermocouples  were  made  from  conseoutlve  sections 
of  the  same  strand  of  wire. 
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The  error  In  Indication  was  considered  to  be  a 
maximum  of  1/4  deg  (See  Section  4.2.4)  the  accuracy  with  which 
the  Portable  Precision  Potentiometer  can  be  read.  The  potenti¬ 
ometer  was  not  calibrated  during  the  present  experiment.  Since 
It  was  not  mishandled  It  is  assumed  that  the  manufacturers 
specified  accuracy  Is  valid. * 

The  error  In  cooling  water  determination  was  then 

1/4  _  _1 

10  40 

The  total  fractional  error  In  ^^is 


+  i- 

l  C±)  320  900 

40 

The  Determination  of 

^*Arotor 

Arotor 

Arotor  “  ^  ^rrotor^  ** 

Measurement  of  rrotor  was  made  to  within  0.001  Inches 
In  1.660  Inches;  L  was  measured  to  within  1/64  Inches  In  4.5 
Inches . 

-  Ar°t.°.r  1.  then 

A 

rotor 

£Arotor  m  1  +  1  +  1  +  1 

Arotor  "  1660  64  x  4.5  1660  288 

1 Verbal  communication  with  Mr.  Bowe,  Leeds  and  Northrup  Company 
(New  York);  stated  accuracy  +0.05J*  of  reading  +0.006  mv. 
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7.1.3  The  Error  In  Determination  of  AT 

The  errors  In  the  determination  of  AT  can  be 
classified  as 

1.  Error  In  determination  of  stator  temperature 

a.  Conduction  error 

b .  Placement  error 

c.  Calibration  error 

d.  Indication  error 

and 

2.  Error  In  determination  of  rotor  temperature, 
composed  of  the  error  In  the  thermocouple  and  the  error  In 
slip  rings. 

7. 1.3.1  Error  In  determination  of  Stator  Temperature 

a .  Conduction  error 

This  was  found  to  be  0.13  F  (Section  4. 2. 5. 1.1) 

b .  Placement  error 

This  error  was  found  to  be  maximum  of  0.4  P 
a  minimum  of  0.2  P. 

The  correction  for  this  has  been  made  in  the 

calculation  of 

c.  Calibration  Error 

Prom  the  monograph  by  P.  H.  Dike,  entitled 
"Thermocouple  Thermometry"1  the  nominal  limits  of  error  for 

1P.  H.  Dike,  "Thermocouple  Thermometry,"  (Leeds  and  Northrup 
Company,  Philadelphia,  195^) >  P*  18. 
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stock  selected  chromel-alumel  thermocouple  wire  Is  5 
degrees  In  the  range  32-600 Accordingly,  since  the  stator 
tenperature  was  about  100  degrees  the  error  Is  5  P. 

d.  Indicating  error 

As  with  the  cooling  water  thermocouples,  this 
error  is  0.25  P. 

The  total  error  Is  0.13  +  0.25  +  5  ■  5. 38  P. 

^tator)  .  .0538 
/stator 

7. 1.3. 2  Error  In  Determination  of  Rotor  Temperature 

a.  Conduction  Error 

The  error  was  found  to  be  negligible. 

(Section  4.2. 5. 1.2) 

b .  Placement  Error 

This  error  was  found  to  be  1.5  P  under  maximum 
conditions  of  heat  transfer  (Section  4. 2. 5. 2)  and  0.7  under 
minimum  conditions  and  the  correction  has  already  been  made 
In  the  calculations. 

c.  Calibration  Error 

This  error  Is  given  by  Dike  as  Vjt>  of  Indication 
above  500  P.  Since  the  millivolt  Indication  at  the  highest 
temperature  difference  was  of  the  order  of  10#this  error  Is 
of  the  order  of  0.1  millivolt  or  3  P  at  800  R  and  2  F  at 
600  R. 
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d.  Indicating  Error 

As  with  the  cooling  water  thermocouples  this 
error  is  1/4  P. 

e.  Slip  Ring  Error 

As  stated  In  Section  4. 2. 2.1  this  error  was 

found  to  be  3 


The  total  fractional  error  In 


maximum  values  Is 


or 


0.0538  +  0.0329 

0.0238  +  0.0131 


using  the 

0867  max 
0669  min 


The  error  In  A>  Is  then  determined  from 


- 

1  + 

320 

JL  + 
900 

hm 

,  _1 _ + 

1  _ 

0.004 

1660 

288 

*6  AT 
T 


0.0867  max 
0.0669  min 


and  Is  calculated  as 

-a  8.7$  *ax  to  6.7#  min 


Note  that  a  large  part  of  the  error  calculated 
results  from  the  error  In  the  stator  thermocouple.  Thermo¬ 
couple  readings  taken  before  startup  of  the  equipment  after 
overnight  shutdown  gave  Indications  for  all  the  thermocouples 
(cooling  water,  rotor  and  stator)  which  were  within  one  degree 
P  of  each  other. 


The  error  given  here  therefore  assumed  to  be  on  the 


high  side 
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7*2  Velocity  Distribution 

7*2.1  The  Error  In  Determination  of 


In  accordance  with  the  discussion  of  Section  4. 3.1. 3*3 
the  error  caused  by  neglecting  turbulent  velocity  fluctuations 
(See  Section  6.2.1)  cannot  be  estimated. 

The  accuracy  of  velocity  determination  depends  on 
the  accuracy  of  the  manometer  used. 


The  velocity  Is  related  to  the  reading  on  the 


manometer  by 


M  ./■ 
v,  * 


Wttl 

(G~p) 


The  quantity 


Is  constant  and  Is  estimated  to  be  of  the 


order  0.003  inches  of  fluid  as  a  maximum,  not  because  of  any 
limitations  of  the  micromanometer  but  because  of  the  method 
used  to  bracket  the  reading  (See  Sections  4.3.3,  and  5*3). 

The  minimum  value  of  the  quantity ^^"^^occurred  at 

yr9  Stator  -  5  at  3000  rpm,  £7~  -  330  %  and  was  of  the 

order  0.020  Inches,  so  an  error  of 

jjk  _  1  0.00  7  ^ 

^  '  2  0.020  “ 


1  0.00 

2  0.020 


-  7.5* 


Is  Indicated  at  that  point. 


At  Vr,  stator  ■  9,  for  the  same  conditions, 


■  0.037  for  an  error  of 
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1  x  0.003 

2  0.037 


4# 


At  7000  rpm,j7-.33o^stator 


and  the  error  Is  2.4#. 


5,  fa-p) 


0.064 


At  5^7"  -  9  for  the  same  run,^/£-^)  -  0.140  and  the  error 

Is  1.1#. 


7.2.2  The  Error  In  Determination  of  Radial  Position,  r 


The  error  in  radial  position  r  is  estimated  to 
be  of  the  order  of  0.0005  inches,  by  the  method  used  for  probe 
location  (.see  Section  4. 3. 1.2  and  Fig.  4-£). 


7.2.3 


The  Error  in  Determination  of 


The  quantity 


dr  r  ' 


can  be  written 


dr 


If  r-^,  r2,  and  r^  are  three  evenly  spaced  radii  then  a  good 
approximation  for  the  present  purposes  is 

V&  3  ~  (/&■) 

~7hT  ~  ** 


The  error  in 


is  then 


dr 


iVti  JH,i  S(rj-n) 

%-  til,  (n-o 

The  errors  in  V^.  are  generally  small  in  the 

area  of  greatest  rate  of  change  in  ^ (i (at  about  •  20) 

dr 

and  are  taken  as  1#. 
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Purther,  the  accuracy  with  which  the  probe  oan 
be  positioned  over  an  Interval  la  of  the  order  0.0003  Inches, 
which  Is  the  accuracy  with  which  a  micrometer  can  be  read, 
admittedly  with  sowe  effort  .However,  If  the  traverse  Is 
always  made  In  the  same  direction,  eliminating  the  problem  of 
backlash  In  the  micrometer  screw  and  barrel,  this  estimate  of 
the  location  error  Is  not  unreasonable.  The  Interval  (r^  -  r^) 
is  taken  as  0.010  Inches. 


Is 


The  accuracy  with  which 


can  be  determined 


0.01  +  0.01  +  0.03  -  5 $ 


The  error  In  the  quantity  U  is 

k  ^  .  £^£ 

£  r  V+/r 

Again  ^  Is  taken  as  1#.  The  quantity  €v  Is 
0.0005  Inches,  r  *»  1.680. 


The  experimental  error  In  determining 


dr 


This  value  seems  somewhat  low,  but  even  with  the  approximations 
made  It  does  not seem  that  It  could  be  Incorrect  by  any  large 
amount  for  the  region  considered.  This  Is  borne  out  by  the 
results  for  discussed  In  the  next  chapter. 
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7.2.4  The  Error  In  the  Determination  of 


The  quantity  (1  +  ^5^/)/f  is  correlated  with  y1* 
In  the  next  chapter.  It  is  given  by  the  equation 


Since, as  shown  above,  the  fractional  errors  In  the 
determination  of  r  are  small,  the  error  In  [1  +  ]/f  are 

c/&A- 

those  resulting  from  errors  in  the  velocity  and  the  term  r  ^ — 


The  fractional  error  In 


/(0.02)2  +  (0.05)2 


/  A; 

-  0.054 


Is 


5.4* 


7.2.5  The  Error  in  Determination  of  f 


Because  of  the  manner  in  which  the  shear  stress 
at  the  wall  was  determined  (see  Sec.  6.4)  It  is  difficult  to 
specify  by  equation  the  error  In  f .  An  estimate  based  on  the 
results  obtained  Is  556.  The  results  are  presented  In  Chapter 
VIII,  Section  8.4. 
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7*2.6  The  Error  In  Determination  of 

Combining  the  results  of  Seotlone  7*2*4  and 

7.2.5 

/  (Iggf  .//»*# *74ii 

X  * 

Again,  this  la  In  the  region  -  20.  At  y  * ' 
50  the  error  la  presumed  to  be  greater. 

7.3  The  Error  In  Determination  of 

7.3*1  The  Error  In  Determination  of  dT 

dr 

The  error  In  cffar  can  be  eatlaated  In  the  same 
manner  as  the  error  In  dfc/dr. 

Since  the  variation  In  ~p>  across  the  gap  la 
negligible  so  far  as  a  computation  of  T  and  7  concerned 

(e.g.  1/2  Inch  water  above  atmospheric),  dZ df 

Por  three  evenly  spaced  radii,  r^,  rg  and  r^,  an 
approximation  for  d/&  is 

&/m  &£ 

cfrA  r3-n 

The  error  In  ~/f  results  from  errroa  In 
measurement  of  fringe  shift, 4. 
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Prom  equation  (6-8 )c 

Aj>  9  J3./A/4*  /A  A  / A fir/rye 

The  measurement  of  fringe  shift  A  ,  involves  the 
x  and  y  coordinates  on  the  comparator  used  to  evaluate  the 
photographic  plates. 

Taking  as  typical  the  experimental  data  for  5000 
rpm,  a£s£  160R,  the  fringe  shift  between  the  points  3  and  1  Is 
0.021  Inches  where  f  A3  can  be  expressed  as 

,A3  *  Afa-Zt)  +  Sty*-#)  +  corutant 

and  where  x  and  y  are  the  coordinates  measured  on  the  comparator. 
The  values  for  (x^  -  x1)  and  (y^  -  y^  are  0.0143,  and  0.01 
respectively.  Given  x^  and  y^  the  location  of  x^  and  y^  on 
the  comparator  can  be  determined  to  0.0002  Inches  without 
tiring  the  operator. 

The  fractional  error  in  fAj  is  then 

0.0002  .0002 

0.0143  0.010 

-  0.014  +  0.02 

The  measurement  of  A  «  4  Involved  only  one 

Fringe 

coordinate  on  the  comparator,  and  was  taken  over  a  span  of 
about  ten  no -flow  fringes.  The  error  Is  considered  negligible. 
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Therefore 

$ &&/LT  -  0.034 

The  error  in  r  -  r  Is  determined  from  similar 
j  X 

considerations.  However,  since  the  magnification  factor  for 
the  Interferometric  plates  was  about  one -half  the  error  In 
determination  of  coordinates  X  and  ^  for  the  present 
calculations  Is  then  0.0004  Inches 

tGi-rj/t-n  *  Sfo-xd/Xi’X,  * 

0.0004  0.0004 

“  0.0143  0.010 

-  6.8* 

The  fractional  error  In  c/T/dr  Is  then 

(0.034)2  +  (0.068)2  -  7.6* 

The  error  In  aC  can  be  determined  from  equation 
(2-30)  which  Is  repeated  here. 

{/-jc}  +t)/x 

where 

^  .  fa/rf 

$c/T/d ( Vslr)}^,! 

or  from  computation  directly  from 
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The  latter  la  the  somewhat  simpler  procedure  since 
the  quantities  necessary  have  been  computed. 


From  equation  (2 -26) 


The  fractional  error  ln^  Is 

.S**i  .  Jfc  . 

'+/*  dVdr 


Is 


From  Seotlon  7.1.1 


the  fractional  error 


«  .006  +  .0011  +  .025 
-  .032 


From  the  previous  section 

S  &T/dr} 

~mr  'oor6 

The  fractional  error  ln^t^a.  Is  then 


0.0825 

8-1/^ 


The  error  In  a/  negleotlng  the  error  In 
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determination  of  the  properties  Is _ 

l/ {&*/o()X  =  +($**  /6m)* 


*  f  (0.074 )2  +  (O.O825)2  -  11^ 

7.4  The  Error  In  Determination  of  Taylor  Number,  73 

73  «  T^/fotor  ^rbtord/lt  /?'»&WJ/Z>6>r/6oz 

(&Ta/la)  + 

7.4.1  The  Error  In  Rpm 

The  rpm  as  measured  with  a  chronometric  tachometer 
and  as  kept  under  continuous  observation  with  a  Strobotach  was 
never  observed  to  vary  from  a  given  setting  by  more  than  10 
rpm  over  the  6  second  time  Interval  of  operation  of  the 
ohronometrlc  tachometer.  The  chronometric  tachometer  had 
been  previously  calibrated  and  showed  no  dlscemable  error. 

The  calibration  stand  consisted  of  a  shaft  directly  driven  by 
a  synchronous  motor.  The  fractional  error  In  determination  of 
rpm  Is  then 

10  1 
5ZJTO  "  JGU 

7*4.2  The  Error  In  Determination  of  d 
The  gap  width  d  is 

rstator  “  rrotor 

£  ^ rotor  and  £ r stator  are  taken'  aa  befora'  aa 


0.001  Inches 
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Concentricity  of  the  rotor  and  8tator  were  estab¬ 
lished  during  assembly  with  "go"  and  "no-go"  gages  for  which 
the  "no-go"  gage  was  0.0015  inches  greater  in  diameter  than 
the  "go"  gage.  The  concentricity  was  therefore  within  0.002 
inches  and  Sd  is  taken  as 


Therefore, 


Sd  »  0.002" 

Sd/d  m  0.053 


The  error  in  Taylor  number  is  therefore 

/(S&A S'  if  (0.053) f  +  if  * 

-  6.5# 
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Chapter  VIII 
Experimental  Results 

In  this  Chapter ,  the  experimental  results  are 
presented  and  compared  with  other  Investigators  where  possible. 
The  principal  results  are  as  follows: 

1.  That  turbulent  flow  can  exist  In  the  gap  without 
an  axial  flow  of  air  in  the  range  of  Taylor  number  3000  to 
9000. 


2.  That  the  correlation  offers^  by  Kays  and  BJorklund 

/  ,1/2 

for  the  heat  transfer  in  the  gap,  Nu/NuCond  -  0.175  (Ta) 
can  be  extended  for  low  temperature  differences  from  Taylor 
numbers  of  2000  to  Taylor  numbers  of  9000.  At  higher  temperature 
differences  the  heat  transfer  is  reduced  but  still  varies  as 


(Ta) 


1/2 


3.  a.  That  vs  (j  ^  follows  the  universal  curve 

for  pipe  flow  out  to  ■  10,  and  levels  off. 

b.  That  r  In  the  central  portion  of  the  gap 

Is  not  constant  but  Increases  outward  toward  the  stator. 


4. 


That  the 


ctlon  factor  follows  a  relation 


5,  a.  That  near  the  walls  and  out  to  y4,  ■  50, 
two  curves,  one  for  the  rotor,  and  one  for  the  stator,  suffloe 
to  describe  vs.y*  for  3000  <  Ta  -s.  9000,  and  for 

temperature  differences  of  fro®  60  P  to  330  P» 
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b.  That  the  eddy  diffuslvity  In  the  middle 
of  the  channel  can  be  expressed  by  the  relation 

*  0,227  "7#  ^ 

6.  That  the  ratio  equals  1.4  for  the  condi¬ 

tions  studied. 

and  7.  The  findings  can  be  used  to  determine  friction 

factors  and  overall  heat -transfer,  given  T  .  ,  T  ^  ,  and 

rotor*  stator* 

rpm.  This  Is  done  In  Chapter  IX. 

8.1  Hot-wire  Measurements  to  Establish  That  Turbulent  Plow 
Exists  in  the  Gap 


In  order  to  determine  whether  secondary  flows  existed 
In  the  gap,  axial  traverses  were  made  In  the  center  of  the  gap 
for  two  rotational  speeds;  1200  and  3000  Rpm.  The  results  of 
these  tests  are  shown  In  Pig.  8-1  in  which  are  plotted  7/cwF/ 
for  1200  and  3000  rpm.  Shown  in  Pig.  8-2  are  photographs  of 
oscilloscope  traces  for  these  runs  in  which  differences  in 
the  character  of  the  flows  is  apparent. 


The  data  demonstrate  that  two  kinds  of  flow  can 
exist  in  the  gap,  those  for  1200  rpm  showing  the  existence  of 
secondary  flows,  and  those  for  3000  rpm  showing  the  random 
fluctuations  of  turbulent  motion. 


for 


It  Is  Interesting  that  the  peaks  ln^ 

1200  rpm  occur  at  a  spacing  not  too  different  from  3/8  Inches, 
that  predicted  by  Taylor  for  the  Taylor  vortices. 


The  photographs  for  1200  rpm  are  similar  to  those 


y/f/Ajl 


VN 

o 


o 


o 
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8-2* 

Vortioal  Flow,  1800  RPH 


*  „  V*~*v 


8- 2b 

Random  Turbulence,  3000  Rpm 


Figure  8-2 

Hot-wire  anemometer  traces  at  1200  and  3000  rpm 


< 


Shown  by  Kaye  and  Elgar*  for  turbulent -plua -vortices  flow. 

The  value  for  the  turbulent  intensity 

for  3000  rpm  in  the  center  of  the  gap  was  found  to  be  3.5£. 

This  compares  with  the  value  of  656  found  by  MaoPhall,  quoted 
2 

by  Townsend. 

The  overall  velocity  in  the  center  of  the  channel 
was  checked  with  the  hot  wire  and  agreed  to  wlthwl.5#. 

Further  explorations  made  with  the  probe  at  a  fixed 
axial  position  and  varying  speed  (from  low  to  high  to  low 
values)  Indicated  that  at  about  1500  rpm  a  changeover  occurs 
from  the  secondary  type  of  flow  to  the  turbulent  flow. 

It  appears  that  the  line  of  demarcation  between  the 
turbulent-plus -vortices  and  turbulent  flow  regimes  as  originally 
hypothesized  by  Becker^  (see  Fig.  1-3  and  Section  1.2. 1.1)  is 
correct,  and  turbulent  flow  can  exist  In  the  gap  with  zero 
axial  velocity. 


*J.  Kaye  and  E.  C.  Elgar,  "Modes  of  Adiabatic  and  Diabatic  Fluid 
Flow  in  an  Annulus  with  an  Inner  Rotating  Cyliner,"  ASME  Paper 
No.  57-HT-14. 

2D.  C.  MacPhail,  Turbulence  in  a  Distorted  passage  and  Between 
Rotating  Cylinders,"  Ph.D.  Dissertation,  University  of  Cambridge, 
1941  as  reported  by  A.  A.  Townsend,  The  Structure  of  Turbulent 
Shear  Flow,  p.  303>  Fig.  12-3. 

•%.  M.  Becker,  Sc.D.  Dissertation,  MET,  1957. 
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8.2  Overall  Heat  Transfer 

In  Pig.  8-3  are  shown  the  overall  heat  transfer 

results  for  the  present  experiment.  Drawn  also  Is  the  line 

1/2 

Nu/NUcond  ■  0.175(Ta)  '  which  Kays  and  BJorklund  suggest  as 
representing  their  data  and  those  of  Becker  for  90  Ta  2000. 
The  data  of  Gazley  are  also  shown. 


Gazely's  results  were  obtained  from  Fig.  8  of  his 


paper1  noting  that  Ta 


and  also  that  the  Nusselt 


number  used  there  Is  twice  that  of  the  present  work.  Becker, 
and  Kays  and  BJorklund  have  not  found  any  reason  for  the 
discrepancy  between  Gazley 's  results  and  theirs  although 
Becker  suggests  that  the  fact  that  Gazley  did  not  account  for 
radiation  between  the  two  cyllners  may  be  the  cause  of  the 


discrepancy. 


*0.  Gzley,  Jr.,  "Heat  Transfer  Characteristics  of  the  Rotational 
and  Axial  Plow  Between  Rotating  Cylinders,"  ASME  Paper  No.  56- 
A-128. 

2 

No  data  are  known  concerning  the  emisslvlty  of  the  cylinders 
of  Gazley.  If  an  emisslvlty  of  O.36  Is  assumed  and  also  that 
the  temperature  difference  was  of  the  order  of  100  P,  then 
for  a  run  at  a  Taylor  number  of  600,  neglecting  radiation  would 
account  for  a  discrepancy  of  1(#  as  compared  to  a  dlscrepency 
of  605&  as  reported  (e.g.,  at  Ta  «  600,  the  line  Nu/Nuqq^  ■ 

.175(Ta)V2  yields  a  value  of  about  4.4  for  Nu/Nucon(j  as  compared 
with  7  as  reported  by  Gazley). 

Thus  it  does  not  appear  that  under  the  assumptions  made.  Including 
the  effects  of  radiation  heat  transfer  can  account  for  the  difference. 
It  Is  Interesting  however  that  as  the  Taylor  number  Is  reduced  the 
discrepancy  Is  more  acute,  which  Is  the  effect  one  would  expect 
If  radiation  effects  were  the  cause,  since  the  radiation  heat 
transfer  Is  constant  at  constant  temperature  difference. 


Fy  0-3 

Over*//  Feat  Transfer 

Nu/NuiCond  i/s.  Ta 

o  Trvitr  ■  boo  ^  <l7»  80 
*  TToior  3  /&£  AT*  //S 
+•  Troktr  &  8do0  AT*  3 io  ^ 

Extension  of  Fu/f/u^^  *0.(JsTa 
Trend  of  Fesent  Work  at  higher  A.T 
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As  oan  be  seen  from  Pig.  8-3  the  results  of  the 
present  experiment  follow  the  trend  of  the  curve  Nu/Hueon(i  - 
O.r^Ta)1/2  but  the  effects  of  property  variation  are  clearly 
In  evidence. 

At  low  temperature  differences  the  agreement  with 
the  straight  line  representation  of  Kays  and  BJorklund  is 
good. 


However  at  temperature  differences  of  330  7,  the  heat 
transfer  at  a  given  Taylor  number  Is  reduced  by  30£,  and  at 
a  given  rpm  by  20 $  (e.g.  for  7000  rpm,  at  about  60  7  temperature 
difference,  Btu/hr.  transferred/degree  temperature  difference 
■  3,  whereas  at  about  330  7  temperature  difference  the  value 
Is  2.4) . 


Heat  transfer  balance  between  electrical  Input  and 
heat  output  In  the  cooling  water  was  only  to  within  50£.  Study 
of  the  problem  has  shown  that 

1.  A  negligible  amount  of  heat  Is  leaking  from  the 
rotor  to  the  cooling  water  through  the  stator  end  plates. 

This  Is  an  experimental  fact. 

2.  The  difference  between  electrical  energy  Input 
and  heat  abstracted  from  the  cooling  water  oan  pass  directly 
from  the  rotor  heater  to  the  stainless  steel  shaft.  This  Is 
a  calculated  result. 

Other  Investigators  were  not  subjected  to  the  same 
problem  by  virtue  of  lower  temperature  differences.  Consequently, 


In  their  designs  they  were  able  to  use  materials  of  low  thermal 
conductivity  such  as  plastic  spacers. 

Agreement  with  their  results  Is  noted  In  the  range 
In  which  overlap  occurs  (Tajs  3000). 

8.3  Velocity  Distributions 

8.3.1  V#  /r  versus  r 

Shown  In  figures  8-4  to  8-9  are  typical  velocity  dis¬ 
tributions,  / r  versus  r,  near  the  rotor  and  stator  obtained 

for  a  range  of  rotational  speed  3000  rpm  to  7000  rpm  (2500^  Ta 
^9000)  and  temperature  difference  (  0 ^  (Trotor  “  Tstator)^ 
330  P).  The  distances  from  the  wall  have  not  been  adjusted 
for  the  displacement  effect  discussed  In  Section  6.4.  These 
are  the  data  from  which  As  was  calculated. 

8.3 .2  r  versus  r 

Shown  in  Fig.  8-10  and  8-11  are  typical  distributions 
of  r  versus  r.  As  can  be  seen  In  the  oenter  of  the  gap  r 
appears  to  follow  a  relation  r  ■  w  +  b,  Increasing  with  r. 
The  variation  Is  small  but  does  exist.  This  result  agrees  with 
Taylor's  original  data1  but  not  his  final  result  obtained  from 
extrapolation  of  the  data  to  a  probe  of  zero  thickness  (see 
Section  1.2. 1.3).  Pal's2  experiments  were  conducted  to  test 

*0.  I.  Taylor,  Proc.  Royal  Society  of  London,  ser  A,  vol.  157, 
1936,  pp.  492-512. 

2Pal  Shlh-I,  NACA  TN  No.  892,  p.  2. 
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the  validity  of  Taylor's  extrapolation  but  In  his  report  he 
mentions  this  and  then  makes  no  further  statement  concerning 
the  problem. 
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Pal's  data  however  show  a  variation  of  r  Increasing 
with  r.  See  his  Pig.  7  and  9» 

1  M 

Townsend  states  "A  close  study  of  the  velocity 
distribution  shows  that  Ur [V^r  In  the  present  notation]  Is 
slightly  less  than  Its  central  value  near  the  Inner  boundary 
layer  and  slightly  greater  near  the  outer  one.w 


Accordingly,  the  result  for  the  center  of  the  gap, 
V^r  -  ar  +  b  Is  accepted  as  correct  and  was  used  to  determine 
values  of  In  that  region. 


8.3.3  v/  vs.  <4+ 

Shown  In  Pig*  8-1*21  to  8-l4-are  curves  of  Ve  vs.cJT 
for  the  present  experiment. 


^  t- 

The  close  agreement  to  the  curve  of  4,0  vs .  JJ 

J.  ' 

presented  by  Slelcher  for  the  lower  jjTis  understandable 
(see  Sec.  6-4).  However,  agreement  exists  beyond  the  lower 
out  to  10.  At  high  values  the  behavior  with  heat 

2 

transfer  Is  In  agreement  with  the  trend  predicted  by  Deissler  j 


*A.  A.  Townsend,  The  Structure  of  Turbulent  Shear  Plow  (Cambridge 
at  the  University  Press,  19 56}/  p.  29a. 

2r.  q.  Deissler,  "Investigation  of  Turbulent  Plow  and  Heat 
Transfer  In  Smooth  Tubes,  Including  the  Effects  of  Variable 
Fluid  Properties,"  Trans.  ASME,  vol.  73*  1951*  PP.  101-107. 
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namely,  a  divergence  between  the  curves  for  the  stator  and 
rotor  caused  by  the  effects  of  cooling  on  the  stator  and  heating 
on  the  rotor. 

Wattendorf '^reported  results  for  V q*  vs.  r  over 
the  rotor  only  approach  the  rotor  up  to  a  value  of  y  ^  of 
about  20,  so  no  comparison  in  that  region  is  possible.  Purther 
out  in  the  stream  Wattendorf  notes  a  divergence  in  vj"  for 
the  rotor  and  stator.  His  Pig.  23  is  reproduced  as  Pig.  8-  15. 


Such  a  divergence  does  not  occur  in  the  present 
experiment.  No  satisfactory  explanation  can  be  offered  for 
the  difference.  A  curious  fact,  however,  is  that  on  page  582 
Wattendorf  shows  two  successive  equations  for  the  shear  stress 
distribution  as 

rv  -  constant 


and 

1=1 


The  omission  of  the  power  "2"  on  the  radii  may 

be  a  typographical  error.  But,  if  not  since  Wattendorf 

measured  T  >  his  computation  for  7^  would  be  in 
outer  inner 

error  by  router  ^  25  „  1.25.  Thus .  ^  g would  be 

rinner  '  f3 

low  by  255S.  Multiplying  Wattendorf  *s  # inner  an<* 

his  y^by  1.25  reduces  the  discrepancy. 


1  P.  L.  Wattendorf,  Proc.  Aoyal  Soo.  London,  ser  A,  no.  865, 
vol  •  148,  Peb.  1935  >  PP*  565-598. 
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8.4  Friction  Factor 


The  friction  factors  found  In  this  experiment  have 
been  determined  from  the  values  of  obtained  from  the  velocity 
distributions  as  outlined  In  Section  6.4. 


For  purposes  of  correlation  of  the  G.'tfxs  data  the 
friction  factor  has  been  defined  as 


and  not  as 


r 

**“  x  a* 


Outt 


7is^- 


as  by  Kays  and  BJorklund 


Thus  neglecting  property  variation 

Jl  m  5  2 

-/rotor  Kays  and  BJorklund. 

A  plot  of  /"rotor  Stator  ver8us  Ta  18  Sh0"n 

In  Fig.  8-/£.  A  least-squares  straight  line  has  been  calculated 

for  these  data  and  yields  the  relationship 

J- rotor  -  (8-1) 

Ta°*358 

Apparently,  the  manner  in  which  fluid  properties 


1W.  M.  Kays  and  I.  S.  BJorklund,  ASME  Paper  No.  58-A-99 

o 

The  factor  5  results  from  the  fact  that  V*c  -  0.45  V*  rotor 
for  the  dimensions  of  the  present  apparatus.  * 
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used  In  the  dimensionless  ratios  are  evaluated  Is  beginning 
to  have  some  Influence  at  ■^rotor  “  R* 


In  Fig.  8-/7,  7^rotor  18  Plotted  against  Ta,  with 
properties  based  on  rotor  conditions.  Similarly  r  B^at0r  *s 
plotted  against  Ta  with  properties  based  on  stator  conditions. 
For  the  rotor,  it  is  seen  that  the  evaluation  of  the  Taylor 
number  with  the  kinematic  viscosity  based  on  rotor  temperature 
has  succeded  in  placing  the  points  representing  high  temperature 
difference  on  the  correlation  curve. 


For  the  stator, the  effects  of  property  variation 
on  the  correlation  have  been  reduced,  since  the  stator  was 
always  at  about  the  same  temperature ,  80  to  100  F.  The  obvious 
conclusion  is  that  the  effect  of  heat  transfer  has  been  to 
reduce  the  friction  by  about  20$  for  the  temperature  range 
considered. 


Two  other  efforts  in  the  determination  of  friction 
factors  have  been  by  Wagner^  and  Taylor .2 


Their  data  were  obtained  by  torque  measurements  on 
the  stator,  Wagner  using  oil  as  the  fluid  and  Taylor  using 


^E.  M.  Wagner,  "Frictional  Resistance  of  a  Cylinder  Rotating 
in  a  Viscous  Fluid  Within  a  Coaxial  Cylinder,"  Thesis  for 
Degree  of  Engineer,  Stanford  University,  June  1932. 

2G.  I.  Taylor,  'Fluid  Friction  Between  Rotating  Cylinders." 
I-Torque  Measurements,  Proc.  Roy.  Soo.  (London;  sen  A,  vol.  157 > 
1936,  pp.  546-564. 
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water,  pentane,  glycerine  and  water,  and  aviation  spirit 
(18°C)  [gasolene?]. 

Wagner's  data  have  been  correlated  by  Kays  and 
BJorklund,  who  have  found  the  relation 

o-w**0-5®8  <8-2> 

Agreement  of  Wagner's  data  with  the  correlation  Is 
excellent,  especially  when  viewed  in  the  light  of  the  wide 
variation  of  frotor  which  is  d/rrotor  dePendent*  The 
correlation  is  achieved  by  the  factor  1/  (f  Ta)^am. 

A  comparison  between  the  Days  and  BJorklund  correlation 

for  Wagner's  data  and  the  correlation  for  the  present  data  Is 

shown  in  Pig.  8-18.  Also  shown  are  Taylor's  data  for  d/r  .  * 

rotor 

0.267  and  Wagner's  experimental  points  for  his  runs  d  »  O.OI563  ft, 

d/rrotor  ■  0.0807.  This  was  the  largest  gap  he  used  and 

d/r  .  was  the  closest  to  that  of  the  present  experiment 
rotor 

(for  which  d/r  ^  -  0.227). 

rotor 

The  present  data,  for  d/r  .  ■  0.227,  and  those 

rotor 

of  Taylor  for  d/r  ^  =*  0.267  are  lower  than  the  Kays  and 

rotor 

BJorklund  correlation  of  Wagner's  data. 

However,  this  correlation  does  well  for  Taylor's 

data  for  d/r  „  -  0.0278,  0.0588,  and  0.0108. 

'  rotor 

No  firm  conclusion  can  be  drawn  concerning  the 


195 


adequacy  of  the  correlation,  but  it  appears  that  at  large 
d/r^tor  it;  ls  Inadequate. 


Pig.  8-/0  is  presented  to  show  the  dependence  of 
/rotor  on  d/rrotor*  In  this  figure,  Taylor's  data  for 
d/rrotor  “  0.028,  and  0.267  are  included  together  with  Wagner's 
data  for  d/rrotor  =  0.0808. 


It  should  be  noted  that 

a.  The  circled  points  are  not  Taylor's  original 
data  but  points  taken  at  convenient  places  from  a  smoothed 
curve  of  his  data  as  presented  in  his  Pig.  10,  pg.  560  in  the 
reference  cited,  and 

b.  The  shear  stress  used  to  compute 
in  that  figure  is  the  stator  shear  stress. 


between 

i.U  . 


For  small  valuesfd/rrotor  there  is  little  difference 

^rotor  and  ^stator  but  tor  hls  value  rotor 
This  is  mentioned  because  Gazley*  in  his  Pig.  4, 


"Skin  friction  on  outer  cylinder  for  rotation  of  inner  cylinder 
and  no  axial  flow'^has  apparently  reported  friction  factors 
for  Taylor's  work  based  on  the  shear  stress  at  the  rotor  rather 
than  the  stator. 


8.5  Eddy  Diffusivity  for  Momentum, 


The  quantity  can  be  determined  from  knowledge 


*C.  Gazley,  Jr.,  ASME  Paper  No.  56-A-128,  p.  8. 
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of  the  shear  stress,  temperature  for  determination  of 
and  as  shown  In  equation  (2-20). 

<?*>/&) 

Thus  (1  +^»^)  -  — 

/ -S  0/n 

The  determination  of  %  requires  a  knowledge  of 
the  shear  stress  at  the  wall.  Before  the  method  for  determining 
u*  at  the  walls  was  developed  it  was  decided  to  present  the 
values  in  a  form  which  included  only  experimentally 
determined  quantities.  This  is  accomplished  as  below: 


/s  atr 


Since 


/•  _  /*/ _ 


and 


%r2- 

*.  <2 


0+^A>)/j£  -  — /m#  ) _ L 

3*3*/  ( yU  / 


8.5*1  Grijb/  Near  the  Wall 


c/r 


(8-3) 


In  Fig.  8-20  are  presented  (1  +^j<^)/y7  versus  y*  for 
the  various  rpm  and  temperature  differences.  The  correlation 
is  good  especially  when  it  is  considered  that  the  variation 
of^  for  the  runs  presented  is  from  3S  =>  3.8  x  10“4  at 
890  R  to  ^  «  1.6  x  10“4  for  540  or  2,4  times. 


hence  the  divergence  of 


As  y^°j  —j  j:  9 
the  lower  ends  of  the  curves  In  Pig.  8 -20 1 


In  an  attempt  to  arrive  at  a  universal  correlation 
It  was  decided  to  correlate  These  correlations 

are  shown  in  Pig.  8-2/. 

In  Fig.  8-22,  are  shown  the  curves  finally  selected 
S3  representing  for  the  range  of  speed  3000  rpm  to 

7000  rpm  and  the  range  of  temperature  of  the  rotor  890  R  to 
600  R. 


8. 5. 1.1  Comparison  with  other  Investigators 

8. 5. 1.1.1  Experimental  work  of  Slelcher 

Por  purposes  of  comparison  with  the  work  of 
Slelcher,  values  of  €#  vs.  have  bean  plotted  against 
together  with  the  curve  suggested  by  Slelcher  for  pipe  flow 
at  temperature  differences  0  -  20  P  and  a  Reynold rs  number 
range  of  14,200  to  80,500.  This  Is  shown  In  Pig.  8-23. 

The  eddy  diffusivlty  for  heat,  ^  has  been  plotted 
rather  than  ^tf  since  It  was  this  quantity  that  Slelcher 
evaluated  from  his  experimentally  determined  temperature 
distributions.  For  the  present  date,  use  was  made  of  the 
experimental  fact,  not  yet  discussed,  that  *  *»/&,  .  1.4. 

Values  of  were  then  multiplied  by  1.4  to  yield  , 
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4 


as 


Shown  dashed  In  the  figure  are  lines  representing 
the  range  of  Sleicher's  experimental  points.  The  lower 
dashed  line  results  mostly  from  data  from  two  runs,  those 
at  Re  ■  80,500,  the  highest  experimental  value  of  Re  at 
which  Slelcher  reports  data.  The  Reynolds  number  for  his 
other  runs  varied  from  14,000  to  about  38,000,  and  these 
data  are  closer  to  his  suggested  curve. 


Referring  to  the  figure,  note  that  a  correlation 
for  the  present  data  of  is  not  as  effective 

as  the  one  finally  selected,  of  vs.  # 

If  Sleicher's  data  are  presented  on  an  vs.  ^ 

basic  as  in  Pig.  8-24  a  better  correlation  is  obtained  for 
his  data  out  to  y ^”of  about  20.  Shown  also  Is  the  relation 
for  the  present  work  shifted  in  the  6a/M^dlrection 
so  as  to  lie  with  Sleicher's  data.  The  agreement  between  the 
shapes  of  the  curves  in  this  region  is  rather  good. 
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8. 5. 1.1. 2  Comparison  With  Formulations  of  Delssler  and  Van 


For  the  value  of  eddy  dlffusslvlty  near  the  wall 

1  *  . 

Delssler  has  suggested  v  . 


u  <sT 


(8-4) 


2 

and  Van  Driest  has  suggested 

where  and  A*  are  experimentally  determined 
constants  with  the  values  K  -  0.4,  and  A  -  26. 

These  relations  are  plotted  In  Fig.  8 -25”  together 
with  the  range  of  values  for  the  present  experiment.  That  Is, 
the  universal  curve  of  Fig.  8-22.  has  been  multiplied  by  the 
range  of  friction  factors  encountered  during  the  experiment, 
the  higher  ones  for  3000  rpm,  and  the  lower  ones  for  7000  rpm. 

It  Is  seen  that  the  correlation  used, 
results  In  a  variation  In  eddy  diffuslvity  for  momentum  with 
Reynolds  number.  Since  the  original  velocity  data  near  the 
wall  were  placed  on  a  single  curve  of  vs.  a  small 

error  results  in  the  use  of  the  present  correlation.  For  the 


^■C.  C.  Lin,  Turbulent  Flows  and  Heat  Transfer  (Princeton 
University  Press,  Princeton,  New  Jersey,  1959)*  Section  E,  by 
R.  G.  Delssler,  p.  230. 

^E.  R.  Van  Driest,  "On  Turbulent  Flow  Near  a  Wall,"  Heat 
Transient  and  Fluid  Mechanics  Institute  Preprints  of  Papers, 
Paper  12,  University  of  California,  1955. 


7}  -  y,  (8-5) 
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range  of  friction  factors  found  In  the  present  experiment  the 
error  in  is  small  at  low  being  of  the  order  of  i,  10$<. 

y* 

This  is  in  a  region  of  y  where  the  eddy  dlffuslvlty  is  very  . 
small.  Therefore,  small  error  is  made  in  the  computations 
for  overall  heat-transfer  and  friction  which  follows. 

It  is  interesting  to  note  that  Sleicher's  experiment¬ 
al  data  show  a  Reynolds  number  dependence  close  to  the  wall. 

8.5*2  In  the  Center  of  the  Gap 

Because  the  velocities  in  the  center  of  the  gap 
were  not  determined  at  sufficiently  close  Intervals  so  that 
eddy  dlffuslvltles  for  momentum  in  the  center  of  the  gap  could 
be  determined  by  the  same  methods  as  in  Section  8. 5. 1.1  these 
quantities  were  determined  as  below: 

was  almost,  but  not  quite  constant  in  the 
middle  of  the  gap.  As  stated  in  Section  8.3*2;  was  found 

to  be  a  function  of  r  according  to 

V^.  r  m  ar  +  o 

Values  of  vs.  r  for  5000  rpm  are  given  in 

Table  “WL-f. 

Table  VIII-I 

in  the  Central  Region  Expressed  as  a  Function 

of  r,  for  5000  rpm,  0.146^  r  O.I56 
rpm  /  t 

5000  0  r  -  24.9066r-  +  1.1494 

60  r  -  17*9372r  +  2.1223 

160  r  -  I6.8492r  +  2.34945 

360  r  -  I8.9400r  +  2.0811 


I 
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Prom  these  relations  can  be  found  and  hence 

dr 

Typical  values  f or  vs.  r  for  various  rpm  are 

given  in  Table  WT-jT 


In  order  to  develop  a  relation  for  in  the 

center  of  the  gap  which  did  not  depend  on  a  knowledge  of  the 
variation  of  with  r^  was  assumed  constant  in  that 

region. 

Thus 

-  A 

and  ... 

_ 2/4 

dr  ~  r^ 


and  in  the  center  of  the  gap 

=  -2^c/fc  _  _  2  ify c 

Using  the  following  £ 

2/Z 

Ta  - 
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and  assuming 


The  value  of  GrjJzJ  obtained  from  this  equation  is 
incorrect  because  even  though  the  variation  in  V^r  is  small 
an  error  amounting  to  approximately  60#  is  introduced  by 
assuming  V^r  =  constant  rather  than  the  correct  formulation 
Var  a  ar  +  b.  Accordingly,  since  the  correct  values  were 
known  the  equation  above  was  multiplied  by  a  constant  factor 
adjusted  to  give  the  correct  value  of  in  the  center  of 

the  gap.  This  works  rather  well  and  does  so  because  the  product 
of  rpm  and  the  slope  of  the  curve  V^r  vs.  r  is  about  constant. 

The  final  result  for  the  eddy  diffusivity  in  the 
central  portion  of  the  gap,  for  1.75  <  ** "<  1.95 »  is 

0  (8"6a  ) 

Making  use  of  the  relationship  between  fpotor  and  Ta, 

0.33^2 

frotor  -  Ta0.358 

and  the  dimensions  of  the  present  apparatus,  there  results 

(1  =  0.227  Ta0*642  [1  +^£tator  j  (g_6b) 

1.75  <  r"  ^  1.95  1  rotor 

A  comparison  of  this  formulation  with  the  values  obtained  from 

the  velocity  distribution  is  shown  in  Table  VIII-II. 
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Table  8-2 

Values  of  £*/x/  In  Center  of  Gap 


Rpm 

Value  in  center  of  gap 

computed  from 

V^/*  •  ar  +b 

Computed  from 

eq.  (8-6 )b 

5000 

0 

138 

124.3 

6o 

113.5 

113.8 

16  0 

86 

94.6 

360 

75 

75.9 

7000 

0 

162 

157.1 

60 

360 

141.5 

141.4 

Agreement  Is  acceptable. 

8.6  Temperature  Measurements 
8.6.1  T  vs.  r 

Representative  curves  of  T  vs.  r  are  shown  In 
Pig.  8--Z£. 

Because  refractive  effects  over  the  cooled  stator 
resulted  in  an  interference  similar  to  Lloyd's  mirror  inter¬ 
ference  (see  Sec.  5.4.1)  data  could  not  be  obtained  near  the 
rotor.  This  interference  effect  is  clearly  shown  in  Pig.  8-Z7*-, 
which  is  an  interferogram  from  a  typical  run. 


Further  although  the  glass  correction  was  small  for 
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Figure  8-27*5 
((park  Interferogram 


Figure  8-27a 

Typical  Interf erogram 
showing 

Refractive  Effeote  of  Cooled  gtator 
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small  change  In  r,  the  cumulative  effect  of  errors  In  the 
correction,  which  was  applied  starting  at  the  rotor,  over  the 
entire  gap  width  Is  clearly  In  evidence. 

Also  shown  Is  an  lnterferogram  taken  with  a  single 
flash  of  the  mercury  vapor  lamp.  This  technique  was  used  at 
high  speeds  and  density  gradients  because  rotor  run-out  of 
0.0015"  caused  blurring  of  the  fringes  during  exposure  times 
required  for  continuous  operation  of  the  lamp. 

8.6.2  T  vs.  V  /r;  The  Determination  of  o< 

Of  Interest  Is  the  variation  of  T  with  /r.  Pig. 
8-28  shows  the  variation  of  T  versus  Vj  /r  In  the  gap  for 
several  runs. 


The  graphs  were  prepared  In  the  following  manner: 


On  the  rotor  surface  Trotor  is  known  from  thermo¬ 
couple  measurements,  and  V^/r  from  the  speed  of  the  rotor. 

The  straight  line  drawn  Is  (  dJ  Jd^e/r  )wan  and  Is 
obtained  from  (  dT/jr  )*all*  known  from  overall  heat  transfer 
measurements  and  {dltlrjdr  )waii  known  through  .  The  straight 
line  drawn  is  Independent  of  Interferometer  measurements  used 
to  obtain  T  and  total  pressure  measurements  used  to  obtain 

V,/r. 

The  converse  however  Is  not  true  so  far  as  temperature 
is  concerned  because  of  the  manner  In  which  the  temperature  of 
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the  first  measured  point  on  the  interferogram  was  determined 
(see  Section  6.3*1.2).  This  Is  considered  to  have  a  negligible 
effect  on  the  results  for  two  reasons.  First,  although  the 
temperature  distribution  was  "tied"  to  at  the  wall  the 

first  measured  velocity  occurred  at  a  distance  of  the  order 

Ifr  -  5  to  8. 

On  the  Interferograms  however,  data  Is  obtained 
somewhat  closer  to  the  wall  than  this  so  that  the  first  point 
shown  on  the  graph  T  vs.  /r  over  the  rotor  Is  separated 
from  the  first  measured  Interferometric  point  by  two  or  three 
measurements  of  temperature.  Influence  not  only  by  (  )¥in 

but  also  by  the  Intervening  fringe  shifts. 

Secondly,  the  distribution  measured  from  an  Inter¬ 
ferogram  Is  vs.  r.  An  Incorrect  choice  of  a  reference 

temperature  will  Influence  the  distribution  of  4  T  vs,  r  by 
a  negligible  amount.  Thus,  If  an  Incorrect  reference  temperature 
Is  chosen  the  level  of  the  T  vs.  V$/r  curve  Is  affected  In 
T,  but  not  the  shape  of  the  curve.  As  can  be  seen  In  Fig. 

8-28,  T  has  a  straight  line  relationship  with  /r  near  the 
walls. 


From  equations  (2-30)  and  (2-31),  repeated  here 


(2-30) 

(2-3D 

Is  very 
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nearly  unity1  and  since  In  that  region: 


Prom  equation  2-30,  for  x  -  1 

tfot  «r/ 

/fa 

and  for  air 

1.4. 

< 


For  example  at  f/r  ■  60,  corresponding  to  i/+  of  the  same 
order  r'-  1.720,  and  (r/r^^f/^/.cuA &°ss}*/.o7 

so  that  for  x  -  1,  In  fact,  this  Is 

the  means  whereby,for -  1.4, the  transition  is  made  from 
the  slope  at  the  rotor  to  that  at  the  stator. 


i 
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Chapter  IX 

Reconstruction  of  Experimental  Data  and  Prediction 
of  Overall  Heat  Transfer  and  Friction  Factors 
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9 . 1  Introduction 

In  this  chapter  It  will  be  demonstrated  that  the 
experimental  data  can  be  reconstructed  using  the  correlation 
for  VS.  r.  and  that  the  correlation  can  be  used  to 

predict  values  of  overall  heat  transfer  and  friction  faotor, 
knowing  Tstator,  Trotor,  rpm  and  f<  -  l A. 


9.2  Reconstruction  of  Experimental  Data  Using  the 
Correlation  vs . 

9.2.1  Derivation  of  Equations 


From  equation  (2-26) 

which  Is  written 
and  equation  (2-27) 

Integrating,  and  taking  88 


Therefore 


(2-26) 


(2-27) 


lSee  Appendix  V 
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Substitution  Into  equation  (9*1)  yields 

'  /— 

—  /+&•/*() 

_/**+ 


(9-2) 


(9-3) 


Prom  equation  (9-2),  given  Tro^or  and  Tstator*  and 
the  variation  of  ^j^rwlth  r,  the  overall  heat  transfer  can  be 
calculated. 


From  equation  (9-3),  knowing  the  same  yields  the 
variation  of  T  with  r. 


Comparisons  can  now  be  made  with  two  experimentally 
determined  quantities;  as  measured  from  overall  heat 

transfer  results  and  T  vs.  r  as  determined  from  the  Inter¬ 


I 


ferometer. 
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The  same  can  be  done  for  the  velocity  distribution.  Prom 
equations  (2-20)  and  (2-21) 

(2-20) 

and 

^  foo//  ~  tfa'Ar  (2-21) 


we  find 

A 


—  — .  *?-' 

-W,  ***rj /■** 

'Sri erf, 


(9-4) 


The  velocity  distribution  can  be  reconstructed  from 
this  equation. 


However,  to  achieve  a  more  tractable  form  for 
computation  purposes,  consider  the  following 

C/T/C 7/*  =  ~~  £/  ty /^r(^  ^9^0 

and 


Defining  a  Reynolds  number  based  on  the  rotor  diameter, 

and  as 

rotor 
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_  £  'Zlo&i 

J/u4*v  yO  u  * 

equation  (9-5)  can  be  put  into  the  following  form 

['-( Gt/J** 

The  term  in  the  bracket  is  a  constant  for  a  given  apparatus 
speed  and  temperature  difference.  The  factor  1.75  however  results 
from  the  variation  of  ^  with  temperature  for  air. 

If  the  term  on  the  right  is  equal  to  unity  then 
the  derived  velocity  distribution  matches  the  experimental. 


I 


''kr  /  df 


/l~£ri  J  ¥■£*•/*>) 


& 


fy/lctfel 


and  noting  that 


r 


9.2.2  Reconstruction  of  Experimental  Data 

The  run  chosen  for  study  was  one  at  5000  Rpm, 

T rotor  *  722,4  R>  Tstator  “  549,4  R»  the  measured  valu®  of 
overall  heat  transfer  was^^«  1250  Btu/hr-Pt2. 

The  value  of  the  friction  factor  determined  from 
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The  same  can  he  done  for  the  velocity  distribution.  Prom 
equations  (2-20)  and  (2-21) 

J ''r—  /// 

and 

^  ^  /*W'  **■  ^W^r- 


(2-20) 


we  find 


(2-21) 


(9-^) 


The  velocity  distribution  can  be  reconstructed  from 
this  equation. 

However,  to  achieve  a  more  tractable  form  for 
computation  purposes,  consider  the  following 

and 


where ^  =  constant 

Defining  a  Reynolds  number  based  on  the  rotor  diameter. 


and  <2J  as 

rotor 
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p 

velocity  measurements  was  frotor  **  1.589  x  10  . 

Prom  this,  Stator  can  be  calculated  as  fstafcor  -  0.805  x  10"2, 

Other  constants  required  are 

r  =  1.6633  inches 
rotor 

rstator  =  2'037  lnches 

rc  =  £..t237..+  ,1.6633  =  1#8502  inches 


3/  rotor  "  2,64  x  10"9  T1*75  =  2.657  x  10"4ft2/sec 
•^stator  55  1  *^4^  x  10  ft2/sec 
Calculated  from  the  information  above  are 


V©, rotor  =  rp"  x  55  x  1#^  “  72 '58  ft/se0 

ir  _  72.58  X  1.6633  __ 

Vc  “  ""S~x  "1.8502 -  =  32.63  ft/sec 


„  _  1 rotor 

w*, rotor  “  - 2 — 


ve,o  *  2'9°7 


from  which 


and 


u  =  2.0698 

*,  stator 


yT,rotor  ’  1*°968y't 


T,  8  ta  tor  =  o.95Jt3y 

The  first  step  in  the  calculation  is  to  determine 

w 
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Prom  the  curves  for^»»/^vs.  y  +  In  Pig.  8-22 
values  of  are  read  off  at  evenly  spaced  intervals  of 

y  +.  A  convenient  interval  may  be  A  y  +  «  5,  for  example. 
These  curves  can  be  used  for  values  of  y  +  up  to  50. 


Two  integrals  can  then  be  evaulated  for  example 
by  Simpsons'  Rule,  as  tf^sc  and 

J fr/fMntyr  J dr//t+ 

The  equation  (8-6 )a  for  GoA)  in  the  center  of  the  channel  Is 

J'*? 

used  for  1.75  <T  r"  <  1.95  to  evaluate  /  / 

/  ar/f/-t  Gn/v)r 


Two  areas  remain  to  be  evaluated;  those  for  (r  at  y  *  ,  =*  50)  < 

stator 

r"  C  1.95  and  1.75  r"  <.  (r  at  yrQtor  *  50).  No  expression 
has  been  derived  for  these  areas;  however,  a  plot  of  the  function 
5»>/^vs.  r,  as  shown  for  the  present  example  in  Pig.  9-1, 
reveals  how  small  the  contribution  these  areas  make  to  the 


total  integral.  To  evaluate  the  unknown  integrals  the  function 
is  plotted  in  the  known  areas  and  for  the  unknown  ones  the 
curve  is  "faired".  The  unknown  areas  are  then  numerically 
evaluated  using  values  read  from  the  curve. 

The  sum  of  these  five  areas  Is  the  integral 
Prom  equation  (9-2),  can  now  be  computed. 

The  variation  of  T  with  r  can  also  be  calculated. 


9-7a 
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r 

/  dr 

The  integral  / can  be 

4- 

desk  calculator  at  each  value  of  y  as  the  ev 
the  total  integral  proceeds. 

Equation  (9-3)  is  then 

-7-/75  — .  /7sr 

'  -  ^ 

-r-  /  7S~  y-  /7r  “  /&*2n. 

This  can  be  put  into  the  following  forms : 


can  be  read  from  a 


as  the  evaluation  of 


/  ^ 


(OVER  THE  STATOR) 


(9-8) 


~%*s4n. 


(9-9) 


(OVER  THE  ROTOR) 


Thus,  T  vs.  r  can  be  evaluated  and  the  integral 
f  Or  computed. 

J  r*~ 

fath,  All  the  terms  in  equation  (9-6)  are  now  known  and 

the  final  check  can  be  made  concerning  across  the  channel. 

For  the  example  chosen  the  calculation  proceeded  as  follows, 
for  the  area  over  the  rotor: 
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Vr 

A 

1  /W+£~A/) 

Cum 

0 

0 

0 

0 

1.6633 

1.6633 

0.6012 

5 

15 

0.238 

5.48 

1.6688 

2.0660 

0.4840 

2.796x10”3 

10 

47 

0.7468 

10.97 

1.6743 

2.9247 

0.3419 

5.2630x10"3 

15 

127 

2.018 

16.45 

1.6797 

5.0693 

0.1973 

20 

294 

4.672 

21.94 

1.6852 

9.5585 

0.1046 

7.52l8xl0"3 

25 

532 

8.453 

27.42 

1.6907 

15.9822 

0.0626 

30 

840 

13.348 

32.92 

1.6962 

24.3370 

0.0411 

8.2460xl0"3 

35 

1170 

18.591 

38.39 

1.7017 

33.3380 

0.03 

40 

1570 

24.947 

43.87 

1.7071 

44.2941 

0.0226 

8.58l6xlO~3 

45 

1980 

31.462 

49.36 

1.7126 

55.9027 

0.0179 

50 

2410 

38.295 

54.84 

1.7181 

65.5127 

0.0148 

8.781x10'3 

Thus 

r=/7'* 

r> 

1 

^  cl r 

==  <T7#/*/o' 

r=/.ff9 


/  dr  __  7/^^  'l 

Similarly  over  the  stator^  /  /%'/■£ ‘  ^  X/° 


region  r  = 

1.75* to  r  - 

1.95, 

( 

1.75 

167.30 

0.0060 

1.85 

176.9 

0.0057 

1.95 

186.4 

0.0054 

/V*f 


f 

\j-V  / 


-  //<£o  x/o 


.3 


and 
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YrO+4*) 


A  plot  ot  /•(  vs •  r  Is  made,  as  In  Fig.  9-1 

for  the  regions  In  which  Is  known.  For  the  regions 

0  „ 

1.989  r  1.95  and  1.7181  r  ^  1.75  the  curve  Is 

faired  as  shown  In  the  figure.  The  areas  usder  these  portions 
of  the  curve  are  0.415  x  10"^  and  0.322  x  10“^  respectively. 


The  total  area  under  the  curve  Is 


^4/-  0.01708  x  104 


and  si  tl  o 

&  .  *  0.01708  x  10  x  12  -  1232  Btu/hr-ft* 

9  1.6633 

p 

The  measured  value  was  1250  Btu/hr-ft  ,  a  difference 

of  3.256. 

The  temperature  distribution  can  now  be  calculated 
over  the  rotor  the  results  are 


r 

TAk>hi 

T 

/S 

1.6633 

0 

1 

722.4 

0.3615 

1.669 

0.1671 

0.9632 

695.8 

0.3589 

1.674 

0.2956 

0.9341 

674.8 

0.3568 

1.685 

0.4224 

0.9047 

653.5 

0.3522 

1.696 

0.4631 

0.8951 

646.6 

0.3477 

1.707 

0.4891 

0.8906 

643.4 

0.3431 

1.718 

0.4931 

0.8880 

641.5 

0.3388 
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Over  the  stator  the  results  are 


r 

r 

2.037 

549.4 

.2410 

2.027 

590.7 

.2434 

8.018 

690.6 

.2456 

2.008 

617.2 

.2480 

1.999 

620.9 

.2503 

1.989 

623.2 

.2528 

1.95 

627.3 

.263 

1.85 

632.7 

.292 

1.75 

638.3 

.327 

A  plot  Is  now  made  of  T  vs.  ,/r  as  In  Pig.  9-2.  The  area 


under  this  curve  can  be  Integrated  numerically  to  give 


The  quantities  needed  to  determine  G  In  equation 
(9-7)  are  all  known.  Recall  that  G  must  be  unity  for  agreement 
between  the  experimental  and  calculated  velocity  distributions. 


Substitution  Into  equation  (9-7)  results  In 
? 

1  *  G  ■  0.9938 

a  difference  of  0.62# 


In  figure  9-3  Is  shown  the  temperature  distribution 
over  the  rotor,  as  measured  and  as  computed. 
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9.3  The  Use  of  the  Correlation  of  ^*^/v s.  y ^ ,  to  Predict 
Overall  Heat  Transfer  and  Friction  Factors 


The  purpose  of  this  section  is  to  demonstrate  the 
use  of  the  correlation  for  vs.  y+  and  the  fact  that 

^  =  1.4  to  solve  the  problem  "Given:  Rotor  Temperature, 
Stator  Temperature  and  rpm.  What  is  the  overall  heat  transfer 
and  friction  factor?" 


The  proceedure  is  the  same  as  previously,  except 
that  in  this  case  a  trial  friction  factor  is  used  until  the 
right  hand  side  of  equation  (9-6)  is  equal  to  unity. 


This  calculation  has  been  made  for  the  present 

example.  Friction  factors  frotor  1#2^  x  10  2  and  1.9  x  10 
were  selected  to  obtain  G.  Use  was  made  of  the  result  in  the 


foregoing  section.  The  resulting  values  for  G  with  various 
^rotor  v,ere  f  df  f  d *T  S' 


Sous*  0 

/  7^ 

& 

1.25  x  lo"2 

21.998  x  10" 3 

51.439 

.95190 

1.59  x  10'2 

17.806  X  10~3 

52.162 

.99368 

1.9  x  10“2 

15.132  x  10 

52 .293 

1.0123 

The  graph  of  Fig.  9-^was  then  prepared. 

_2 

The  quantity  Gis  unity  at  fro^«1.67  x  10 
Accordingly,  the  computation  was  made  with  frotor  *  x  10 

The  value  G  was  found  to  be  0.9987  and  is  plotted  in  Fig.  9-4, 
This  is  considered  to  be  of  sufficient  accuracy. 


,'ff  Jl 
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The  value  of  frotor  dlff«rs  from  that  found  from 
the  experimental  velocity  distributions  by  5.1*. 

The  overall  heat  transfer  computed  with 

frotor  ■  1.67*  10“2  Is,  a  difference  of  8 *  from  the  experl- 

2 

mental  value  of  1250  Btu/hr-Ft  . 

The  temperature  distributions  calculated  with  the 
trial  values  of  friction  factors  used  In  the  computation  are 
plotted  In  Fig.  9-3.  Better  agreement  with  experiment  Is 
achieved  with  friction  factor  1.67  x  10"2  than  I.589  x  10"2. 

Similar  calculations  have  been  made  for  5000  rpm 
for  temperature  differences  ^  t  ■  70,  170,  and  320  F. 

The  overall  results  compared  with  the  experimental 
values  Is  shown  below. 

5000  rpm 

"“/"“eona  ^rotor 


rotor*** 

Tetator* 

R  At 

calc. 

exp. 

%  dlff. 

calc. 

exp. 

604.7 

535.7 

69 

12.8 

14,6 

12* 

0.0139 

0.0148 

722.4 

549.4 

173 

12.3 

11.5 

7* 

0.0167 

0.0159 

882.7 

564.3 

318.4 

10.4 

9.7 

7.2* 

0.0195 

0.0175 

The  differences  between  the  calculated  and  experimental  values 
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of  the  heat  transfer  are  about  equal  to  the  percentage  error 
calculated  for  in  Chapter  VIII  for  the  rune  at  fit*  17oFand 
320^’but  is  higher  for  the  lower  temperature  run.  For  this 
run,  agreement  between  the  calculated  and  experimental  Nusselt 
Number  with  the  extrapolated  correlation  of  Kays  and  BJorklund 
is  to  within  +6.5$  for  the  experimental  point  and  -6.5$  for 
the  calculated  point. 

Fig.  9-^* Is  a  plot  of  Nu/Nucond#  vs.  Ta  for  the 
calculated  points  and  the  experimental  points  together  with 
the  line  Nu/Nucond>  -  0.175  Ta1/2, 

The  values  of  the  calculated  friction  faotors  have 
been  plotted  in  Fig.  9-^  against  Taylor  number  together  with 
the  experimental  friction  factors. 

A  trend  toward  higher  friction  factors  with  increased 
heat  transfer  is  now  in  evidence  for  the  lower  temperature 
runs  as  well  as  for  the  runs  at  a  temperature  difference  of 
320  F  (see  Fig.  Q‘/£).  Plotted  also  in  Fig.  9-^  are  the 
experimental  points  for  3000  rpm  and  7000  rpm  at  this  temperature 
difference.  The  experimental  value  of  frotor  for  5000  rpm, 

&£  =  320  F,  appeared  to  be  too  low  compared  to  fr0fcor  f°r 
3000  and  7000  rpm.  The  calculated  value  of  frotor  f°r  5000 
rpm  is  11$  higher  than  the  experimental  value,  and  now  appears 
to  be  in  better  relation  to  the  other  two  points. 

The  difference  between  the  calculated  and  experimental 
friction  factors  for  the  lower  temperature  runs  is  within  5-1/2$. 
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The  conclusions  drawn  from  the  foregoing  are  that 
1.  the  correlation^^  vs.  y*  ,  will  reconstruct  the  experimental 
data  from  which  it  is  derived  and  2.  the  correlation  can  be 
used  to  predict  overall  heat  transfer  and  friction  factors 
for  the  present  apparatus  under  conditions  of  property 
variation  vjhen  the  cylinder  temperatures  and  the  speed  of  the 
inner  cylinder  are  given. 


247 


Appendix  I 

I. 1  Description  of  operation  of  a  Mach-Zehnder  Interferometer 

A  schematic  diagram  of  a  Mach  Zehnder  Interferometer 
Is  shown  In  Fig.  1-1.  Monochromatic  light  from  a  point  source 
Is  placed  at  the  focal  point  of  a  collimating  lens  b,  resulting 
In  a  beam  of  parallel  light.  Beam  splitter  c,  a  half  silvered 
mirror  allows  one  half  the  Incident  light  to  pass  along  path  I 
and  reflects  the  other  half  along  path  II.  Plates  a  and  e  are 
totally  reflecting  plane  mirrors.  The  light  proceeding  along 
path  I  Is  reflected  from  the  mirror  d  and  on  striking  beam 
splitter  f  Is  again  divided.  One  half  Is  discarded,  and  the 
other  half  reunites  with  the  light  which  proceeded  along  path 

II. 


If  one  of  the  mirrors,  say  d.  Is  rotated  an  angle 
about  an  axis  perpendicular  to  the  plane  of  the  paper  by  the 
laws  of  Incidence  and  reflection  those  rays  from  path  I  will 
intersect  those  of  path  II  at  an  angle  2  f  ,  as  shown  In 
Fig.  1-2.  In  this  figure  the  straight  lines  drawn  beyond  the 
splitting  plate  f  mark  the  wave  crests  of  the  Individual  rays. 
Each  point  of  Intersection  between  the  lines  for  the  light 
from  path  I  and  those  for  path  II  marks  a  line  perpendicular 
to  the  plane  of  the  paper  along  which  the  rays  from  each  path 
are  In  phase. 

On  a  screen  placed  beyond  spllttlng-plate  f ,  the 
line  of  intersection,  appears  as  a  bright  band,  and  On.  Inter¬ 
ference  pattern  Is  produced  as  shown.  If  a  region  of  uniform 
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density  different  from  that  outside  the  region  Is  plaoed  In 
path  II  a  light  ray  along  that  path  will  be  accelerated  or 
retarded  depending  on  whether  the  air  In  the  Interposed  region 
Is  less  dense  or  more  dense  than  that  outside  the  region. 

This  results  In  a  shift  of  the  wave  front  for  these 
rays  by  an  amount  and  a  shift  of  the  Intersecting  lines 
of  the  wave  fronts  for  the  two  paths  by  an  amount  feS 
That  is,  there  would  be  a  band  shift  on  the  screen  an  amount 
CXtfzS  .  This  Is  shown  In  Pig.  1-3. 

In  general,  during  flow  the  density  in  the  test 
section  will  not  be  uniform  In  a  plane  perpendicular  to  the 
light  path.  Consequently  the  band  shifts  will  be  non-uniform, 
and  the  density  field  must  be  evaluated  on  a  point  by  point 
basis. 

1.2  Relationship  Between  Density  and  Pringe  Shift 

The  relationship  between  density  and  fringe  shift 
is  developed  as  follows: 

Let  t^  be  the  time  for  a  light  ray  to  traverse  the 
length  L  of  the  test  section  during  no-flow  conditions  and 
tg  the  same  during  flow  conditions. 

If  Ve  is  the  velocity  of  light  at  room  air  conditions, 
and  A»  the  wavelength  at  room  air  conditions,  the  displacement 
of  the  wave  fronts  for  the  light  through  the  test  section  is 

^  C't’i-'t- ,) 

where  £  is  the  displacement  measured  In  wavelengths. 


w 

/  *4“ & 
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The  difference  in  time  to  traverse  the  test  section 
between  flow  and  no-flow  conditions  is 

‘2  -  *i-  t/ff-  --fc-) 

The  index  of  refraction  n  is  defined  as 

n  m  Vvacuum 
V 

and  the  frequency  of  the  light,  a  constant,  is  given  as 

'■t 

Therefore 

^  vacuum  “  L<n2  "  nl> 

Note  that  is  the  difference  in  optical  path  length 

vacuum 

between  conditions  1  and  2  and  can  be  Interpreted  not  only  as 
the  change  in  optical  path  length  between  no -flow  and  flow 
conditions,  but  also  as  the  difference  in  optical  path  length 
between  any  two  points  1,  and  2  under  flow  conditions. 

Prom  the  Gladstone -Dale  Equation 

n  -  1  -  kp 

where  M  is  constant  for  a  given  gas  and  wavelength  of  light 
over  a  large  range  of  pressures.* 


1  Values  of  can  be  obtained  from  GL.  M.  Edelman  and  M.  H. 

Bright,  "The  Specific  Refractivity  of  Gases  for  Various  Wave¬ 
lengths  of  Light,"  Gas  Turbine  Laboratory  Report  No.  6, 
Cambridge  Massachusetts:  MIT,  May  1948. 
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Thus  ■  *T,  +l 

and  rig  » 

Therefore 

*  A..0UU.  - 


Knowing  ^  allows  determination  of  the  density  field, 

1.3  Adjustment  of  the  Interferometer 

Adjustment  of  the  Interferometer  Is  aooompllshed  In 
the  following  steps: 

1.  Obtaining  a  parallel  beam  of  light  Into  the 
Interferometer 

2.  Adjusting  the  plates  and  mirrors  to  obtain 
Interference  fringes 

3.  Obtaining  Interference  fringes  of  maximum 
contrast  In  focus  at  the  test  section. 

For  alignment  of  the  Interferometer  with  the 
test  section,  see  Section 

1,3.1  Parallel  light  Into  the  Interferometer 

In  order  to  obtain  a  parallel  beam  of  light  Into 
the  Interferometer  the  small  opening  serving  as  the  effective 
source  must  be  placed  In  the  focal  plane  of  the  collimating 
lens.  This  Is  accomplished  by  focusing  a  telescope  on  a 
distant  object,  sighting  through  the  collimating  lens  with 
the  telescope  at  the  small  opening  whose  position  Is  adjusted 
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along  the  lens  axis  until  it  is  in  sharp  focus. 

1.3.2  Adjusting  the  plates  and  mirrors  to  obtain  interference 
fringes . 

For  the  initial  adjustment  of  the  Interferometer 

1  2 

use  was  made  of  the  method  of  Price.  *  The  method  results 
in  visible  interference  fringes  with  monochromatic  light. 

This  phase  of  the  adjustment  could  be  the  most  time  consuming 
but  for  this  method  which  is  described  below. 

Consider  the  M-Z  arrangement  shown  in  Fig.  1-4. 

For  fringe  formation  the  plates  and  mirrors  must  be  very  nearly 
parallel  and  the  two  available  path  lengths  different  by  not 
more  than  a  few  wavelengths  of  light. 

To  accomplish  this,  place  a  small  plane  mirror  at 
B,  and  a  pinhole  light  at  C. 

A  viewer  at  the  telescope  lens  sees  four  primary 
images  of  the  pinhole  light,  one  for  each  of  the  four  possible 
light  paths  that  the  light  from  the  pinhole  can  traverse  in 
proceeding  from  A  to  B  and  back  to  the  telescope  14ns. 

To  accomplish  the  adjustment  one  rotates  the  mirrors 
1E.  Price,  NAVORD  Report  No.  1182  (1949). 

2E.  Price,  Review  of  Scientific  Instruments,  vol.  23,  Number  4, 
April  1951,  P*  162, 
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and  plates  until  the  four  Inages  are  superimposed. 

One  complication  arises  because  of  secondary  Images 
resulting  from  multiple  reflections  In  the  splitting  plates. 

In  the  present  apparatus  two  secondary  Images  were  seen  for 
each  primary  Image.  The  secondary  Images  unfortunately  do  not 
differ  greatly  In  brightness  from  the  primary  Images  but  one 
finds  In  rotating  the  splitting  plates  that  the  Images  move 
In  sets. 


Once  the  superposition  Is  accomplished,  filtered 
light  from  the  mercury  vapor  lamp  Is  passed  through  the  Inter¬ 
ferometer.  If  no  fringes  are  visible  with  only  the  eye,  rotation 
of  one  of  the  plates  and/or  motion  of  the  translation  mount 
will  usually  result  In  visible  fringes. 

1.3.3  Fringes  of  Maximum  Contrast 

The  fringes  now  visible  are  probably  not  of  maximum 
contrast.  If  the  filter  Is  removed  the  fringes  may  disappear. 
Motion  of  the  translation  mount  will  bring  Into  the  field  the 
few  fringes  about  the  central  band. 

To  focus  the  fringes  In  the  plane  of  Interest  In  the 
test  section,  an  ordinary  60  watt  light  bulb  can  be  used  In 
place  of  the  mercury  vapor  lamp. 

Focus  a  telescope  In  the  plane  of  Interest.  Tilt 
the  splitting  plates  until  the  fringes  come  Into  focus  In  the 
telescope.  When  the  bulb  is  replaced  by  the  mercury  lamp,  the 


fringes  now  appearing  are  of  maximum  contrast.  The  inter¬ 
ferometer  is  ready  for  use. 
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Appendix  II 

Experimental  Determination  of  Temperature 
Gradient  In  the  Windows  at  the  Ends  of  the  Gap 

To  determine  the  temperature  gradient  In  the  windows, 
narror,  thin  stlps  of  Templlaq  of  different  temperature  cali¬ 
brations  were  applied  to  the  quartz  windows  as  shown  In  Pig. 
II-l. 


Templlaq  was  selected  as  the  most  practical  means 
for  determining  the  temperature  gradient.  There  were  two  other 
alternative  available.  One  was  drilling  holes  Into  the  windows 
and  applying  thermocouples  at  selected  points.  This  was  rejected 
as  impractical  and  of  questionable  accuracy.  A  second  was 
applying  thermocouples  at  the  contact  surface  between  the 
quartz  window  and  the  end  plates.  This  was  rejected  because 
only  two  temperatures  on  the  glass  would  be  known,  assuming 
the  temperatures  determined  were  of  sufficient  accuracy. 

Templlaq  Is  a  solution  of  a  substance  whose  melting 
point  is  known  and  a  solvent,  the  compound  being  different 
for  different  temperature  ratings.  When  the  temperature  of 
the  Templlaq  reaches  the  melting  point  of  the  compound,  melting 
occurs  and  on  drying  a  sharp  line  can  be  seen  between  the 
unaffected  Templlaq  and  that  which  has  melted  and  resolidified 
on  cooling. 


In  general,  the  applied  Templlaq  Is  slightly  opaque, 
whereas  the  melted  Templlaq  is  transparent.  There  is  thus  no 
mistaking  the  line  of  demarcation. 
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The  apparatus  was  operated  exactly  as  during  an 
Interferometric  run  and  the  pattern  of  the  Isotherms  emerged 
on  the  windows  as  circles  concentric  with  the  rotor. 

Measurements  were  made  on  the  Jones  and  Lamson 
Pedestal  Comparator  to  determine  the  gradient  but  the  size  of 
the  correction  did  not  warrant  that  kind  of  accuracy  so 
ultimately  the  isotherms  were  traced  onto  a  piece  of  graph 
paper  and  the  gradients  determined. 

A  diagram  of  the  isotherms  for  a  typical  result  is 
shown  in  Pig.  II-/. 

In  Pig.  II-2  is  shown  the  temperature  difference  in 
the  glass  over  the  gap  width  versus  rpra  at  different  temperature 
differences  across  the  gap. 

Listed  below  are  the  experimental  results 
Table  II-l 

Degrees  p/inch  versus  Rpm  at  Different  Trotor-Tstator 


60 

160 

360 

RPM 

3000 

51.2 

90 

185 

5000 

57.6 

101.4 

187 

7000 

59.2 

117.5 

202 

A  simple  fin  analysis  on  the  window  shows  that 
convective  effects  on  the  surface  of  the  glass  play  only  a 
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small  part  in  determining  the  gradient  in  the  glass. 


< 


< 


i 


261 


868 


Appendix  III 

Experiments  Performed  to  Determine  the  Accuracy 
of  the  Slip -Ring  Assembly 

The  sequence  of  tests  performed  to  determine  the 
accuracy  of  indication  of  the  slip-ring  assembly  was  as 
follows  s 

Test  #1.  To  determine  the  error  caused  by  frictional 
heating  of  the  rings  when  the  brushes  were  in  continuous  con¬ 
tact  with  the  rings. 

Test  #2.  To  determine  the  reduction  in  error  when 
cooling  air  is  passed  over  the  slip  rings. 

Test  #3.  To  determine  the  reduction  in  error  when 
cooper  lead  wire  is  introduced  between  the  rotor  thermocouples 
and  the  slip  rings. 

Test  #4.  To  determine  the  slip  ring  error  with 
the  circuit  arranged  as  in  Test  #3,  with  cooling  air  passing 
over  the  slip  rings,  with  the  brushes  in  momentary  contact 
during  a  reading  and  with  current  in  the  heater. 

This  comprises  the  final  slip-ring  arrangement.  A 
description  of  these  tests  together  with  the  results  obtained 
follows. 

Test  #1. 

a.  Arrangement  as  shown  in  Pig.  4-2.  except  that 


iron-constantan  leads  from  the  rotor  thermocouples  are  soldered 
directly  to  the  silver  slip  rings. 
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b.  Brushes  in  continuous  contact  with  the 

si ip -rings. 

c.  No  heater  current. 

With  no  heater  current  the  error  in  eraf  is  the 
difference  between  the  emf  indicated  with  no  rotation  and  the 
eraf  indicated  with  rotation  (since  frictional  heating  effects 
on  the  rotor  are  negligible). 

Result:.  The  slip-ring  error  after  steady  state 
conditions  are  reached  at  4000  and  5000  rp»  is  of  the  order 
of  60  P. 


Test  #2. 

Conditions  the  same  as  in  Test  #1  except  that 
cooling  air  was  passed  over  the  slip  rings.  The  effect  of 
passing  compressed  air  over  the  slip-rings  was  to  reduce  the 
error  of  indication  to  an  order  of  10  P  for  the  thermocouple 
slip  rings  closest  to  the  heater  slip  rings. 

The  errors  resulting  were  considered  to  be  too 
high,  and  the  slip  ring  assembly  was  altered  so  as  to  interpose 
copper  wire  between  the  rotor  thermocouples  and  the  silver  slip 
rings.  Only  three  thermocouples  were  connected  in  this  manner. 
The  copper  wires  issuing  from  the  two  slip-rings  closest  to 
the  heater  were  Joined.  Any  indication  from  this  circuit 
was  due  mainly  to  slip  ring  error. 

Two  further  tests  were  then  conducted. 
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Test  #3. 

This  test  was  conducted  under  exactly  the  sane 
conditions  as  Test  #1,  except  for  the  copper  leads  between 
the  rotor  thermocouples  and  the  slip-rings.  See  Fig.  4—2. 

With  the  brushes  in  continuous  contact,  no  current 
in  the  heater,  and  no  cooling  air  over  the  slip-rings,  the 
errors  in  indication  were  of  the  order  of  1  to  2  deg  F  for 
speeds  in  steps  of  1000  rpm  up  to  6000  rpm  for  20  minutes  at 
each  speed. 

Test  #4. 

This  test  was  conducted  with  the  slip  ring 
assembly  in  its  final  arrangement. 

Cooling  air  was  passed  over  the  rings,  the  brushes 
were  only  in  intermittent  contact  with  the  rings,  and  heater 
current  was  passed  through  the  rings  as  during  the  subsequent 
experiments . 


The  reading  obtained  from  the  set  of  rings  adjacent 
to  theheater  rings.  Which  were  joined  as  in  Test  #3#  should 
have  read  in  the  absence  of  ring  error,  the  temperature  of  the 
compensating  Junction  in  the  circuit. 

An  independent  thermocouple  was  placed  near  the 
stationary,  compensating  thermocouple. 

The  indication  from  the  independent  thermocouple 
agreed  with  the  indication  from  the  two  rings  adjacent  to  the 
heater  rings  to  within  3  F  at  conditions  approximating  those 
met  during  the  subsequent  experiment. 


The  slip  ring  assembly  is  thus  considered  to  be 
capable  of  indicating  rotor  temperature  to  within  3  P  under 
conditions  of  the  experiment. 


Overall  Heat  Transfer 
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Appendix  IV  B 
Velocity  Distributions 


(0.  versus  r 


Run  54  V 

3000  rpm 

Adiabatic, 

t  »  8o  F 

Stator 

r 

,  inches 

/^ft/sc 

0 

2.037 

0 

5 

2.032 

6.68 

9 

2.028 

9.85 

12 

2.025 

11.67 

15 

2.022 

13.94 

20 

2.017 

16.08 

30 

2.007 

18.23 

40 

1.997 

19.02 

65 

1.972 

18.56 

90 

1.947 

18.30 

14  0 

1.897 

18.37 

Rotor 

v,  inches 

^  ft/sec 

0 

1.660 

(43.45) 

8 

1.668 

34.60 

10 

1.670 

32.11 

15 

1.675 

28.13 

20 

1.680 

25.94 

25 

1.685 

24.53 

35 

1.695 

23.14 

45 

1.705 

22.19 

60 

1.720 

21.40 

85 

1.745 

20.74 

185 

1.845 

18.89 
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Appendix  IV  B 
Velocity  Distributions 
versus  r 


Run  55  V 

3000  rpm 

^ rot or  “ 

606.1  R 

T 

stator  " 

>  538.2  R 

AT  - 

■  67.9  R 

Stator 

r,  inches 

ft/sec 

0 

2.037 

0 

5 

2.032 

6.68 

9 

2.028 

9.83 

12 

2.025 

11.41 

15 

2,022 

13.51 

20 

2.017 

15.75 

30 

2.007 

18.00 

40 

1.997 

18.63 

65 

1.972 

19.12 

90 

1.947 

18.90 

14  0 

1.897 

18.52 

Rotor 

r,  inches 

I/s,  ft/sec 

0 

1.6613 

43.49 

8 

1.6693 

34.91 

10 

1.6713 

32.72 

15 

1.6763 

28.56 

19.5 

1.6808 

26.38 

24.5 

I.6858 

24.85 

34.5 

1.6958 

23.24 

44.5 

1.7058 

22.11 

59.5 

1.7208 

21.51 

84.5 

1.7458 

20.69 

184.5 

1.8458 

18.77 

* 
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Appendix  IV  B 
Velocity  Distributions 
^  versus  r 


Run  58  V 

3000  rpm 

^rotor  * 

725.73  R 

Tstator  ■ 

546.22  R 

A  T  - 

179.51  R 

Stator 

r,  inches 

ft/sec 

0 

2.037 

0 

5 

2.032 

5.94 

9 

2.028 

8.26 

12 

2.025 

9.91 

15 

2.022 

11.15 

20 

2.017 

12.35 

30 

2.007 

14.01 

40 

1.997 

15.31 

65 

1.972 

17.63 

90 

1.047 

18.11 

140 

1.897 

18.48 

Rotor 

fy/T) 

r,  inches 

lb  ft/sec 

0 

1.6634 

(43.54 

8 

1.6714 

36.22 

10 

1.6734 

34.61 

15 

1.6784 

29.91 

20 

1.6834 

27.28 

25 

1.6884 

25.90  * 

35 

1.6984 

23.92 

1.7084 

22.86 

60 

1.7234 

21.72 

85 

1.7484 

20.54 

185 

1.8484 

18.73 
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Appendix  IV  B 


Velocity  Distributions 
)/q  versus  r 


Run  59  V 

3000  rpra 

Trotor  -  847.9  R 

Tstator  '  555.4  S 

/J  T  .  292.5  R 

Rotor 

!Jrno 

r.  Inches 

(4  ft/sec 

0 

1.6654 

43.6 

8 

1.6734 

38.71 

10 

1.6754 

36.27 

15 

1.6804 

31.87 

20 

1.6854 

28.81 

25 

1.6904 

26.83 

35 

1.7004 

24.37 

45 

1.7104 

22.86 

60 

1.7254 

21.66 

85 

1.7504 

19.90 

185 

1.8504 

17.55 

273 


Appendix  IV  B 


Velocity  Distributions 
versus  r 


Run  46  V 

5000  rpm 

Adiabatic, 

Stator 

J h*  r 

t  =  81  P 

,  inches 

\^ft/sec 

Rotor 

r,  inches 

14  ft/sec 

0 

2.037 

0 

0 

1.660 

(72.4) 

5 

2.032 

13.5 

8.3 

1.6683 

51.73 

9 

2.028 

19.25 

10.3 

1.6703 

47.57 

12.4 

2.0246 

23.56 

15.3 

1.6753 

41.18 

17.9 

2.019 

26.53 

20.3 

1.6803 

38.38 

27.9 

2.009 

29.28 

25.3 

1.6853 

37.32 

37.4 

1.9996 

30.73 

35.3 

1.6953 

35.58 

62.4 

1.9746 

31.58 

45.3 

1.7053 

34.79 

86.4 

1.9506 

32.03 

60.3 

1.7203 

34.06 

136.4 

1.9006 

32.18 

85.3 

1.7453 

32.96 

188.9 

1.8481 

32.09 

185.3 

1.8453 

32.11 
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Run  47  V 
5000  rpm 


T 

rotor 

=»  606.2  R 

^•stator 

=  526.4  R 

A  T 

Stator 

-  79.8  R 

<k» 

r,  inches 

0 

2.037 

5 

2.032 

9 

2.028 

12 

2.025 

15 

2.022 

20 

2.017 

30 

2.007 

40 

1.997 

65 

1.972 

90 

1.947 

140 

1.897 

Appendix  IV  B 
Velocity  Distributions 
V^versus  r 


i^ft/sec 

0 

13.08 

16.93 

20.21 

22.77 

26.01 

28.19 

29.63 

30.04 

31.08 

31.37 


inches 

ft/sec 

1.6613 

(72.49) 

1.6693 

50.75 

1.6713 

47.92 

1.6763 

41.60 

1.6813 

38.30 

1.6863 

37.18 

1.6963 

35.  *7 

1.7063 

34.63 

1.7213 

33.90 

1.7463 

32.74 

1.8463 

31.44 

Rotor 

0 

8 

10 

15 

20 

25 

35 

^5 

60 

85 

135 
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Appendix  IV  B 
Velocity  Distributions 
j^-versus  r 


Run  48  V 

5000  rpm 

Trotor  = 

725.8  R 

^stator  = 

549.9  R 

4  T  - 

179.5  Jt 

Stator 

tyrn 

r,  inches 

l^ft/s< 

0 

2.037 

0 

5 

2.032 

11.30 

9 

2.028 

15.71 

12 

2.025 

18.99 

15 

2.022 

22.06 

20 

2.017 

25.34 

30 

2.007 

28.30 

40 

1.997 

29.62 

65 

1.972 

30.98 

90 

1.947 

31.86 

140 

1.897 

31.75 

Rotor 


r,  inches 

(^ft/si 

0 

1.6632 

72.49 

9 

1.6722 

52.29 

10 

1.6732 

51.18 

15 

1.6782 

44.07 

20 

1.6832 

40.64 

25 

1.6882 

38.66 

35 

1.6982 

36.85 

45 

1.7082 

35.81 

60 

1.7232 

34.62 

85 

1.7482 

33.30 

185 

1.8482 

32.74 
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Appendix  IV  B 
Velocity  Distributions 


Run  49  V 

5000  rpra 

Trotor  * 

880.5  R 

Tstator  ■ 

559-7  R 

A  T 

320.8  R 

Stator 

Jr,™  1 

?,  Inches 

ft/sec 

0 

2.037 

0 

5 

2,032 

11.48 

9 

2.028 

14.81 

12 

2.025 

18.08 

15 

2.022 

20.75 

20 

2.017 

24.22 

30 

2.007 

27.74 

35 

2.002 

28.67 

40 

1.997 

29.28 

65 

1.972 

31.28 

90 

1.9^7 

31.69 

140 

1.897 

32.21 

Rotor 

%/n 

r.  Inches 

Vg,  ft/sec 

0 

1 • 6656 

(72.7) 

9 

1.6746 

56.2 

10 

1.6756 

54.6 

15 

1 .6806 

47.1 

20 

1.6856 

43.1 

25 

I.6906 

4o.4 

35 

1.7006 

38.1 

45 

1.7106 

36.6 

60 

1.7256 

35.1 

85 

1.7506 

33.8 

185 

1.8506 

32.4 
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Appendix  IV  B 
Velocity  Distributions 
V#  versus  r 


Run  51  V 

7000  rpm 

Trotor  m  606.7  R 
^stator  -  5*0.2  H 
a  T  -  66.5  B 

Stator 

v>  in°hes 

j^ft/sec 

Rotor 

5t>” 

r.  Inches 

ft/sec 

0 

2.037 

0 

0 

1.6614 

(flOl.49) 

5 

2.032 

21.74 

8 

1.6694 

73.21 

9 

2.028 

29.71 

10 

1.6714 

66.33 

12 

2.025 

34.52 

15 

1.6764 

57.52 

15 

2.022 

37.09 

20 

1.6814 

54.03 

20 

2.017 

39.67 

25 

1.6864 

52.17 

30 

2.007 

41.95 

35 

1.6964 

50.13 

40 

1.997 

43.11 

45 

1.7064 

48.96 

65 

1.972 

43.90 

60 

1.7214 

47.75 

90 

1.947 

44.23 

85 

1.7464 

46.40 

140 

1.897 

44.62 

185 

1.8464 

44.81 
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Appendix  IV  B 
Velocity  Distributions 
\Jq  versus  r 


Run  50  V 
7000  rpm 

Adiabatic,  t  =»  80  P 


Stator 

if^ro 

r,  inches 

)^ft/sec 

0 

2.037 

0 

5 

2.032 

24.53 

9 

2.028 

32 .66 

12 

2.025 

36.29 

15 

2.022 

38.48 

20 

2.017 

40.75 

30 

2.007 

42.71 

40 

1.997 

43.87 

65 

1.972 

44.46 

90 

1.947 

44.66 

140 

1.897 

44.65 

Rotor 

r,  inches 

if*,  ft/se 

0 

1.660 

101.4 

8 

1.668 

73-53 

10 

1.670 

66.40 

15 

1.675 

57.26 

20 

1.680 

57.77 

25 

1.685 

51.88 

35 

1.695 

49.95 

45 

1.705 

48.74 

60 

1.720 

47.58 

85 

1.745 

46.58 

185 

1.8455 

45.34 
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Appendix  IV  B 
Velocity  Distributions 
^versus  r 


Run  53  V 

7000  rpm 

T  *s 

rotor 

888.7  R 

T  _ 

stator 

564.8  R 

d  T 

323.9  R 

Stator 

!frrr, 

r.  Inches 

l^ft/sec 

0 

2.037 

0 

4.5 

2.032 

17.36 

8.5 

2.028 

24.34 

11.5 

2.025 

28.52 

14.5 

2.022 

32.34 

19.5 

2.017 

36.14 

29.5 

2.007 

40.55 

39.5 

1.997 

42.55 

64.5 

1.972 

43.60 

89.5 

1.947 

44.51 

139.5 

1.897 

44.36 

Rotor 


t,  Inches 

ft/sec 

0 

1.6654 

0 

8 

1.6734 

80.99 

10 

1.6754 

74.49 

15 

1.6804 

63.24 

20 

1.6854 

57.50 

25 

1.6904 

54.35 

35 

1.7004 

50.98 

45 

1.7104 

49.53 

60 

1.7254 

48.23 

85 

1.7504 

46.46 

185 

1.8504 

44.78 

Appendix  IV  C 
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Appendix  IV  D 

1  2 
Friction  Factor,  f  and  Taylor  Number 


rpm 

A  T 

frotorxl° 

2  2 
^stator*10 

Taylor  Number 

3000 

0 

1.665 

1.105 

3829 

5000 

0 

1.430 

0.949 

6385 

7000 

0 

1.299 

0.870 

9092 

3000 

67.9 

1.735 

1.027 

3621 

5000 

79.8 

1.479 

0.872 

6035 

7000 

66.5 

1.332 

0.793 

8400 

3000 

179.5 

1.728 

0.916 

3018 

5000 

179-5 

1.589 

O.805 

5030 

3000 

292.5 

2.291 

2496 

5000 

320.8 

1.747 

0.748 

4003 

7000 

323.9 

1.85 

0.77 

5450 

1 

2 

-71*7 
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Appendix  IV  -  E 

(1  +  ^*40/f,  and  (  6»/p)/\ f  ,  versus  y+ 
3000  rpm,  Adiabatic 


f 

rotor 

=  I.665  x 

CXI 

1 

O 

r-l 

f 

stator 

=.  1.105  X 

CVJ 

1 

O 

rH 

Rotor 

y+ 

(1+  4^)/f 

( )/f 

4.8 

75.52 

15.45 

6.98 

86.63 

26 .56 

9.16 

112.37 

52.30 

11.34 

149.53 

89.46 

13.52 

200.58 

140.51 

15.7 

254.96 

194.89 

20.07 

406.55 

346.48 

24.43 

606.74 

546.67 

28.79 

719.3 

659.23 

33.15 

854.7 

794.63 

37.52 

1039 

978.9 

Stator 

y+ 

(  */v)/t 

4.83 

110.95 

20.43 

6.61 

122.28 

31.76 

8.39 

133.23 

42.71 

11.94 

176.06 

85.54 

15.5 

282.81 

192.29 

19.05 

447.14 

356.62 

22.61 

676.15 

585.63 

26.16 

1109 

1018.5 
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Appendix  IV  -  E 

(1  +  and  (  ,  versus  y+ 


3000  rpm 

T 

rotor 

-  606.1  R 

^stator 

-  538.2  R 

A  T 

«  67  *  9  R 

> 

rotor 

-  1.735  x 

10"2 

^stator 

-  1.027  x 

10"2 

Rotor 

Stator 

y+ 

(1+  a.mJ{  ( )/f 

y+ 

(1+  f 

3.85 

63.54 

5.91 

3-58 

109.77 

12.37 

7.7 

87.13 

29.5 

7.1 6 

128.12 

30.72 

11.55 

148.25 

90.62 

10.74 

161.93 

64.53 

14.25 

208.97 

152.34 

12.17 

182.15 

84.75 

18.1 

315.46 

257.83 

15.75 

269.32 

171.92 

21.95 

453.24 

395.61 

19.33 

433.3 

335.9 

25.8 

582.94 

525.31 

22.91 

768.1 

670.7 

29.6 

726.96 

669.33 

26.49 

1422.3 

1324.9 

33.5 

955.3 

897.67 

30.07 

2286.1 

2188.7 

33.65 

3096.5 

2999.1 
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Appendix  IV  -  E 

(1  +  and  ,  versus  y+ 

3000  rpm 


T  ^ 
rotor 

=  725.7  R 

T 

stator 

»  546.2  R 

A  T 

=  179-5  R 

f  t 

-  1.728  X 

10'2 

rotor 

f 

stator 

-  0.916  X 

10“2 

Rotor 

y+ 

(1+  S"/#)J 

t  ( )/t 

2.84 

58.74 

.87 

4.26 

60.36 

2.49 

5.68 

67.24 

9.37 

7.10 

74.74 

16.87 

8.52 

93.17 

35.30 

9.94 

122.33 

64.46 

10.79 

150.41 

92.54 

13.63 

237.38 

179.51 

16.47 

293.23 

235.36 

19.31 

383.81 

325.94 

22.15 

517.69 

459.8 

24.99 

599.92' 

542.05 

Stator 

y+ 

(1+ 

( 

3.27 

113.11 

3.93 

4.9 

123.8 

14.62 

6.53 

138.72 

29.54 

8.17 

162.41 

53.23 

9.8 

200.24 

91.06 

13.07 

333.61 

244.43 

15.68 

426.85 

317.67 

17.32 

474.24 

365.06 

20.59 

565.57 

456.39 

23.86 

647.73 

538.55 

27.13 

745.56 

636.38 

30.40 

795.62 

686.44 
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Appendix  IV  -  E 

(1  +  £nofy/Z ,  and  (  ^n/pi/Z,  versus  y+ 


3000  rpm 


947.9  A 
-55*.  4- 
292.  S/S. 

Rotor 

y+ 

(1+  4&)/t 

if&f 

2.49 

44.5 

0.85 

4.97 

52.6 

9.00 

7.46 

78.5 

34.9 

9.98 

122.6 

78.95 

12.43 

169.8 

126.2 

14.91 

218.2 

174.6 

17.40 

280.34 

236.7 

22.37 

447.73 

4o4.1 

886 


Appendix  IV  -  -E 

(1  +  *"A/)/tt  and  (  6>Jv)/t,  versus  y 
5000  rpra  adiabatic 


f 

rotor  0 

^stator  “ 

Rotor 

y+ 

1.430  x  10"2 
0.949  x  10"2 

(1 

(34/>/f 

Stator 

y+ 

{*p/t 

8.09 

105.85 

35.93 

9.35 

178.5 

73.1 

11.46 

137 

67.1 

12.1 

226.77 

121.4 

14.83 

185.5 

115.6 

14.85 

299.79 

194.4 

18.2 

265.1 

195.2 

17.6 

415.15 

309.8 

21.59 

396.2 

326.3 

20.35 

571.78 

466.4 

24.94 

595.2 

525.3 

23.I 

757.92 

652.5 

28.31 

835.1 

765.2 

28.6 

1116.5 

1011.1 

31.68 

1082.6 

1012.7 

34.1 

1482.44 

1377 

38.42 

1437.7 

1367.8 

39.6 

1927.34 

I821.9 

45.16 

1777.2 

1707.3 

45.1 

2893.28 

2787.9 

58.64 

2761.5 

2691.6 
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Appendix  IV  -  E 


(1  +  and  ( 


5000  rpm 

Trotor  “ 

606.2  R 

J?  s 

stator 

526.4  R 

A  T  - 

79.8  R 

f 

rotor  = 

1.479  x  10"2 

^stator  = 

O.872  x  10’2 

Rotor 

y+ 

( 

8.41 

103.9 

36.3 

11.37 

135.  s 

68.2 

17.28 

247.8 

180.2 

23.2 

478.4 

410.8 

29.12 

918.3 

850.7 

35.04 

1315.2 

1247.6 

40.96 

1748 

1680.4 

46.88 

2157.2 

2089.6 

52.81 

2563.8 

2496.2 

versus  y+ 


Stator 

y+ 

(l+6^)/f 

(V*/)/ 

5.27 

151.1 

36.4 

10.54 

220.7 

106. 

15.81 

341.1 

226.4 

21.08 

611.7 

497. 

26.35 

1152.3 

1037.6 

31.62 

1639.7 

1525. 

36.89 

2105.1 

1990.4 

42.16 

2622.2 

2507.8 

47.43 

3232.5 

3117.8 

52.7 

3953.9 

3839.2 
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Appendix  IV  -  E 

(1  +  ^y^)/ f,  and  (  £<oL>)/t ,  versus  y+ 


5000  rpra 

T  . 
rotor 

725.8  R 

Tgtator  = 

549.9  R 

d  T 

179.5  R 

f  — 

rotor 

I.589  x  10-2 
-2 

f  a 

stator 

0.805  x  10 

Rotor 

y+ 

(1+  ^)/f 

4.52 

76.4 

13.5 

9.04 

96.7 

33.8 

13.56 

168.6 

105.7 

18.08 

307.6 

244.7 

22.61 

515.1 

452.2 

27.13 

854 

791.1 

31.65 

1242 

1179.1 

36.17 

1513.4 

1450.5 

40.69 

1554.8 

1491.9 

Stator 


y+ 

(1+^  )  f 

(^/ 

5.07 

151.5 

27.3 

10.15 

193.4 

69.2 

15.22 

295.7 

171.5 

20.29 

487.5 

363.3 

25.37 

815.2 

691. 

30.44 

1273.2 

1149. 

35.51 

1719.5 

1595.3 

40.59 

2186.3 

2062.1 

45.66 

2537.9 

2413.7 

269 


Appendix  IV  -  E 

( 1  +  ^JJ)/ f ,  and  (  )/f ,  versus  y+ 


5000  rpm 

T  . 
rotor 

880.5  R 

T  m 

stator 

559.7  R 

it- 

320.8  R 

f 

rotor  - 

1.747  x  io~2 

^stator  = 

0.748  x  10"2 

Rotor 

Stator 

y+ 

(l+*%  )/f 

( )h 

y+ 

(1+ 

P  i+/*/t 

3.75 

64.4 

7.2 

5.12 

69.6 

12.4 

5.92 

139.9 

6.2 

6.49 

75.6 

18.4 

7.81 

181.5 

47.8 

9.22 

93.9 

36.7 

9.71 

202.5 

68.8 

11.96 

141.9 

84.7 

13.5 

266.6 

132.9 

14.69 

222.3 

165.1 

17.3 

333.9 

200.2 

17.42 

344.2 

287 

21.1 

465.3 

331.6 

20.15 

5H.7 

454.5 

24.9 

705.3 

571.6 

22.88 

618 

560.8 

28.6 

987.5 

853.8 

25.62 

659.6 

602.4 

32.4 

1277.9 

1144.2 

26.98 

723 

665.8 

36.2 

1464.8 

1331.1 

38.1 

1506.1 

1372.4 
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Appendix  IV  -  E 

(1  +  f,  and  (  €m/z>)/Z,  versus  y+ 

7000  rpm  Adiabatic 


f  m 

rotor 

1.299  x  10"2 

^stator  = 

0.87  x  10"2 

Rotor 

Stator 

y+ 

(1+ 

y+ 

U+<S£)/r 

12.52 

149.5 

72.5 

6.55 

170.1 

52.45 

17. 

233.5 

156.5 

10.22 

204.2 

86.6 

21.5 

393.2 

316.2 

13.87 

267.5 

149.9 

25.9 

617.8 

540.8 

17.52 

371.5 

253.9 

30.4 

894.2 

817.2 

21.17 

537.6 

420 

34.9 

1205.3 

1128.3 

24.82 

788.1 

670.5 

39.4 

1534.5 

1457.5 

28.47 

1115.7 

998.1 

48.3 

2216.6 

2139.6 

32.12 

1530 

1412.4 

57.3 

2777.3 

2700.3 

35.77 

1967 

1849.4 

66.2 

3281.8 

3204.8 

39.42 

2359 

2241.4 

75.1 

3742.3 

3665.3 

46.72 

2919 

2801.4 

54.02 

3963.8 

3846.2 

61.32 

6341.8 

6224.2 

68*62 

8796.3 

8678.7 
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Appendix  IV  -  E 

(1  +67/^)/ f,  and  (  GnJcJ)/ f,  versus  y+ 


7000  rpm 

^rotor  s 

606.7  R 

T  a 

stator 

540.2  R 

A  T 

66.5  R 

f  a 

rotor 

1.332  x  10"2 

^stator 

0.79  x  10~2 

Rotor 

y+ 

(l+^)/f 

( 

7.91 

114.8 

39.7 

13.^5 

163.7 

88.6 

21.36 

361.2 

286.1 

29.27 

775.2 

700.1 

37.18 

1362.6 

1287.5 

45.09 

1953.4 

1878.3 

53 

2422.8 

2347.7 

60.91 

2871.5 

2796.4 

68.82 

3277.9 

3202.8 

Stator 

y* 

(l+^)/f 

7.3 

179.7 

51.8 

14.6 

279.4 

151.5 

18.25 

388 

260.1 

24.1 

681.1 

553.2 

31.39 

1301.7 

1173.8 

38.69 

2092.1 

1964.2 

45.99 

3019 

2891 

53.29 

4186.6 

4058.7 

60.59 

5897.2 

5769 

67.89 

7779.6 

7651.7 

2*2 


Appendix  IV  -  E 

( 1  +  )/f ,  and  ( )/f ,  versus  y+ 

7000  rpra 


Stator 

« 


y+ 

(1+W#)/? 

(^)/f 

6.69 

175.29 

45.29 

13.37 

292.86 

162.86 

20.06 

463.77 

333.77 

26. 74 

778.70 

648.7 

33.43 

1228.21 

1098.2 

40.11 

1642.73 

1512.73 

46.80 

2242.2 

2112.2 

53.48 

3441.6 

3311.6 

56.83 

5120.4 

4990.4 

Rotor 


y+ 

(l+^4/)/f 

<  ®/W/4 

7.1 

82.3 

28.3 

9.56 

107.9 

53.9 

12.03 

143.8 

89.8 

14.49 

198.2 

144.2 

16.96 

279.7 

225.7 

19.42 

376.1 

322.1 

21.88 

491.4 

437.4 

26.81 

711. 

657. 

29.28 

868. 

814. 

34.2 

1153. 

1099. 

39.1 

1643. 

1589. 

44.1 

2018 

1964 
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Appendix  V 

Fluid  Properties  of  Air  Used  In  the  Computations 


For  evaluation  of  the  Integrals  In  Chapter  IX 
containing  as  factors  the  fluid  properties  ^  and  #  it 
was  desirable  to  express  both  these  quantities  as  functions 
of  the  absolute  temperature  T. 


Values  of  *  and obtained  from  Keenan  and 
Kaye,  Oas  Tables1,  Table  2,  Air  at  Low  Pressures,  p.  34 
were  plotted  on  log-log  coordinates  as  In  Fig.  V-l. 


If  constant  Prandtl  number  and  constant  specific 
tssumed, 
with  temperature. 


heat  are  assumed,  then  and i  vary  in  the  same  manner 


As  can  be  seen  from  the  figure,  the  curves  for 
and  JU-  are  slightly  divergent.  A  straight  line  was  drawn 
representing  the  average  slope  between  the  two  curves.  This 


line  varies  as  T 


0.75 


From  McAdams,  Heat  Transmission,2  2nd  edition,  the 
thermal  conductivity  of  air  at  212°F,  from  the  data  of  Moser, 
Dissertation,  Berlin  1913#  is  given  as  0.0183  Btu/hr-ft2-degF/ft. 
A  line  varying  as  T0*?^  was  passed  through  this  point  resulting 


In 


m  1.389  x  10'4  T0*75 


^•J.  H.  Keenan  and  J.  Kaye,  Qas  Tables  (1st  ed.,  John  Wiley 
and  Sons,  Inc.,  New  York,  New  York,  1956),  p.  3^« 

2 

McAdams,  Heat  Transmission,  p,  391* 


1  75 

The  kinematic  viscosity  Ovaries  with  T  as  T  under 
the  assumptions  made. 


of 


Using  a  value  ofyWcentlpolses  from  McAdams,  Pig.  204 


LL  =»  0.021  at/lOO  C  -  212  p) 

/  centipoises  ^  S 

a  relationship  for  is  obtained  as 
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